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VoL. 64, Part 2 1 February 1951 No. 374A 


The Theory of the Propagation in Liquid Helium II of 
‘Temperature-Waves’ of Finite Amplitude 


By H. N. V. 'TEMPERLEY 
King’s College, Cambridge 


Communicated by D. Shoenberg; MS. received 26th June 1950 


ABSTRACT. ‘The two-fluid model of liquid helium II is applied to the problem of 
predicting the behaviour of waves and pulses of second sound of finite amplitude. It is 
found that equations quite analogous to the Riemann and Rankine—Hugoniot equations for 
ordinary sound waves and shock fronts can be established. These equations predict the 
“overtaking effect’ with the establishment of a ‘ temperature-shock’ for pulses of second 
sound of finite amplitude. The temperature-shocks should be of ‘ front-edge’ type (propa- 
gation of an abrupt rise in temperature) if the liquid helium is at a temperature below 
1-96° k., but should become ‘ back-edge ’, corresponding to an abrupt fall in temperature, 
between 1-96° k. and the A-point. There is experimental evidence in support of both these 
predictions. It is shown that, just as for ordinary shock fronts, the irreversible effects due 
to the passage of a temperature shock should be proportional to the cube of the height of 
the shock. 


§1. INTRODUCTION 

ECENT measurements on the propagation of temperature differences in 
R liquid helium II (Pellam 1949, Osborne 1951) have provided evidence, 
first that a pulse of finite amplitude changes its shape as it travels 
forward, and secondly that the front of such a pulse may become effectively 
vertical when it has travelled far enough. This suggests that this phenomenon of 
“second sound’, the propagation of temperature differences without appreciable 
attenuation or dispersion over a wide range of frequencies, may have other 
properties similar to those of ordinary sound. In particular, one is tempted to 
look for something analogous to the Riemann theory of the propagation of sound of 
finite amplitude, and to the Rankine-Hugoniot theory of the propagation of 
pressure discontinuities or ‘shock fronts’. In this paper such results will be derived 
on the basis of the ‘two-fluid’ model of liquid helium I1—that is, on the postulate 
that liquid helium II can be represented as consisting of two mutually inter- 
penetrating liquids each having its own density and velocity-field and being capable 
of relative motion without appreciable friction, at all events as long as the relative 
velocity is not too great. ‘This model has been shown (Gorter and Mellink 1949) 
to give a reasonable quantitative account of many of the experimental facts, but its 
relationship with the ‘ surface creep’ phenomena that are known to occur with 

liquid helium II has still to be worked out. 
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§2. THE EQUATIONS OF MOTION 


Various attempts have been made to derive the equations of motion in such a 
way as to include second and higher order terms in the velocities (for example, 
Landau 1941, Zilsel 1949, Gorter and Mellink 1949). Unfortunately, the higher 
order terms obtained by different workers do not agree, and the matter is therefore 
re-investigated in the Appendix, the equations derived being closely similar to 
those of Zilsel (1949), but differing slightly from them. 

We obtain the following two equations, where p, and p, are respectively the 
masses of normal and superfluid components per unit volume, and u,, u, the 
corresponding hydrodynamic velocities: 


0 7) oP 
at (Prt n a7 pss) an ax (Prt ay psll,”) si Ox = 0, she, etenes (1) 
op , a p 0(pS) 
aa ar | (e+ £rs) Pit + 8p. | +1 
lo 2 2, 10 3 ae 
an } 5p (Palen + pslls ) me 2 ax (Prttn + plts ) =0. ...... (2) 
In addition we have the equation of continuity, 
Op a “a 
an ag ax (Pyln+ Pills) =O ee wine (3) 


and the postulate that entropy is transported by one of the fluids only, the normal 
component, the superfluid component being unable to carry entropy : 


ee + S (pu, 5) 23) eet Were Ux eae (4) 


Equation (4) also states that entropy is conserved, which means that the effect of 
irreversible processes is neglected. Such neglect seems justified by the experi- 
mental fact that both absorption and dispersion of second sound are small over a 
wide range of frequencies, except possibly in the neighbourhood of the transition 
temperature. 

We introduce the two velocities V and W, which we may call the ‘ stream 
velocities’, associated respectively with ordinary sound and second sound, defined 
as follows : 

pV =pat pals Wade eee (5) 
V is a measure of the translation of the liquid as a whole, and W is the relative 
velocity of the two components. In terms of these velocities, equations (1) to (4) 
may be rewritten 


S (pV) + jz (pv2+ Sf we) +2 =0, re (6) 
pSW +5 a (2) oe eae able 
= V5 oe | 0 re (7) 
_ £ (pv) =0, ae (8) 
+= (pVS+pSW)=0. ase, (9) 
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‘| _ In deriving equation (7) from equation (2) use has been made of equations (6) and 


(8) to eliminate'terms involving V and the pressure. Equations (6) and (8) show 
us that, if thermal expansion can be neglected, ordinary sound is excited by second 
sound only through the third term in equation (6), which is of order W2, so that, if 
we have a purely thermal source, the ordinary sound will have an amplitude 
proportional to W?, and we can consistently neglect the third term in equation (9), 
of order W®, and the terms involving VW? in equation (7), of order W4. If the 
effect of thermal expansion becomes important, as it may near the transition 


|| temperature, we can no longer neglect these terms with consistency, as the 


amplitude of V is then proportional to that of W, and the third term in equation (9) 
is then only of order W?, so that we cannot then consistently consider the non- 
linear theory of second sound while neglecting the first sound that it excites. 


§3. THE RIEMANN THEORY OF SECOND SOUND 
Assuming a zero coefficient of expansion and neglecting V, equations (7) and 
(9) may be written in the form 


OT | pnpsOW | 1Wps—pn)Opn | 3 py PnPs(Ps—Pn) OW _ Kars 
Ps oe ars Cet ere p eye Be eee Ot) eee eters (10) 
aw ae as. as 
pS a —-WSa- + We. TD Ap =(). sie witeiteie! (11) 


We now make an assumption, originally due to Tisza (1940), that the entropy per 
unit mass of normal fluid is constant, that is pS =p,S,, where S, is a constant. 
Taken together with equation (4), this would imply that the total masses of normal 
and superfluid are separately conserved, and thus neglects the possibility of 
transitions between the two fluids, so that it cannot be more than an approximation. 
However, our experimental knowledge of the behaviour of p, and S as functions of 
temperature indicate that the accuracy of the assumption is probably reasonable 
down to about 1°k., so that we put the consequences of the assumption on record, 
until the behaviour of S,, as a function of temperature is accurately known, when 
the terms involving S|, can easily be calculated. Making this assumption, and 
introducing the velocity C, of second sound, defined by the equation 


7 a Ps S? La PsSn 1 » 
oe PROS Oe Op OL ee I a ee) 
and a function o, analogous to the Riemann function, defined by 
PE Pole mea 4" L895 (13) 
Pn Ps 


we may write equations (10) and (11) in the form 


do = W(ps—pyn) OW poll 
Cra, + creer eyo or == (ieee delaeorce ct (14) 
OW W(ps— Pn) Gon Coa 15 
C255 at soir aiag ap =) eee eee a oh (15) 


which imply that the quantity o + W is propagated with a velocity C, + W(p,—p,)/p, 
while the quantity o— W is propagated in the opposite direction with a velocity 
C,—W(p,—p,)/p. (To see this, add and subtract equations (14) and (15).) 

We have thus arrived at a theory quite analogous to the Riemann theory of 
ordinary sound in one dimension, but there are some very important differences : 


The Riemann theory of first sound is exact, whereas we have neglected terms of the 
8-2 
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order W? in addition to those involving S, in deriving equations (14) and (15). 
There is no particular difficulty in taking these terms into account, but, as we have 
seen, we ought then logically to consider also the first sound excited by the second 
sound. Secondly, the presence of the term W(p,—p,)/p in the propagation 
velocity is, in general, much more important than the fact that C, itself varies with 
temperature, whereas, with ordinary sound, the particle velocity is usually small 
compared with the change of sound velocity with pressure. Thirdly, the theory 
predicts that the front edge of a temperature pulse should progressively become 
steeper only if p,>p,, that is, if T <1-96°x. For temperatures between this 
and the transition temperature, the theoretical prediction is that the steepening 
should occur at the rear edge of the pulse. However, the coefficient of expansion is 
becoming appreciable in this region of temperatures, and the excitation of first 
sound and the consequent cross terms can probably be no longer neglected. 


§4. SHOCK-FRONT THEORY 


We attempt to find a solution of equations (6) to (9), representing a disturbance 
of constant form travelling forwards with a constant velocity U, by imposing a 
backward velocity U upon everything and then setting all time derivatives equal to 
zero. Equation (8) then becomes (0/0x)[p(U— V)]=0 or p(U—YV) is a constant, 
which constant we may put equal to p)U, where the suffix zero always refers to 
conditions in the undisturbed liquid, so that V,) and W, are always zero. Thus we 
may write equation (6) in the form 


Py +p,U2=P+p(U—V)?+ a W?. 


Equations (8) and (9) may also be handled in this way, but equation (7) is not quite 
so simple, because, when we put all time derivatives equal to zero, the remaining 
terms do not form an exact differential in the variable x. However, equation (2), 
from which equation (7) was derived, does have this property, and so we may use it 
instead. Equations (7) to (9) are thus replaced by 
(So + ToSo)poU + 2poU? =glpn( U — uy) + p(U —u,)] + TpS(U — uy) 
+4p,(U—u,)tip(U— a) (17) 
pol =p(U SV)" ee (18) 
poUS) =pS(U— V)—p.SW. 
We now suppose that, just as in ordinary shock-front theory, these equations 
can be applied when the differences of temperature and pressure are finite but 
steady, so that the only effect of viscosity, thermal conduction and mutual friction 
between the fluids is to prevent the shock front from becoming completely abrupt. 
Since, in our model, areversible mechanism exists by which temperature differences 
can give rise to relative motion of the two fluids, there may not always be the 
irreversible transformation of potential or kinetic energy into heat that occurs at 
shocks in other fluids. The interesting question therefore arises whether a 
sudden drop in temperature can ever be propagated asashock. A negative pressure 
shock is excluded because the compressibility decreases and, therefore, the 
velocity of ordinary sound increases with pressure. The Riemann theory shows us 
that, above 1:96°K., the front edge of a positive temperature pulse should be 
spreading out and the back edge sharpening up, so that this region of temperature 
seems to be a possible one for ‘back-edge’ or negative shocks, and this will be 
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confirmed by a calculation of possible shock-front velocities. This effect was 
confirmed experimentally by Osborne (1951), but there are disturbing effects, 
such as the thermal expansion and the possible attenuation of second sound in 
this region of temperature (Pellam 1949), which may prevent a quantitative 
verification. The exact causes of this attenuation are not known, but one of them 
may well be a conversion of second sound into ordinary sound by one of the two 
mechanisms already mentioned, both of which should theoretically become 
more and more effective as we approach the A-point. 

From equation (17), written in terms of the velocities V and W, by expanding 
&—& in ascending powers of P— Py and T— To, using equation (16) to eliminate 
the term in P— P), and neglecting the thermal expansion, we get 
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From equations (18) and (19) we get 
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In the absence of any source of first sound, equation (16) shows that P— P, and 
therefore V and p—py are all of the order W?. 
If in equation (20) we neglect all terms of order W4, we obtain for the velocity of 
a ‘temperature shock’ the possible values 
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For a front-edge shock to maintain itself, its velocity must be numerically greater 
than C,, the velocity of propagation of a disturbance of very small amplitude, so 
that we are led to expect front-edge shocks for p, greater than p,, but otherwise 
back-edge shocks. The velocity of propagation for a shock of height T— Ty is thus 
equal to the Riemann velocity associated with a wave of height $(T— T)). 

In deriving the velocity of propagation of a shock we have not had to use Tisza’s 
assumption. Indeed, by a process of successive approximation, we could derive, 
from the above equations, the velocity of propagation of a shock to any desired 
degree of accuracy, even if we did not neglect thermal expansion. 

As a check on these equations we may notice that if we put wu, =u, or W=0 
in equations (16) to(18) they become precisely the same as the shock-front equations 
in an ordinary fluid, equation (18) being the equation of continuity, equation (16) 
that of conservation of momentum, and equation (17), after the use of equation (18), 
becoming just the conservation of enthalpy plus kinetic energy. Equation (19) 
is inapplicable, except for very weak shocks, in an ordinary fluid, and for a shock to 
be physically possible there must be an increase of entropy in the fluid after the 
shock has passed. In this way it is possible to rule out the possibility of negative 
pressure shocks (subject to reasonable assumptions about the equation of state). 

Equation (19) when applied to liquid helium apparently indicates that any 
‘spare ’ heat can be taken up by a change in W, the relative velocity of the two 
fluids, instead of being dissipated as in an ordinary fluid. ‘The experimental fact 
that second sound can be propagated without appreciable attenuation if the pulses 
are not steeply sloped indicates that this thermodynamically reversible mechanism 
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predominates for waves and gently sloped pulses, and one would expect the 
minimum possible thickness of a shock front to be determined by the condition that 
the rate of transfer of energy by conduction or viscous damping becomes 
comparable with the rate of transfer by the reversible conversion of temperature 
difference into relative motion of the fluids. 

Consider a pulse of height AT and maximum mutual velocity W, and suppose 
that it rises to its peak in a distance Ax. The rate of transfer of thermal energy per 
unit area is of the order of magnitude CypSAT ~pC,?W. If heat conductivity were 
present alone, the rate of transfer of heat across the shock front would be of the 
order of magnitude KAT/Ax per unit area. Taking K as of the same order as for 
helium I, that is 6 x 10-5 cal/deg/cm/sec., p as 0-15 gm/cm3, C, as 2 x 10? cm/sec. 
and S as 10’ ergs/gm/degree, we obtain Ax~10-*cm. If viscosity were present 
alone, the rate of loss of kinetic energy would be of the order of magnitude 
1WC,/Ax per unit area. Taking 7 =30 x 10-* poise, we obtain as an estimate 
for the thickness of the shock front Ax~10-7cm. ‘These estimates are of the 
same order of magnitude as those for the thickness of an ordinary shock. 

Another way in which we might expect a lower limit to be set on the thickness of 
a shock front is through the mutual friction between the two liquids (proportional to 
W*) inferred by Gorter and Mellink (1949) from flow experiments in slits and tubes. 
By itself this mechanism does not seem to constitute a dissipative process, even if, 
as is not at all certain, the W? form applicable to processes of steady flow can be 
applied to second sound. For mutual friction between the fluids indeed provides a 
mechanism by which kinetic energy can be converted into heat, but since this 
heat will nearly all appear in the places where the amplitude of the mutual motion is 
a maximum, it can again be converted into kinetic energy unless heat conductivity 
or the viscosity of the normal fluid come into play. It therefore seems that this 
mutual friction does not, by itself, limit the thickness of a shock. 

There is one point which needs discussion in connection with the existence of 
shock fronts associated with asudden dropintemperature. At first sight they might 
appear to contradict thermodynamics, in that the passage of such a shock leaves the 
liquid behind it at a reduced temperature and, therefore, with a smaller entropy 
(since the density is practically unaltered). We consider the conditions that 
would have to be satisfied by a source capable of maintaining such a shock at a 
steady amplitude as it travels down a tube. (This is analogous to the problem of 
maintaining a steady pressure shock in a gas by means of the motion of a piston 
behind it.) If the shock is travelling in the positive direction with a velocity U, 
and the source is supposed to be at the origin, uw, will be in the negative direction 
while u, will be in the positive direction. If these conditions are to persist right 
up to the source, it can only mean that there is a continual conversion of normal 
fluid into superfluid at the source, that is to say, that the source must contain some 
mechanism for withdrawing entropy from the liquid. (It would be possible to do 
this, for instance, by means of a paramagnetic salt in a magnetic ficld that was 
steadily diminished with time.) The necessary rate of withdrawal of entropy 
would be p,uv,,S,, which, by equation (19), exactly balances the rate of change of 
entropy content of the system U(p»S,)—pS) due to the advance of the shock. 
However, even in ‘ideal’ circumstances, this balance between the rate of with- 
drawal of entropy at the source and the rate of loss of entropy by the liquid holds 
only in the limit of very small shocks. For shocks of finite amplitude we have also 
to take account of the changes in kinetic energy that must occur at thesource. The 
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rate at which kinetic energy is brought up to the source by the normal fluid is 
—p,u,°, while it is being withdrawn from the source by the superfluid at the rate 
$pu,°. If no ordinary sound is present p,u,+p., vanishes. The rate at 
which kinetic energy is having to be created is thus proportional to u,(u,’— 1,2), 
which is positive if |u,| > |u,|, that is, if p,.<p,. In these conditions we have the 
following state of affairs. Work has to be done at the source (e.g. by pulling the 
magnetized salt into a region of weaker field) in order to withdraw entropy at the 
rate p,u,S,. Work hasalso to be done to create the extra kinetic energy. Thus, if 
Ps< Pn the work done at the source is greater than that corresponding to the entropy 
removed from the liquid even if nothing but reversible processes occur at the 
source, and the regime is therefore thermodynamically possible, whereas a 
negative shock is not thermodynamically possible for p, greater than p,. An 
exactly similar argument shows that steady positive shocks are impossible for 
p; less than p,. ‘The ‘ dissipative’ terms which arise through the existence of 
the two velocities wu, and u, are thus of the order (AT)? or W?, just as the dis- 
sipation in pressure shocks is of the order (AP)? or V. 

In practice, steady negative shocks have not yet been observed, but positive 
pulses of temperature have been produced, and negative shocks have appeared at 
their trailing edges when p,<p, (Osborne 1950). It seems that we could 
calculate their attenuation by an argument exactly analogous to that used to calculate 
the dissipation of energy in explosion pulses, where we have a shock followed by a 
Riemann wave of decreasing amplitude. ‘he Riemann wave travels faster than 
the shock and is ‘eaten up’ by it, the energy in the portion of the wave thus destroyed 
corresponding just to the energy dissipated by viscosity and thermal conduction 
in the limiting case when these are very small but the shock front very steep. In 
our case the shock at the trailing edge travels faster than the Riemann wave (at a 
velocity of approximately C,—4(p,—p,)W compared with C,—(p,—p,)W, and 
the latter is again gradually ‘eaten up’. Physically this means that the heat 
originally contained in the positive pulse is eventually dissipated into the whole 
region traversed by it. When p, >p, dissipation occurs via a front-edge shock, the 
analogy with an ordinary pressure pulse then being almost exact. 

The assumption implicit in equation (19), that the rate of creation of entropy at 
the shock front is negligible, is thus seen to lead to thermodynamically consistent 
results. This arises because of the existence of two velocity fields. 
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APPENDIX 


Derivation of Equations (1) and (2) of the Text 


We apply the principles of conservation of energy and momentum to a small 
volume fixed in space, enclosed between parallel planes at distances x and x+ dx 
from the origin, of unit cross section. The extension to three dimensions is 
obvious. Applying first the principle of conservation of momentum, the net force 
onthe volume is — (0 P/dx)dx, which must be equal to the total change of momentum 
within the volume. This change arises from the changes in both density and 
velocity of both normal and superfluid within the volume, together with the net 
momentum carried into the volume by the motion of the two fluids. Since the 
rate of transport of momentum by a fluid of velocity wu and density p is pu? per unit 
area, we have 


7) 0 oP 
at (Pry i pels) ils ax (ppt, 1 pelts”) ae Ox =0, seeeee (Al) 


which is equation (1) of the text. 

We apply the principle of conservation of energy inthe following form. Since 
we are supposing that no heat enters the volume by conduction or radiation, the 
sum of the rates of change of internal energy and kinetic energy due to the changes. 
of densities and velocities within the volume and the net rates at which these 
energies are being brought into the volume by convection must be zero, that is 


Opn) (BE) Sn, (8ps\ (9E a OE OE 
BE NON Js (Lt ot | NOE AON a 0 Ox lett? NOI) ie selene 


10 pe ia ; 
+ 5 55 (Pata + sts) + 5 5 (Prtln® + Pats’) =. eae 


The first and third terms in equation (A2) represent the respective changes in 
internal energy due to variations in the masses of normal and superfluid within the 
volume, while the second term allows for a possible specific heat of the normal 
- fluid, S,, being the entropy per unit mass of the normal fluid, which may be a 
function of pressure and temperature. We assume that the entropy of the super- 
fluid is always zero, so that there is no corresponding term for the superfluid. 
(dE/ON), and (@E/dN), represent respectively the changes in internal energy due to 
the transfer at constant volume of unit mass of normal fluid and superfluid 
respectively. 

These quantities can be derived thermodynamically as follows. For a homo- 
geneous one-component assembly 


dh gdN¢ Td = PaV> =) ee (A3) 


where & is the entropy of the whole assembly and g is the Gibbs function per unit 
mass of liquid. Assuming that a transfer of superfluid involves a strictly zero 
change of entropy, we have for transfers at constant volume and entropy 
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while for transfers at constant temperature and volume 
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_ In order to maintain the temperature constant it is necessary to maintain within 
the volume the same value of p,/p,, which quantity is a function of the temperature 
only provided that we neglect thermal expansion. Thus we have 
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Gorter (1949) has called attention to the fact that a definite assumption is 
involved, for instance in replacing (0X/0N), by S, the entropy of unit mass of the 
whole liquid. ‘This assumption is not nearly so drastic as that made by Tisza 
(1940), in which S is assumed to be always proportional to p,/p, so that the quantity 
S,,=pS/p, is constant. Making Gorter’s assumption, and substituting from 
equations (A6) and (A7) into (A2), we obtain equation (2) inthe text, the term 
p,/ (dS, /0t) already mentioned being a correction for the fact that S,, may not be’ 
entirely independent of p,, in which case 0X/0N, taken at constant volume, is not 
quite equal to S. 

The equations of motion may also be written, after using the equations of 
continuity, 
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For comparison the equations obtained by Zilsel (1949) by a variation method are 
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If S,, is constant these become equivalent. 

Landau (1941), in his equation of motion for the superfluid, obtains an extra 
term described as the “ kinetic energy of mutual motion ”’, which is, apparently, 
based on the idea that the Gibbs function has to be corrected if the two fluids are 
in relative motion. The investigation given above seems to show, however, 
that the kinetic energy of the fluids is already allowed for through the Bernoulli 
terms in equations such as (A8) and (A9), and that the Landau correction is not 
needed. In any case, the introduction of the kinetic energy of mutual motion 
would only be correct in a frame of reference at rest relative to one of the fluids, and 
the remaining equations would have to be referred to a similar frame. 

It is quite easy to modify the Riemann and shock-front theories to take 
account of such a term, if it really exists. 


Note added in proof. While this paper was in the press, Zilsel published a revised 
derivation of the equations of motion using his variation method (Phys. Rev., 1950, 79, 309). 
He now obtains the term 4x(v,—v;)”, as was found by Landau, but otherwise his equations 
are the same as (A10) and (A11) above. In his notation x=py/p and U is the internal 
energy, while I measures the rate of transition between normal and superfluid, and the 
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terms involving S, in equations (A8)—(A11) above can also be expressed in terms of I. 
His present paper is open to the following criticisms : 

(a) It is not mathematically consistent. Zilsel’s equation (2.12), dUjox=H(Un—Us)*, 
implies the presence of a term }x(v,—vs)? in U, yet, in his equations (2.8) and (2.9), no terms 
arising from @U/év, or 0U/évs appear. 

(b) The term—div {3(pspn/p)v2v} occurring in his energy equation (3.12) is difficult to 
interpret physically. The kinetic energies of the normal and superfluid components are 
already taken fully into account by the other terms in this equation. 

(c) Zilsel introduces the very interesting idea that the ‘mutual friction’ introduced by 
Gorter and Mellink (1949) is a manifestation of the transition terms involving I or Sy 
together with the creation of entropy by the viscous mo.ion of the normal component. 
Gorter and Mellink had to postulate a mutual force acting on both normal and superfluid, 
in order to account both for the departure of the fountain effect from the theoretical value 
and for the linear relationship between fountain effect and heat conductivity. Such an 
apparent mutual force can be obtained from my equations, in which I appears in both 
equations of motion, but it is hard to see how Zilsel’s equations, where F appears only in 
the equation of motion of motion of the normal fluid, can give it. 
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ABSTRACT. A pulse technique for the propagation of second sound is described and 
the expected shape of the received pulse is discussed. In this connection, experimental 
evidence has been obtained for a temperature discontinuity at the interface between liquid 
helium II and a solid medium at second sound frequencies. ‘The pulse method has been 
used to study the attenuation of second sound with distance and the formation of shock 
waves. ‘The observations on the latter are shown to be in agreement with the predictions 
of the equations of motion of liquid helium II if second order terms are not neglected. 


§1. INTRODUCTION 


T is now well known that temperature variations in liquid helium II are 
propagated according to a wave equation, and these temperature waves 
have received the name of second sound. The two principal quantities 

characteristic of any wave motion, namely the velocity and the impedance, have 
already been measured for second sound and their values may be considered 
well known at temperatures between 1-1° k. and the A-point (2:19°k.), (Peshkov 
1944, 1946, Lane, Fairbank and Fairbank 1947, Pellam 1949, Osborne 1948). 
This being so, interest naturally begins to centre on secondary effects, and this 
paper describes measurements on the attenuation of second sound with distance 
and on the propagation of waves of large amplitude. The pulse technique, 
briefly mentioned in an earlier note (Osborne 1948) and developed also by 
Pellam (1948, 1949), has been used throughout and a short description of the 
method will first be given, followed by a discussion of the principles from which 
it is possible to predict the general shape of the received pulse. 
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§ 2. GENERAL DESCRIPTION OF THE TECHNIQUE 
A general scheme is shown in Figure 1. The essential features are a heater, 
usually of fine wire, from which the second sound is radiated, some kind of tube 
filled with helium down which it is propagated, and a resistance thermometer 
at which it is received and converted into an electrical signal. The auxiliary 
€quipment must include a means of supplying power to the heater and a means 
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Figure 1. General scheme of pulse technique. 


of amplifying and displaying the received electrical signal. In the case of the 
pulse technique, the requirements are very similar to those of a radar installation 
except that the pulse supplied to the heater is a pulse of direct current, not of 
alternating current, and it is the directly transmitted signal in which we are 
primarily interested, although echoes can readily be detected and do provide 
additional information. 

At the heart of the pulse circuits is a multivibrator which has a natural 
recurrence frequency of 0-5 to 50 cycles per second. ‘This device is loosely 
synchronized with a crystal oscillator which provides the standard of time in 
the form of calibrator pulses at intervals of 20sec. ‘The first duty of the multi- 
vibrator is to initiate a pulse for the heater, and this is accomplished by means 
-of a conventional circuit incorporating a variable delay so that the pulse may 
occur at some convenient time shortly after the beginning of the multivibrator 
“waveform. 

The next requirement is that a time-base be started at the beginning of each 
cycle, and for this purpose the multivibrator drives a Miller time-base which is 
connected, through an electronic switch whose function will appear below, 
to the X plates of the cathode-ray tube. 

There is also set in motion at the beginning of each cycle a mechanism for 
producing a short length of expanded time-base later on in the cycle. ‘This is 
managed via a variable delay which provides a ‘ strobe’ pulse after any desired 
time up to about 40msec. This strobe drives a fast time-base and it also 
operates the electronic switch so that, during the period of the strobe pulse, 
the main time-base is replaced by the fast one. ‘The strobe pulse fulfils one 
other function, namely the provision of a deflection voltage on the Y plate of 
the cathode-ray tube so that the fast trace shall appear below the main one. 

For timing purposes the calibrator pulses from the crystal oscillator are 
supplied to the grid of the cathode-ray tube so that bright spots are seen on the 
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trace at intervals of 1 msec., together with fainter ones at intervals of 4 msec., 
100 psec., and 20 psec. 

The phosphor bronze resistance thermometer which is used for detecting the 
signals carries a direct current of the order of a few milliamperes and the variations. 
of resistance consequent upon the arrival of a second sound signal at the ther- 
mometer cause corresponding variations of the voltage across it. These variations. — 
are amplified by a three-valve amplifier of conventional design with an input trans- 
former having a step-up ratio of 15 : 1, so that the overall voltage gain, including 
a small amplifier built into the oscilloscope, is about 10’, at frequencies from 
about 300 c/s. to 300 ke/s. 

The temperature of the helium bath is obtained from the vapour pressure 
of the helium, using the 1949 scale (van Dijk and Shoenberg 1949). 


§3. DISCUSSION OF THE EXPECTED SHAPE OF THE RECEIVED 
PULSE 

Several authors (Dingle 1948, Osborne 1948, Pellam and Scott 1949) have 
commented on the analogy between the propagation of second sound and of 
electrical waves on wires, and Peshkov (1948a) has given a formally similar 
analysis although without actual reference to electrical waves. ‘The analogy is a 
useful one, provided that all the relevant circuit elements are taken into account, 
but there is one important factor that has not been mentioned in previous treat- 
ments. ‘This is the observation of Kapitza (1941) that heat transfer from a solid 
surface into liquid helium II is accompanied by a temperature discontinuity at 
the surface proportional to the heat flow per unit area. The constant of pro- 
portionality K is of the order of a few degrees per watt cm~? and decreases towards 
the lambda point. Since in the analogy heat flow corresponds to current and 
temperature corresponds to voltage, K will appear as a resistance. 

The liquid helium itself is equivalent to a transmission line of characteristic 
impedance y given by 


y= (pCa). tia Wa as 6 oe a Gere (1) 


where pce is the specific heat of helium II per unit volume, and w, the velocity of 
secondsound. ‘This fact, which has been verified experimentally (Osborne 1948), 
may be deduced from the equations of motion which lead to the wave equation 
for second sound, but, as will now be shown, it is not in fact dependent upon this 
detailed theory. Let the assumption be made (say as an experimental fact) that 
in liquid helium II the propagation of temperature waves is possible without 
appreciable attenuation or dispersion. ‘Then consider in one dimension the 
progress of a sudden step of temperature which in such a medium will travel 
without change of shape at a velocity u,. Such a step will sweep out a volume 
Us per unit time per unit area of wave-front and within this volume all the material 
will rise in temperature by the amplitude AT of the step. The quantity of heat 
required to bring this about will be pcu,AT. With the step, therefore, there 
must be associated a flow of heat per unit area W in the direction of propagation, 
where W=pcu,AT. Therefore y=AT/W=(pcu,)1, as stated above. 

The behaviour of the solid media constituting the apparatus is to be obtained 
from the well-known equation of heat conduction 
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where 6(x, y, 2, t) is the temperature at any point at time ¢, s the specific heat per 
. t be Wo . . . 

unit volume, and k the thermal conductivity. For a semi-infinite solid undergoing 

simple harmonic temperature variation at its plane surface, the solution is 


0 =0, exp [ — (sw/2k)/2x] cos {wt — x(sw/2k)?}, 


where 0, cos wt is the temperature at the surface x=0. This is a wave attenuated 
to e+ of its original amplitude in a distance 
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‘The input impedance at the surface (i.e. temperature amplitude divided by the 
heat flow amplitude) is given by 
Zee escent ae OL! rahe (4) 

The representation of the propagation tube of Figure 2(a) by the equivalent 
circuit of Figure 2(b) can now be justified. The input power appears as a 
current J supplied to a capacity C, representing the heater wire. Since the 
penetration depth for the heater wire is large compared with its radius, this is a 
sufficient approximation. Following Kapitza’s observation, the heater is shown 
coupled to the helium through a resistance K,. The heater is usually stretched 
across a solid backing, but is not in good thermal contact with it, and the backing 
is therefore shown coupled to the helium through a separate resistance K,. 
‘The impedance Z, of the backing is to be calculated from equation (4). The 
column of helium, in so far as wall effects may be neglected, is imagined as a. 
transmission line of characteristic impedance y. At the thermometer end we 
have circuit elements similar to those at the heater end, and the temperature 
amplitude which is finally detected is represented by the voltage across C, the 
thermal capacity of the thermometer. 

The conditions that a rectangular pulse of heat of duration 7 supplied to the 
heater shall be radiated and received without appreciable distortion are 


|K,+2Z,| and |K,+2Z,|>y 
and C\(K,+y) and Ci\(K,+y)<r. 
At 1-7°K., with a heater and thermometer each consisting of 10 cm. of wire 
80 in diameter stretched (zig-zag) across a tube of 1 cm? cross section, 
y~10°3 deg. watt 7 
K,~K,~3 deg. watt 
Z,~Z~(1 —1)10 deg.watt™ at 1 ke/s. 
K,~K,~10 deg. watt 
: Cy. ~Ci~10 7 joule deg™?. 
‘Therefore |K,+Z,|~|K,+2Z,|~ 1-6 x 10*y 
and Ci (Ky +y)~C (K+ y)~1 psec. 


Since the pulse lengths commonly used are in the range 0-1 msec. <7 <1 msec., 
the above conditions for the absence of distortion are amply satisfied in practice. 
This fact has already been implicitly confirmed by the agreement between the 
measured and predicted values of y (Osborne 1948). It will therefore be 
assumed in $§$5 and 6 that the thermometer waveform is a true representation of 
the temperature variations being propagated through the liquid. With a 
thermometer mounted in front of a backing as illustrated in Figure 2(a), the 
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actual variation in temperature of the thermometer will be just twice the incident 
. amplitude in the helium, i.e. 2/y, because of the reflection at the backing. 

The above treatment has ignored attenuation in the helium and losses at the 
‘walls. The latter will be of the same order as the losses at the ends, i.e. negligible 
in the case considered. The former will be shown to be negligible except very 
near the lambda point (see Figure 5). 

To show that circumstances can arise in which the thermometer does 
introduce distortion, a form of apparatus having a heater and thermometer of 
graphite deposited upon a backing (Fairbank and Lane 1947) is illustrated in 
Figure 3, together with the equivalent circuit. This design has been used by 
several workers, and the present author has employed it in circumstances where 
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Figure 2. (a) Propagation tube with wire Figure 3. (a) Propagation Bike with 
heater and thermometer. (6) Equi- ‘ Aquadag’ heater and thermo- 
valent circuit. meter. (5) Equivalent circuit. 
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pulse length 150 usec. 
Circles : Experimental points (Osborne), 
pulse length 600 psec. 


the shape of the received pulse was unimportant.. The essential difference 
between it and the previous design is that the heater and the backing are now 
coupled to the helium through the common impedance K,. The current actually 
supplied to the line will therefore depend on the ratio of K, to Z, and, whilst 
the thermal data is very meagre, it suggests that K, and | Z,| are of the same order 
of magnitude. Similarly the pulse received at the thermometer will depend 
on the relative values of K, and Z,. It follows, since Z, and Z, are complex, 
that not only the amplitude but also the shape of the received pulse will depend 
on the relative values of K and Z at the temperature in question, and that if K is 
not much smaller than |Z], the received pulse will be more or less severely 
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integrated. ‘This in fact has been observed experimentally, the times of rise of 
the received puise being of the order of 100-200 psec. 

With a view to establishing the relevance of the impedance K to experiments 
on second sound, a propagation tube of the type shown in Figure 4 has been used. 
The first thermometer was an open wound one with a ‘shadow area’ of only 
10° of the cross section of thetube. Thecopper foil had a thickness of 0-0025 cm., 
to be compared with a penetration depth 6 of the order of 1 cm. In the absence 
of a boundary impedance, therefore, the foil should have presented no appreciable 
obstacle to the transmission of second sound and there should have been no 
reflection fromit. In fact there was no observable transmission and the reflection 
coefficient was very close to unity. A rough calculation on the basis of this 
result, taking into account the amplitude of the smallest transmitted signal which 
could have been detected, gives a lower limit of 0-1 deg.watt cm? for K at 1-2°x. 
This value is not inconsistent with Kapitza’s measurements using a steady heat 
flow (K~1 deg.watt-1 cm?). 


§4. SECOND SOUND WAVES OF FINITE AMPLITUDE 

It was discovered fairly early in the investigation that, whereas a small pulse 
of second sound was propagated without change of shape, larger pulses became 
seriously distorted. Plate I shows the effect in question. It represents a series 
of oscillograms of the received temperature pulse for increasing amplitudes of 
the input pulse supplied to the heater. The upper trace is the main time-base 
going from left to right and the lower trace presents an expanded version of 
what occurs during the gap in the upper trace. A negative deflection on either 
trace represents a positive temperature excursion. In Plate I(a)* the signal at 
small amplitudes is seen to have a flat top and sloping back and front edges. 
As the input is increased, the front edge becomes much steeper and actually 
occupies a shorter time, whereas the flat top and the back edge coalesce to forma 
long tail. At the same time the front edge not only becomes steeper but also 
occurs earlier. 

Plate II shows that similar phenomena are observed at constant amplitude 
if the thermometer is progressively withdrawn from the heater, provided that 
the initial pulse be large enough. Near the heater the pulse is a fair reproduction 
of the input, whereas at greater distances the characteristic sharp rise and long 
tail appear. 

This behaviour is familiar to students of ordinary sound as the generation 
of shock waves, and it is due to the non-linearity of the equations of motion of a 
compressible fluid. In the case of ordinary sound there are two contributory 
effects. In the first place the velocity of propagation increases with pressure 
so that, taking account of second order terms, the compressions will be propagated 
faster than the rarefactions; and in the second place any disturbance tends to be 
propagated at the velocity of sound relative to the local stream velocity, so that 
those points of the wave front where the stream velocity is in the same direction 
as the propagation will overtake the places where the stream velocity runs counter 
to the direction of propagation. Both of these effects are second order ones 
and are rarely observed for ordinary sound except in the very intense waves 
generated by explosions. Both effects cause the points of high pressure to 
overtake those of low pressure so that, for a positive pressure excursion, the 
leading edge becomes steeper and the trailing edge becomes a long tail. The 


* For Plates see end of issue. 
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process is limited when the leading edge becomes very steep by the dissipation 
due to the action of viscosity and thermal conductivity in the very steep pressure 
and temperature gradient at the wave front. This limiting condition is known 
as a shock wave or shock front, and the antecedent condition of a wave which is — 
becoming progressively more distorted is referred to as the propagation of a 
wave of finite amplitude. 

It is possible to obtain some idea of the physical meaning of the processes 
occurring in second sound if it is assumed that 


PoVn HP Ve=U 0 2 ee eee (5) 
and that the equation of motion is to be written down including the necessary 
second order term in the acceleration. ‘The conventional 


OVs of 
at vial G= 0! 2" a ee ee (6) 
OvV~. 

then becomes rea (v,. grad) v,+ grad G=0, 


where G is the thermodynamic potential per unit mass of helium. A con- 
ventional Riemann treatment for plane waves leads to the expression 

UP Stlgch Os.) oo T= eee (7) 
where w’ is the velocity of propagation of any plane in the wave front for which 
the superfluid velocity is v,, and u,, the velocity of propagation for very small 
amplitudes, is given by SESS 

uo? = —— —. 

c Pn 

c is the specific heat of helium per unit mass, S the entropy per unit mass, 
p, and p, the superfluid and normal parts of the density respectively, and T the 
absolute temperature. 

As far as second sound is concerned, however, and probably in some other 
applications, the motion of the normal fluid is as reversible as that of the superfluid 
component and a gradient of G (which may be written (1/p) grad p—S grad T) 
acts as a force upon both components. Equation (6), therefore, must be replaced 
by a form more symmetrical in v, and v,, and a convenient starting point is found 
in the formulation by Gorter and Mellink (1949), originally used in connection 
with their theory of mutual friction between the two components. If the friction 
and viscosity terms are omitted, there remain ‘ 


—Pigrad p+ WPS, grad T =p, 
p aie 
~<a lt Gl eae (9) 
_ Pn _ Pos ar 
grad p 3 S, grad T=p, at | 


where S, is the value of S at the A-point. As these equations stand, the use of 
equation (5) will lead straight to equation (6), but, if the second order terms 
(v,.grad)v, and (v,,.grad)v, are included, an equation different from (6) is 
obtained. If the Riemann derivation is now carried through we arrive at 


U = Uy bono bo te ees (10) 
This is the result given by Temperley (1951) in his rigorous derivation in which 


he proves that equation (5) is in fact sufficiently accurate for the purpose in hand 
provided that thermal expansion is negligible. 
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The first experimental results were obtained at a temperature OP e2Stk 
where p,>p,, and therefore, by equation (5), ¥,>v,. In this region, wu’ ~ uy + on 
and it should be possible to verify this formula by studying the rate at which 
the top of the leading edge of a pulse overtakes the foot. v, can be calculated 
from the measured power input W using the relation . 


ees Taree 1) Oe ge ae (11) 


Data taken from Plate IT give a value of wu’ —u, equal to + 140 cm. sec™, whereas 
v, calculated from W is about +155cm.sec™t. In view of the uncertainty of 
the thermal data and the difficulty of measuring the rates of rise of the pulse, 
this agreement is satisfactory and confirms the theory quite well. It shows 
- clearly the inadequacy of equation (7) which would predict a value of u’ — u, equal 
to —5cm.sec7. 

A further test is provided by the predictions at temperatures near the )-point 
for which p,>p,. For a positive temperature excursion 2, is always in the opposite 
direction to u, so that, for p,>p,,¥, + v,<0, ie. u’<uy. In these circumstances 
the front edge of a positive pulse will be stretched out and it will be the back edge 
that is steepened to form a shock wave. At the temperature where p, =p,, about 
1-96°., uw’ is equal to uw, so far as second order terms are concerned, and any 
distortions can only be calculated by taking into account third and higher order 
terms. 

These predictions are on the whole borne out by experiment. Plate III shows 
the effect, at 2:12°k., of steadily increasing the input power, the thermometer 
being held at a fixed distance of 15-6cm. from the heater. In Plate III (a), (0) 
and (c) the power flows are 0-65, 1:7 and 3-6 watt cm respectively, and it is clearly 
visible that for the larger power flows it is the leading edge that becomes elongated 
and the trailing edge that is steepened. (In Plate III alonea positive temperature 
excursion is shown by a positive deflection.) Plate [V shows the signal for a 
fixed power amplitude (3-5 wattcm™* at 2-12°k.) as the thermometer is pro- 
gressively withdrawn from the heater; at long ranges the behaviour again 
corresponds to theory but at short ranges (<5cm.) there appears to be an 
anomaly in that the leading edge is steepened as though the temperature were 
below 1:96°K 

Significant quantitative measurements cannot be made for these very small 
signals, so that it is impossible to say whether the observed distortions agree 
quantitatively with the theory, but qualitatively and in,order of magnitude the 
agreement is satisfactory. 

In the neighbourhood of 1:96°K.+0-1° (i.e. where p,=p,) a square pulse 
input results in the propagation, more or less without change of shape, of a signal 
with both its front and back edges elongated. ‘This phenomenon is clearly not 
covered by the second order theory. 


§5. SHOCK WAVES 
The preceding section has dealt only with the progressive distortion of a 
wave of finite amplitude up to the stage at which the wave front reaches its 
limiting steepness. Temperley (1951), in addition to giving the Riemann 
treatment of this problem, gives also the Rankine-Hugoniot derivation of the 
velocity of propagation of a fully developed shock wave. ‘This velocity u” is 
orn by ja ES re (12) 
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where wv, and v, are the normal and superfluid velocity amplitudes at the top of 
the shock front. Relation (12) has not been investigated in detail experimentally, 
but one observation at 1-:05°K. gives a qualitative confirmation of it. The 
apparent velocity u’’ for a pulse having a power flow amplitude of 1-3 wattcm™ 
and propagated over a distance of 4 cm. was 20-7 m.sec!, as opposed to the value — 
18-8 m.sec-! obtained for u, by extrapolating to zero amplitude the velocities _ 
obtained with smaller power flows. vw, for this power flow should be about 
180 cm.sec~1, so that the difference of 190cm.sec-! corresponds rather more 
closely to v,, than to $v, (v, is negligible). On the other hand, the uncertainty of 
the thermal data at 1:05°K. does not permit accurate calculations, and the — 
calculation of v, from the input power flow is in any case a somewhat doubtful 
procedure for a fully developed shock front. 

This effect is clearly of practical importance in connection with velocity 
measurements at the lowest temperatures. In the instance quoted the apparent 
velocity differed from u, by some 8%, and at still lower temperatures the effect 
will be even more marked for a given power flow, since (from (11)) v, = W/pST. 
It is essential therefore that measurements below 1°x. should be carried out with 
the smallest possible signals. 


§6. ATTENUATION MEASUREMENTS 

Both Peshkov (1948b) and Pellam (1948, 1949) have undertaken 
measurements of the attenuation of second sound with distance. Peshkov’s 
measurements were confined to one temperature, namely 1-63° K., and he deduces 
that here the most important dissipation effects occur at the walls, due to viscosity 
and thermal conduction, rather than in the bulk of the liquid itself. 

Pellam measured the attenuation of a 150 usec. pulse and found that near the 
lambda-point the amplitude varied exponentially with distance. His values of 
the quantity « in the equation 

Ge) =0(0) exp (—ax) eee (13) 
where 6(«) is the temperature amplitude at a distance « from the heater, are given 
Figure in 5. Most of these data are drawn from experiments in which the 
thermometer was withdrawn progressively from the heater and the reduction in 
signal size noted, but a few of the points, especially very near the lambda-point, 
were taken from experiments showing multiple reflection. For these points the 
reflection coefficient was taken to be unity, and the reduction in amplitude of 
successive echoes was ascribed to attenuation en route. 

In the present investigation these measurements have been repeated, using a 
600 usec. pulse, with the results shown by circles in Figure 5. From the 
experimental errors shown by vertical lines in the Figure, it will be seen that these 
results differ from those of Pellam, and they do not conclusively show the existence 
of any attenuation at all at this particular pulse length. The important Fourier 
components of a 150 sec. pulse are in the region of 3kc/s., while those of a 
600 usec. pulse are in the region of 750 c/s., so that it is just possible that the 
greater attenuation measured by Pellam was due to the higher frequencies 
involved in his experiments. ‘To test this, a small 500 sec. pulse was propagated 
at 2:16°K., and its rise and decay times noted both near the heater and 2cm. 
further away. If the rise and decay times increased by 7, because of the extra 
2cm. path length, the conclusion would be that frequencies from about 1/277, 
upwards were being attenuated appreciably within the 2cm. The actual 
figures for W=0-64 watt cm™ and 6 =0-67 x 10-3 deg. are as follows: 
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Range (cm.) Rise time 7p (sec.) Decay time 7) (sec.) 
0-7 130 200 
if) 260 170 
Arp = +130 Aryp=—30 7 


On account of the effect described in §4, it may be calculated from the values 
of W and @ that the rise time ought to increase with distance, and the decay time 
ought to diminish, by an amount 7, = 85 usec. for 2cm. It is impossible to avoid 
this effect without using imperceptibly small signals. If it is assumed that Ar, is 
composed of +, +75, and Ar, is composed of 7, — 72, we have 7, +7, = + 130 psec., 
7, —T,= — 30 usec., whence 7, =50 usec. and 7,=80 psec. The value of 7, is in 
good agreement with the calculated one and it seems possible therefore that 
7, does indeed represent an integration of the pulse, in which case it follows that 
frequencies upwards of 3 kc/s. (=1/277,) are appreciably attenuated within 2 cm. 
at 2-16" K. 

Although this treatment is only tentative pending further experimental 
evidence, it suggests that the short pulse length used by Pellam should have 
caused some genuine attenuation due to the removal of high frequencies, but it is 
possible that at least part of the attenuation observed by him (Figure 5) was due 
to finite amplitude effects. ‘hese are much more serious in a short pulse than 
in a long one, because the back edge may become a ‘ shock front’ (see §5) and new 
dissipative effects will then come into play which may make the attenuation 
much more severe. From the values of W and @ quoted above, it may be seen 
that very small signals indeed (~10~*deg.) would have to be used in order to 
avoid finite amplitude effects at all safely. 

This same kind of treatment can probably account for one other phenomenon, 
namely the difficulty of observing signals at temperatures between 2-16°K. and 
2-185°K., even though in this region the calculated value of the characteristic 
impedance y, 1.e. (pcu,)~! is increasing apparently without limit. The value of w, 
for a given temperature amplitude @ is becoming very large and the value of 1, is: 
approaching zero, so that shock wave effects will become more and more important 
as the A-point is approached. Signals of any observable size will therefore be 

- subject to very large dissipative processes and will quickly be attenuated. 
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ABSTRACT. Measurements of the magnetic susceptibility of copper sulphate penta- 
hydrate have been made at temperatures between 20° kK. and 1°k. on a single crystal 
specimen. It is found that the susceptibilities measured in different directions with respect 
to the crystallographic axes follow a Curie—Weiss law, with different Curie constants but a 
fixed Weiss constant of —0-6°. The discrepancy between this result and that of Krishnan 
and Mookherji is attributed to temperature-independent paramagnetism. 


§1. INTRODUCTION 

HE first measurements of the susceptibility of copper sulphate pentahydrate 
| as a function of temperature were made by de Haas and Gorter (1930), who 
found that the paramagnetic susceptibility per mole of a powdered specimen’ 
could be expressed as y =0-457/(T+0-7°) between room temperature and liquid 
hydrogen temperatures. From a series of very delicate measurements of the 
anisotropy of the susceptibility of a single crystal of copper sulphate between 
room temperature and liquid air temperature, Krishnan and Mookherji 
(1936, 1938) deduced that copper sulphate is approximately magnetically uniaxial 
with a minimum susceptibility y, along the magnetic axis and a maximum 
susceptibility y, perpendicular to the axis. Combining their own results with 

those of de Haas and Gorter, they deduced that 


¥,=0-399)/(T—2:0°) and x, =0-480(f also) Soeee (1) 


While the variation of the Curie constant with direction can be related to the 
crystal structure of the salt, the concomitant variation of Weiss constant from 
+2-0° to —1-8° is a surprising result, which would be difficult to account for 
theoretically. In this paper we describe a more direct measurement of the 
magnetic susceptibilities and Weiss constants in different directions of a single 
crystal, at temperatures between 20°K. and 1°x. The measurements show no 
significant variation of Weiss constant with direction. The discrepancy between 
this result and that deduced by Krishnan and Mookherji is attributed to 
temperature-independent effects. 


§2. CRYSTALLINE AND MAGNETIC PROPERTIES OF 
COPPER SULPHATE 


Crystallographic examination by Beevers and Lipson (1934) has shown that 
the unit cell of copper sulphate pentahydrate contains two ions, each of which is 
surrounded by six negatively charged oxygen atoms forming a nearly regular 
octahedron. Four of these belong to water molecules and are arranged in a 
square coplanar with the Cutt+ion; the other two oxygens, belonging to the SO, 
group, complete the octahedron and are at a slightly greater distance. ‘The 
electric field acting on the Cut* ion is then of roughly cubic symmetry, but with 
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a component of tetragonal symmetry whose axis is the line joining the last- 
mentioned oxygen atoms. According to the crystallographic measurements. the 
angle between the tetragonal axes of the two Cut+ ions of the unit cell is 
approximately 80°. 

For each class of ion, the departure of the crystalline field from cubic symmetry 
gives rise, through spin-orbit coupling, to different susceptibilities parallel (x,) 
and perpendicular (x,) to the tetragonal axes. Taking the angle between the 
tetragonal axes to be 90° and combining the effects of the two sets of ions, the 
magnetic susceptibility of the crystal has an axis of symmetry normal to the 
plane of the two tetragonal axes. In this direction the crystal has the axial 
susceptibility y,=x,, while in any direction in the plane of the tetragonal axes. 
the susceptibility has the equatorial value y,=4(x+,,). This picture, due to 
Krishnan and Mookherji (1938), is in agreement with the paramagnetic resonance 
experiments of Bagguley and Griffiths (1950), which show that the equatorial 
susceptibility is constant to within 1°%, suggesting that the angle between the 
tetragonal axes is very nearly 90°, rather than 80°. The values of the axial and 
equatorial Curie constants derived from their measurements are respectively 
C,,=0:407 and C,=0-480. We take these figures as the most accurate available 
values. Figures derived from this type of experiment do not include 
temperature-independent magnetic effects. 


§3. EXPERIMENTAL METHOD 

Our experiments were made on a single crystal ground roughly to a sphere, 
and also on a powder specimen. Measurements on the sphere were made along 
three directions, denoted L,,; Lz, L,. ‘The direction L, (154°, 64°, 51° with 
reference to the crystallographic axes) was the direction of axial susceptibility ; 
L, (70°, 41°, 69°) was along one of the tetragonal axes; L, was a bisector of the 
angle between the tetragonal axes. Susceptibility measurements were made by 
an A.C. mutual inductance method at temperatures from 20°K. to 1°K. The 
apparatus and method used have already been described elsewhere (Benzie and. 
Cooke 1950). The measurements from 20° kK. to 8°K. were made by condensing 
hydrogen over the specimen and pumping to various temperatures, and 
measurements from 4°K. to 1°K. were made in a similar way using helium. The 
temperatures were measured by the vapour pressures of the hydrogen and of the 
helium. Temperature and susceptibility measurements were much more 
reliable in the liquid helium region, and most weight has been given to these. 
The measurements at liquid and solid hydrogen temperatures were used as. 
little more than checks on the liquid helium observations. 

In each experiment the sphere was mounted with one of the chosen directions. 
parallel to the axis of a mutual inductance coil surrounding the specimen. 
Changes in the susceptibility y in this direction give rise to proportional changes 
in the mutual inductance M, and if m is the difference of the value of M at a 
particular temperature and M,, the infinite temperature or zero susceptibility 
value of M, we have m=ky, where k isa constant. (Strictly, m is a measure only 
of the temperature-dependent part of the susceptibility, but this distinction can 
be neglected in our measurements. In the liquid helium range the diamagnetic 
susceptibility, for example, is less than 0-1% of the total susceptibility.) If the 
susceptibility y= C/(T—A), the value of the Weiss constant A can be obtained 
by plotting a graph of 1/m against T and extrapolating to 1/m=0. The values 
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of M, and A were in fact found by successive approximation, by extrapolating 
graphs of M against 1/T to 1/7'=0, and of 1/m against T to 1/m=0. The final 
graphs of 1/m (or 1/x) against T are shown in the Figure. It will be seen that the 
graphs for the three directions L,, L, and L, all tend to a value of Weiss constant A 
very near to —0-6°. The diagram also shows the result for a powdered sample 
for which the Weiss constant appears to be slightly greater. The value of A 
obtained with the powder (—0-65°) is very nearly the same as that found by 
de Haas and Gorter (1930) and is probably more reliable, as their lowest 
temperature of measurement was 14°k. A value of —0-65° also fits Reekie’s (1939) 
observations as well as the value of —0-7°. However, the difference between our 
powder and crystal results is within experimental error, and the results can best 
be stated as A=—0-60+0-05 for all directions of the crystal. The law 
y = C/(T'+0-6) holds down to approximately 1-8° k., and between this temperature 
and 14°x. the extreme scatter, from the mean, of y(7+0-6) is less than 1% in 
each of the four series of measurements. 
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In-order to rotate the crystal, the apparatus had to be dismantled and 
reassembled between the measurements along the different axes, and as a result 
the slopes of the graphs are not accurately inversely proportional to the 
corresponding Curie constants. Our accuracy from this point of view is not 
better than +1%. (This, of course, does not affect the accuracy of the 
determination of the Weiss constants.) Further, the value for the powdered 
specimen had to be adjusted for the mass and shape of the sample being different 
from that of the single crystal. ‘Taking these factors into account the relative 
slopes of the graphs are in satisfactory agreement with the Curie constants derived 
from the measurements of Bagguley and Griffiths (1950), which are 0-407 for 
the axial Curie constant, 0:480 for the equatorial, and 0-457 for a powder. These 
are in the ratio 1:1:18:1:12, which compares with our corresponding ratio 
1:1-19:1-12. It must be pointed out that this good agreement is to some extent 
fortuitous, as our value for the Curie constant in the equatorial plane is the mean 
of two values (axes L, and L;) which differ by 3%. 


The Magnetic Susceptibility of Copper Sulphate 127 


Combining our measurements with those of Bagguley and Griffiths we 
obtain 
_ 0-480 0-407 
Kes T+06 ’ .C— 706 Oe ae  idecos (3) 


The isotropic Weiss constant is readily attributable to exchange interaction 
(Wright 1949). 
§4. DISCUSSION 

The considerable discrepancy between the above results and those derived 
by Krishnan and Mookherji is attributable to temperature-independent 
paramagnetism. Our measurements give only the temperature-dependent part 
of the susceptibility, while theirs include temperature-independent effects 
which though small are important at the relatively high temperatures of their 
experiments. In the Table the second row of figures gives the values of 
paramagnetic anisotropy at different temperatures, taken from the paper of 
Krishnan and Mookherji, while the third row gives the values calculated from 
equations (3). ‘The fourth row gives the difference of the two sets of figures. 


Table 


Temperature (° K.) DOO OU OO eee ee Om a 200K LOOkm 160s e140) 2001 
Maman (sceer, VIS) Busy 2D Biss BW Bin eee Pe eK) aston ies! 

=, (B. & C.) ZIOM 2252-242, 202.2286. 314 349° «392 -448 694 
Difference 65 66 66 65 66s geil 75 81 88 87 


It will be seen that at the higher temperatures the difference is remarkably constant 
at 66 x 10-® per mole. This is the region of temperature in which the scatter of 
the results of Krishnan and Mookherji is least; at lower temperatures there is 
an appreciable scatter between their experimental results, and the values in the 
table are smoothed values. The values of anisotropy obtained by adding a constant 
66 x 10-6 to our figures of row 3 lie, in fact, within the scatter of their observations 
at all temperatures. We conclude therefore that copper sulphate possesses 
a temperature-independent susceptibility whose anisotropy has a value of 
60 to 70 x 10-®. 

This is in excellent agreement with the theoretical calculations of Polder (1942) 
who calculated the magnetic susceptibility in terms of the splitting of the energy 
levels of the cupric ions. He found that for a copper ion in an electric field of 
tetragonal symmetry the susceptibility along the tetragonal axis, yj, is 

0-376 4A 4-1 x 10-16 
= —— eee 4 
=e (1- He) te “) 
while the susceptibility y, perpendicular to the tetragonal axis is 
0:376 1 A\2 1:04 x 10-18 
err ( Fr a eae 


Here A is the spin-orbit coupling constant, which for the cupric ion is known from 
spectroscopic data to be —852cm™!, and W, and W, are the energy differences 
of higher orbital levels of the Cu** ion from the ground state, produced by the 
crystalline electric field. ‘The two susceptibilities each contain a term following 
Curie’s law and a small temperature-independent term. (Polder’s theory does 
not take account of exchange effects, and so no Weiss constant is obtained.) 
Both the terms depend on the value of the splitting W, or W,, which to some 
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extent must be taken as adjustable constants in the theory, but once the values 
of W, and W, are chosen so as to: give the right values for the Curie constants 
in the main terms, the small temperature-independent terms are fixed. Polder 
chose W,=15,400cm™!, and W,=26,600cm™. ‘Then 

i= +136x10-§ and — AZ XAOS. vise wee (6) 

For the crystal of copper sulphate, with two ions per unit cell having their 

tetragonal axes at 90°, this gives 

— +20x10-® and ——— ES Ao ae es (7) 
The anisotropy of the temperature-independent paramagnetism is thus 58 x 10-® 
on this theory. The agreement of theory and experiment is certainly within the 
margin of experimental error, and affords new evidence of the existence of the 
Van Vleck temperature-independent paramagnetism. 

For a powdered specimen, Polder’s theory gives a mean temperature- 
independent susceptibility of 60 x 10-* per mole. We have examined the effect 
of this term on the experimental results of de Haas and Gorter (1930) and Reekie 
(1939) but it is too small to permit any conclusions being drawn. It is less 
than the diamagnetic correction which has to be applied, and is of the same order 
as the differences between the two sets of experimental results. The effect of 
interpreting the results of de Haas and Gorter in terms of a temperature- 
independent paramagnetism of 60x10-® per mole and a Weiss constant of 
—0-6°, instead of a Weiss constant of —0-7° without temperature-independent 
paramagnetism, is to reduce the Curie constant by less than 1%. 
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ABSTRACT. It has previously been shown by the author that the interaction between 
free electrons and the lattice vibrational field (if strong enough) leads at the absolute zero 
to the formation of a new state. The London equations are now derived in the sense 
conjectured by F. London. 


: Sale 

HE electron theory of metals in the simplest form is based on the 
| assumption that electrons can move freely except for an occasional 
scattering by the lattice vibrations connected with absorption or emission 
of a vibrational quantum. It has recently been shown by the author (Frohlich 
1950a, referred to as I, 1950 b) that this leads to an interaction (the S-interaction) 
between electrons in a similar way as in field theories emission and absorption 
of quanta is closely correlated to the interaction between sources of the field. 
It has been found that the interaction between two electrons can be very simply 
represented in k-space (Wk is the momentum of a free electron): it is repulsive 
when the difference of the absolute k-values of the two electrons is small and 
attractive otherwise; it becomes negligible when this difference is larger than 
20 where o=2ms/h (s=velocity of sound). If o<k the energy range of the 
S-interaction is of the order «a/k for electrons of energy « =fk?/2m. If strong 
- enough this interaction leads in the ground state to a distribution (denoted by f,) 
in k-space which differs from a Fermi f,-distribution where electrons fill a sphere 
of radius ky. In this f,-distribution a shell whose thickness is of the order o has 
been removed by a distance of the same order from the surface of the Fermi 
sphere. In ordinary space this implies that the probability of finding electrons 
with parallel spin at a distance smaller than 1/o =10-°cm. is slightly larger than 
for completely free electrons. Subject to later confirmation this f,-distribution 
has been identified with the superconducting state. It has been shown that a 
condition concerning the strength of the interaction has to be fulfilled in order 
to make the new distribution stable. This condition is not satisfied by monovalent 
metals but holds for most superconductors; also the correct magnitude of the 
energy difference |S| from the normal state was obtained. ‘This energy 
difference | S| should be inversely proportional to the isotopic mass, as has been 
confirmed by recent experiments (Maxwell 1950, Reynolds et al. 1950). The 
methods of perturbation theory were used for the calculations , Bardeen (1950, 
private communication) has now shown that the results can also be obtained with 
the use of a different method. In the present paper it is shown that the magnetic 

properties follow in a simple way. 


130 : H. Frohlich 


§ 2. 

The electromagnetic behaviour of superconductors shows two main striking 
features: the existence of persistent currents in multiply connected bodies, and 
the Meissner—Ochsenfeld (1933) effect in simply connected bodies. Both are — 
described in a satisfactory manner by the phenomenological theory of F. and H. 
London (1935). According to this theory the electric current density j can be © 
written as the sum of a normal (j,) and a superconducting (j,) component, 


j=jhe OS (1) 
The latter component satisfies the London equation 
d* curlj,= —<B--B=curlLA> ) | eee (2) 


Here B is the magnetic induction and A is the vector potential. If the latter is 
gauged in such a way that its normal component A, vanishes at the surface of the 
specimen, then j, = —cA/d? follows from (2). The length d=10-°cm. is denoted 
as penetration depth because a magnetic field can penetrate into a superconductor 
only to this depth. A number », per unit volume of superconducting electrons 
is introduced by the relation 


d= (me|n,e2)}!2, Lee 


Recently in a most interesting analysis, F. London found that it will be sufficient 
to prove equation (2) for a simply connected isolated superconductor “in order 
to obtain the properties of a ring, and presumably also those of an open 
superconducting wire fed at its ends ”’ 

Asa conveniently shaped body we choose a long cylinder of radius 7) with 
the magnetic field parallel to its axis (z-direction). ‘Then B has a z-component 
only, and hence to satisfy 4, =0, A has a d-component A(r) only (fd =azimuth) 
which depends on the distance r from the axis but not on¢org. The electronic 
states are described by the momentum #R, in -the .z-direction, the angular 
momentum fy and a further quantum number v. ‘The electronic wave functions 
are p=4/ Nee", (vr/r9), where N is a normalization factor and J,,, is a 
Bessel function; p is a positive or negative integer. The density of electrons is 


n= Dela *, ae (4) 
where the sum extends over all occupied states. The ¢-component of the 


" momentum is the expectation value of the operator (h/ir)0/06=p,y. Following 
F. London we shall introduce an average momentum field 


=n iP, oh — * Subp. sees (5) 


This expression vanishes when positive and negative p-values have equal 
probability. 

We shall now assume that it is permitted to treat the magnetic field as a 
perturbation given by the operator 


he Ar) &A%r) 
init 7 Ob? Die i) 


For a homogeneous magnetic field this requires that the radius of curvature p 
of most electrons 1s large compared with r, and that A(r)/r+0 as r+0. Ina field 


where Up= 


SR 
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_ of 600 gauss p~10-*cm. for most electrons, and for weaker fields p is corre- 


‘spondingly larger. The density of the perturbed energy can then be written as 


_ the sum of three terms, 


Oe hoa we eh & Gy Ala vee (7) 
sg ANN), ae (8) U,=- AC) Rail (9) 


and U, is the soe in zero order kinetic energy. In the lowest non- 


_ vanishing approximation, p, is (cf. §4) obtained in terms of a redistribution of 
- occupied levels. This is enforced by the magnetic field so as to make Ua minimum. 
_ It is mainly connected with a transfer of electrons from negative to positive 


p-values. Let Az be the number of transferred electrons per unit volume. We 

then expect p, to be proportional to Az, and U,,oc(Az)?; also approximately 
ON a ete ET ET es (10) 

where ¢ is the maximum zero order energy of an electron. The redistribution 


has been studied in detail by Welker (1938), but these details will not be required 
in the present paper. Since Uj is nearly independent of Az we have 


Oe CNS e Ca NS )e orl: wy blaystes sso (11) 
= PA) Pe 
where a A(r) eee hee). hse re (12) 


ONE eta (oa eee ee (13) and Ol Uy Cy] Coy ec ks (14) 
are obtained by minimizing U with respect to Az. Now for free electrons it is 
found that U) vanishes except for an extremely small term which gives rise to the 


normal diamagnetic change in energy and auch, can be neglected here. Then 
using (14) and (8) ne yields 


: n no 15 
a eM ends Ga (<<) ; a Pee eer e (15) 
Hence with (15), (12) and (8), we find from (12) 
Az i \ Ne 
- = (— Dae Fie Re ako 16 
B= 2A) x= = (Fa) de (16) 


It should be noticed that c, is nearly independent of r. 
Since the current density is given by 


a  () 
ee Pa = Ar) id eo 17 
f= = bo— — Alt), (17) 
j=0 follows if Azy=Az,, as must be expected for free electrons in the present 


approximation. 

It has been stressed by F. London (1948 and earlier Sets) that equations (2) 
and (3) with n=n, are obtained from (17) if it can be shown that p,=0, 
i.e. Az=0. ‘More generally, however, if only Az<Az,, with the use of (16) 


we have 5 R 
en a 
ape See stp) ek a 18 
j=-( $2) 40) (18) 
Thus if Az/Az, is independent of r we may define the number n, of superconducting 
rectrons Dy, Re NG) te cients (19) 


For free electrons Az = Az , and hence n, =0. 
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§ 3. 

We consider now the influence of S-interaction on the value of Az, using 
perturbation theory. According to I, at the absolute zero of temperature the 
distribution function f, of a superconductor is obtained from a Fermi distribution fp, 
by shifting a given small number z of electrons in a narrow energy range by an 
energy of equal amount. The energy difference between the two states was. 
denoted by |.S| per unit volume. In this f,-distribution p, =0 as for free electrons. 
Again the magnetic field tends to make this distribution non-symmetrical in 
momentum space by a displacement of, say, Az electrons such that p, 40, because: 
this leads to a depression of energy according to (9). The f,-distribution on 
the other hand was obtained by displacing z electrons from the f-distribution,, 
and z was determined by minimizing the total energy. A further displacement 
of Az electrons must, therefore, weaken the S-interaction energy * and hence lead 
to a further positive term in expression (11) for the energy density. Its order of 
magnitude should be expected to be |S|(Az/z)?. This is confirmed by 
calculations on the energy of some special non-symmetrical states carried out in I, 
equation (4.2) together | with equation (3.30). 

That this expression gives the correct order of magnitude for any sufficiently 
smooth displacement follows immediately from expressions (3.21), (3.22) also. 
Figure 1(b)) of I. According to these equations the energy required to shift a 
single electron by a small energy from the surface of the distribution is larger 
but of the same order as its increase in zero order (kinetic) energy. The total 
energy required to shift Az electrons must therefore be of the same order as U,, 
equation (10). This agrees with the above expression because | S|~¢z?/n apart 
from numerical factors. 

In order to take account of the change of S-interaction we must therefore 
add the above term to expression (11) for the energy density. Thus, since the 
contributions of the two inner surfaces to U,, and U, cancel, 


U=Uy—c,Az+h(cota(Az)?  — avesas (20) 
where toe | Siiet — —  re (21) 


and c, and ¢,y satisfy relation (15). Minimizing U with respect to Az one obtains 
as equilibrium value Az, of Az, 


AS. =Ci/(Cyt 0) Aso.) 9) (22) 


Hence, using (20), (22), (15) and (8) the equilibrium value U, of the energy 
density is given by 


U.=Up- 38 = (2 - ete) =U ue ald A?(r) = 


2 (€, == Xe) Co Co + Xs D Co + Xo 2mc2 : Co + Xs ; 
. (23) 
Also from (12), (15) and (22) we find for the average momentum field 
—, , eA(r) 
(Po). = a1 a, ee ee (24) 


*This makes use of the fact that S-interaction and magnetic interaction can be treated as. 
independent perturbations in the case of our model. In this order of approximation the expressions 
derived in the previous section still hold. 

+ Only the change in S-interaction should be considered here; the increase in zero order kinetic: 
energy is already contained in Ux. 
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Hence, using (17) the current density becomes 


A eAn C e? An, 
J mc ( Co+ =) Now Pre ee i722) 
‘where Lg NES CotOs)? te cts (26) 


Both a, and c, are independent of 7. We have thus derived the London 
equations (2) and (3) which follow immediately from (25) by applying the operator 
curl. From equations (26), (21), (10) and U,,=c,(Az)?/2 we obtain 


ns S 


a Seen coerce 


As mentioned before |S | and ¢z?/n are of equal order of magnitude. ,/n is 
therefore of the order 4, but detailed calculations of U,, and «, would be required 
to establish the exact value. 
It will also be noticed that with (25), (26) and (3) the energy density (23) 
becomes 
U, = 387? |e 


in agreement with the London theory. 

There is no need to prove the existence of a critical magnetic field for 
transition into the normal state, as this follows from the London equations 
with the thermodynamical method developed by Gorter and Casimir (1934). 

Finally it should be remarked that additional terms to S may perhaps be expected 
to occur in the energy difference between normal and superconducting states when 

- the radius 7 of the cylinder is small of the order 1/o or less, for in this case the 
discrete nature of the electron spectrum may be expected to make itself noticeable. 
It happens that 1/o is of the same magnitude as the penetration depth d which 
may make it difficult to separate the effect of finite penetration from those of 
additional energies. Also from the preceding derivation it will be seen that 
small additional terms in (20) do not essentially alter the form of the London 
equations; they might, however, lead to a weak dependence of n, on either 7 or 
field strength. 


§4. REMARKS ON THE METHOD 
The above calculations where the magnetic field as well as the electron lattice 
vibration interaction are treated as a perturbation holds only if 79 is sufficiently 
small, and if A(r)/r remains finite as r>0. The method of minimizing U 
requires that A(7) is a self consistent field due to outside sources as well as due to 
sources in the metal. For, if H, is the density of the Hamiltonian of the 
outside sources plus that of the unperturbed metal, it is required that the integral 


1 ‘ 
{Ho U+ (curl A 


over space is a minimum. Minimizing with respect to A leads to the Maxwell 
equation curlcurl A=47j. Minimizing with respect to Az(r) leads to dU/dAz =0, 
which was used to find Az. 

Above, both current density j(7) and momentum field p(r) were found pro- 
portional to A(r). As mentioned, these quantities are obtained in this order by 
redistribution of electrons in momentum space; perturbation of the wave functions 
would contribute only in higher order. In fact since p,,, has only diagonal 
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elements, and since f,=0, the first order perturbation of the wave function due 
to the first term in (6) (proportional to A(r)) leads to vanishing contributions. | 
to p, or j. There is no perturbation at all if 7) is so small that the magnetic 
field can be considered as homogeneous. For superconductors the extra energy 
term |.S|(Az/z)? is proportional to C?A*(r), where C is the constant measuring 
the interaction between electrons and the vibrational field (cf. I (2.8)). Here, — 
too, redistribution only yields terms proportional to C?A. It should be 
re-emphasized that only for superconductors is this term | .S|(Az/z)? positive. 
It is negative for normal metals, but smaller than U, (otherwise these metals 
would be superconductors). It might lead to interesting magnetic properties. 
of normal metals at very.low temperatures which will be investigated. 

For other models, e.g. for a box with boundary condition 4=0, Le. 
~=sink,« sink,ysink,z, both] and p vanish for the unperturbed wave functions. 
The magnetic energy is then obtained from perturbation only without redistri- 
bution (cf. Papapetrou, Z. Phys., 1937, 107, 387). In this case p,,, has only 
non-diagonal elements. Hence in contrast to the above there are perturbation 
terms in the wave function proportional to A(r) leading to the required order 
in j and p as well as to a A?(r)C? term in the energy density. 

Finally it should be mentioned that for the present method in the case of a 
cylinder it would not be permissible to gauge A by, say, adding a constant. For 
then A(r)/r would diverge as yO and perturbation theory would be excluded. 
In the case of the box, however, A may be regauged in this way because a 
constant A yields no.energy contribution even when perturbation theory is 
applied; second order perturbation of the linear term then just cancels the first 
order perturbation of the quadratic term in A. : 


APPENDIX—NOTES ON PAPER REFERRED TO AS I 
Expression (3.16) for S3 is correct and does not require a further factor } as . 
mistakenly stated in note 2 added in proof. However, a factor 2 is missing on the 
right-hand side of (2.6) defining o as was pointed out to the author by Mr. M. J. 
Buckingham. Hence, factors $ are required in equations (2.12)—(2.15), (2.23), 
(2.24), (3.15). In (3.17) and (6.9) replace [ P(€) — €/F(4v)!3] by 2[ P(é) — 2€/F(4v)48]. 
In inequalities (3.19), (6.7), (6.8) the right-hand side is 2 instead of 1. The 
influence of these corrections on the results is insignificant. Dr. Schafroth 
kindly informs me of a misprint in equation (3.22); the argument of the last 
log term should be multiplied by (y—o)/(y +c). 
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ABSTRACT. Anexpression is derived for the photoelectric current excited at the surface 
of a metal by radiation of frequency below the threshold frequency of the first order photo- 
electric effect. A physical interpretation of the formal results is given and the possibility 
of detecting this current experimentally is briefly discussed. 


ws INTRODUCTION 


N discussions of the theory of the photoelectric effect at the surface of a metal 
the first order approximation alone is normally considered, corresponding 
to absorption of single quanta by the electrons. For light of usual intensities 

and of frequencies greater than the threshold frequency, the higher order 
corrections (corresponding to absorption of two or more quanta) are certainly 
exceedingly small. When the frequency of the radiation is less than the threshold, 
however, the first order current is zero (at zero absolute temperature) and if the 
frequency is greater than half the threshold frequency the second order terms 
give a current proportional to the square of the radiation intensity. 

It is the purpose of this paper to derive a formal expression for this second 
order current and to estimate its order of magnitude. 

These results may be of interest because of the very high intensities now 
available in pulses from flash tubes, also because the radiation field is here treated 
classically, so that some check on the conclusions so reached seems desirable 
either experimentally or from the quantum theory of radiation. ‘The methods of 
that theory do not seem to have been applied yet to processes of this kind. 

The method employed is an extension of earlier work (Makinson 1949) in 
which a formal expression of some generality was derived for the first order 


photoelectric current. 


SOLES COND TOR DER-CURREN T 
(i) The Model 


It will be supposed that the conduction electrons are free in the interior of 
the metal and that near the surface their potential energy may be represented 
by the potential function V(x), x being a coordinate normal to the surface. 
Bardeen has shown that taking account of the exchange energy of the electrons 
in the metal necessitates the use of different effective potential barriers for each 
value of electron momentum. Since in the present discussion this generalization 
merely requires the use of a more complicated notation without altering the 
essential results, it will be assumed that V(«) is the same for all electrons. 

The wave function uw for a conduction electron in the metal satisfies 


(h?/8a?m)V2u + (th/2m)oujot—V(x)u=0. ssa (1) 
For an electron in the state with wave-vector k(h,, kg, ks), 
Dau = 0,0 (x) exp (kay +ik,2—2mibyt/h),. |... +s (2) 


* At present at Physics Department, University College, London | J 
+ At present at Department of Theoretical Physics, University of Liverpool 
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where EF, =hk?/u —hv,, w=82?m/h and hv, = — V(— «) the height of the surface 
potential barrier. ‘Then 4, is the solution of 


ba (2 (a Vie= 0 ee (3) 
for which by =exp (ik,x) +a, exp(—ik,x), «<0, | 
Re ee (3 a) 
=b,,exp(—px), x>0, 


where |a;,,| =1 and p?=pv,—k,?. Without loss of generality k, is restricted to 
positive values, and the normalizing factor «; is supposed chosen so that in the 
interior of the metal 2%,|u,|?=m, the number of conduction electrons per 
unit volume. 

(ii) The Perturbed Wave Equation 


The steady-state equation (1) is perturbed in the presence of incident light 
which we describe classically as plane waves with vector and scalar potentials A 
and ® of frequency v, so that 


(A|p)V2u — (h/2t)du/ot — V(x)u 
= —(the/2mmc){A.Vu+ $uV.A +4 (ime/hc)A*u}—eDu. ...... (4) 
We may take ® =0, but must retain the term in A’. If 9 is the angle of incidence, 
A =.Aa(x) exp { —27ir(t + y sin 9/c} + conjugate, 


where -¥ is selected so that a(x) is dimensionless and has unit magnitude just 
outside the surface, and where the y axis is in the plane of incidence and the 


plane of the surface. 
Let u=u,+¥V,+W,, where uw; is the unperturbed solution (2) and v, is the 


solution of the first order perturbation equation 
(h/p)V?0,,— (h/271)0v,,/0t — V(x)v;, = — (the/2ammc)(A.Vu;, + 4u,V.A). ...... (5) 


The terms in exp {—27i(E),—hyv)t/h} correspond to stimulated emission and 
do not concern us here, so that, putting 


Vp, =Hy( —4rie/hc).0 (x) exp {ikyy + 1kgz — 2nivy sin 9/c — 2ni(E;, +hv)t/h} ...(6) 


and taking v less than the threshold frequency vo, 0,(x) satisfies the equation 


Og" — {7° + (w/h) V(x) }8q = 2b! + Ip i(Redy + Rg.) + 3bxAxs verre (7) 
after neglecting some small terms in v/c, with the boundary conditions 
0,= Ce”, Pak? +py, x<0, 
oe) NEES uate (7a) 
0,=Be', P=—ke+p(y,—v), «>. 


Since, for v<v, 7 is real, there is no net outward electron current corresponding 
to @,, so the first order photoelectric current under these conditions is zero. 
The second order perturbation equation satisfied by w,, is then 


(hp) ey — (h/2m)000,/8t — V(x)20, 
= —(the/2mmc){A.Vv;,+ $0;,V .A + (ime/hc)A2u;,}. 
Rejecting terms corresponding to stimulated emission and putting 
Wy =%,( —Arte/hc)?A? w(x) exp {tkyy + tkgz — 2ni(E;,, + 2hv)t/h —4rivy sin 9/c} 


eee eee 
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__ we find, on using (6) for v;,, that «,(«) is a solution of 


wy” + {8 — (u/h) V(x) }o = a,04' +10, (RoQy + hga,) + 40,0, —,a7/4, ...... (9) 
_ where two small terms in v/c have been neglected*, and 
P=k?+p(2v—v,), P=k?+2ur. 
We require the solution of (9) which represents an outward stream of 
photoelectrons, satisfying the boundary conditions 
w=he™, x<0, 
= (.e"", x>0, 
since the right side of (9) vanishes at large distances both inside and outside the 


surface. 
Let x(x) be the solution} of 


xi +18? — (w/b) V(x) }x0=9, 
which satisfies the conditions 
mae", Reo tae Slay Meer (11) 
= Pe"? .O0¢-"*,  «>0. 
To determine the magnitude of G, and thus the photoelectric current of 
_ electrons which have absorbed energy 2hv from the radiation, we multiply both 


sides of (9) by x, and integrate from — oo to +00. ‘Then from (10) and (11) 
21sG,Q, = L(k, v), where 


L(k, v) = { {4,0 9/ +10 (Rady + h5Q,)+40,09/ —b,a2/4}x)dx. .. 0... (12) 


Now the transmission coefficient of the surface barrier V(x) to the wave 
eee? 15 1s)" 7)/(s|'O,|*)-so that 
Gale Dis) hL(k y\lAi(4is) 0 eee (13) 
The photoelectric current density from the surface corresponding to w,, is 
given, using (8), (10) and (13), by 
aJ , =(2eh/4amti)(w),0w;,/0x — w;, Ow ;,/0x) 
= —e|a;,|2D(s)(47e2/hmc?l)?.A4(4arml/h}| L(k,v) |? wwe eee (14) 
The expression (14) represents the contribution to the current density arising 
from the electrons from the state ¢;, and the total second order current is then 
this expression summed over all the occupied states of the free electron gas in 
the metal. 
§3. DISCUSSION 


(i) Physical Interpretations 


In the integrand in (12) the first three terms involve the wave function 6, 
describing an intermediate singly excited state, whereas the last term does not. 
We may regard the latter (arising from the term in A’ in the perturbation 
equation (4)) as corresponding to double excitation directly. 


* Tne corresponding terms neglected in (7) do not affect the second order approximation, 
since their contribution to the current vanishes when the first order current vanishes. 

+ This procedure is the same as that used by Makinson (1949), cf. Mott and Massey (1949, 
p. 107). It is not applicable here unless the first order current is zero. 
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In the interaction of an electromagnetic field with a charged particle, the 
corresponding A? term leads in classical mechanics to a second harmonic term in 


the motion of the particle. The condition that one may neglect the term in A? in” 


calculation of the classical motion of the particle (as one neglects it in quantum- — 


particle cum classical-field calculations of the first order photoelectric current) 


is that the second harmonic term in the motion shall be much less than that 


with fundamental frequency. This condition may be written eA/c<me, i.e. that 
the maximum velocity acquired in one period is very much less than c. ‘Thus 
relativistic equations would be necessary in classical mechanics for sufficiently 
large A. However, in the present treatment the same limit would appear in the 


failure of the perturbation method. Since we are only concerned here with — 


second order current under conditions where much greater first order currents 
are absent simply because the frequency lies below their threshold, the 
non-relativistic perturbation theory we have used should be adequate. 

Inspection of (7) and (12) indicates that neglect of the term in a? will not 
greatly alter the value of |L|?. The remaining part of L can be expressed 
approximately in terms of first-order excitation functions (Makinson 1949) 
corresponding to two transitions each involving absorption of a single quantum. 
Writing M(k, v) for the excitation integral 


[Cache +b (hatty + hat) + Abus }xud 
where x, satisfies x," —{r?+(u/h)V(x)}y,=0, 
Xq- ae 2 here N< 0, 
=g,e"" x > 0, 


it can easily be shown, by the method used for determining G, in (10), that in (7 a) 
2iqC, = M(k, v). Denote by ¢, the solution of (3) when 7? replaces p? so that 
gq replaces k, in (3a) ; this is what the wave function would be for a conduction 
electron in the metal of energy equal to that of a singly excited electron. Let 
M(q, v) now: denote the integral corresponding to M(k, v) when 4, replaces $, and 
x, replaces yx. 

Comparing the forms of 6, and C,¢, in «<0 and noting that they do not in 
general vary violently in magnitude near the surface, we may assert that in 
computing L(k, v) we shall arrive at a result of similar magnitude if we replace 
8, in (12) by C,d,. 

Thus | L(k, v)|~|C,|-| M(qg, v)| and therefore 

| E(k, v)]~| MCR, ») || MUq »)[I(29). nee (15) 

Introducing an excitation function E(k, v) =(ae2/hmc?q).?| M(k, v)|?, we 
may write as an approximation for the expression (14) for the current density of 
doubly excited electrons : ; 


dJy~ —e. 4|a,|2 « D(s) . E(R, v)E(q, v).4nmf(hg). ....-- (16) 


This may be compared with the corresponding expression for the first order 
current (Makinson 1949), 


aJ,=—e . 4|a;,|®.. D(z) SERS ae eee (17) 


Equation (16) thus describes an excitation from the unperturbed state ¢;, to 
the state ¢, (with absorption of a quantum hy) followed by a further excitation 


fears 


J mere 


ao 


The Second Order Photoelectric Effect at a Metal Surface 139 


_ to the state w, and transmission through the barrier. The rate at which the second 
| excitation occurs should be proportional to the density of singly excited electrons 
near the surface, which is inversely proportional to their speed (hg/(27m) inside 
the barrier), and such a factor is found in (16). 


eX 


(ii) Orders of Magnitude 
| The order of magnitude of the second order current may be roughly estimated 
_as follows. Let dJ,(v) be the contribution to it arising from electrons with some 
- particular wave vector k in the metal, where v<yy<2v. Let dJ,(2v) be the first 
order current which would arise from the same electrons from light of 
frequency 2v but of the same intensity Rergcm~* sect. From (14), (15) and (17), 


dT (v) za (c/2nv)PR® 
dJ ,(2v) (c8av)R 


heq 
where F stands for the dimensionless quantity | M(k, v)|.| M(q, v)|/| M(z, 2r)| 
which is of order unity. Thus 


dd (v) (aI ,(2v)~87e2R| (ch2g?v?)~10-2°R 


when values appropriate to sodium have been used. ‘Taking the whole current 
densities J,(v) and J,(2v) as being in the same ratio, and the first order yield of 
sodium as 10-3 amperes per watt, the ‘second order yield’ of sodium is given 
roughly as 10-16 amperes per cm? per (watt cm~?)?. 

The possibility of measuring such a yield depends on the availability of (i) 
a sufficient intensity of light, (ii) a light filter sufficiently selective to reduce the 
first order yield from frequencies above v) to a magnitude below that of the 
anticipated second order yield. Flash tubes of the gas discharge type are stated 
to give a peak output of about 10° watts in the visible spectrum and we may 
suppose (in the absence of detailed measurements) that half of this lies on either 
side of vy and that 2% of the low frequency output can be focused on the cathode 
(1 cm? in area) of a photomultiplier tube via the filter. The latter is then required 
to reduce the intensity of light of frequency above vy to less than 10~® times that of 
light below vp». If this very severe filtering requirement (which would be more 
severe for lower intensities) can be met, then second order currents from the 
cathode of approximately 10-* amperes are anticipated. This corresponds to, 
the emission of about 10° electrons per flash of two microseconds duration. 

The second-order current should thus be easily detectable in a photo- 
multiplier tube with suitable threshold frequency. If the cathode is free of small 
patches of much lower work function, the whole difficulty lies in adequate 
filtration of the light source without too much reduction in intensity of the 
wanted band. 
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ABSTRACT. A model for plastic deformation in soft metals, i.e. metals which are not 
age-hardened or severely cold-worked, is given. It is based on the assumption that in real 
crystals there is always a number of regions present which, by an unknown process, act as 
sources of dislocations. Furthermore, it is assumed that in soft metals dislocations are very 
mobile except in a limited number of regions which act as obstacles. The model is applied 
to explain a number of facts connected with plastic deformation. 


§1. THE ASSUMPTIONS 


N the course of an experimental investigation of the plastic deformation of 
copper at stresses below the yield point (Kuhlmann and Masing 1948), the — 
following model of the mechanism of creep, involving two assumptions, was 

adopted and subsequently developed in greater detail. 

The first assumption is that in every crystal of a plastic material there is always a 
number of regions which act as sources of dislocations of one type and one sign only, 
probably dislocations in which the edge component is dominant and the screw 
component plays a small réle only.t 

Each source may be characterized by the limiting stress necessary to put it into. 
action, which will be quite low for the most active sources. The process by which 
the dislocations are released is obscure. Possibly one has to think of some kind 
of notch effect, but it is doubtful whether this alone is a sufficient explanation. 

Frank (1948) recently proposed a reflection mechanism for dislocations, 
involving edge and screw dislocations of both signs, which would explain the crea- 
tion of dislocations, the formation of slip bands and the sudden extension on 
applying a stress to a metal. ‘The model considered here, whilst not in conflict 
with this theory, will be applicable at stresses too low to cause reflection, and also 
probably to metals in which Frank’s process has come to an end, 1.e. after some 
degree of cold work. 

The second assumption is that on each active glide plane there are several 
times as many obstacles hindering the movement of the dislocations as there are 
sources, and that the stress required to move the dislocations through these 
obstacles is higher than that required for their creation. 

If the metal is not age-hardened or severely cold-worked, the potential energy 
of one dislocation on an average glide plane P, may then have the form sketched 
schematically in Figure l(a). Mott and Nabarro, on the other hand, worked out 
the theory of plastic deformation in severely disturbed crystal lattices such as are 
probably present in precipitation hardened alloys. ‘The potential energy Py is then 
as drawn in Figure 1(b). In fact there will be a continuous gradation between 
the two extremes. In the present paper, however, only the case of Figure 1(a) is 
treated. 


* Now Mrs. Wilsdorf, N.P.L., P.O. Box 395, Pretoria, South Africa. 
} For the elastic theory of dislocations see, for example, Cottrell (1949). 
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While it has not yet been possible to give a satisfactory theoretical explanation 
for the first assumption, the second one seems to be well justified, as we know that 


_ in real crystals there are, or can be, inclusions, precipitates, internal stresses, dis- 


locations on other glide planes crossing the one considered, grain boundaries, 
mosaic and polygonization boundaries and layers (for example oxide) on internal or 
external surfaces, all of which are obstacles against the movement of dislocations. 
Most of these obstacles act symmetrically, that is, the magnitude of the stress 


necessary to make the dislocations pass them is independent of the side from which 


they are approaching. ‘This property of symmetrical action is of great importance 


in connection with after-effects and fatigue. 
L+n'r/2— 


LAL Om/E 


Pq (potential energy ) 


SE 


Figure 1. Potential energy of a (a) (b) 


ge ge nction: of Figure 2. (a) Elastically bent slab; (b) slab in which the 


its position: (a) in a soft d 4 ‘ 
eS (6) in ses alee same radius of curvature is established by the presence 
Oe age-hardened thetal of a row of dislocations. 


§2. THREE CONSEQUENCES OF THE ASSUMPTIONS 
AND ELASTIC THEORY 


(1) According to the modern theory of dislocations (cf. Cottrell 1949) two like 
dislocations distant d apart on the same glide plane each produce a shear stress 
7 =A/d atthe centre of the other. Here Ad =GA/27(1—,), where G is the shear 
modulus, A is the atomic spacing in the glide plane measured in the glide direction, 
and yu is Poisson’s ratio. These stresses produce a force 7A on each dislocation, 
tending to drive them apart. If one of them is fixed, an applied shear 7=A/d ts 
thus necessary to keep them at a distance d from one another. ‘Therefore we 
expect that by the action of the shear stress 7, m dislocations can be piled up 
between two obstacles separated from each other by a distance /, where n is given 
approximately by n~/r/A. 

The exact evaluation of is difficult and has not yet been solved completely. 
With the assumptions that the first dislocation is held by an obstacle which does not 
exert any forces on the following dislocations, that there is a repulsive shear stress 
A/d between two dislocations tending to repel them, and that the external shear 
stress + acts on each dislocation, Heilbronn (private communication) found 
n>Dz/2A, where D is the distance between the first and the last dislocation. 
n~ Dz/2A is a good approximation for largen. Nabarro (1950) found the upper 


limit n<1+Dr/A. 
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Thus so far there is uncertainty about the exact value of the number of 
dislocations N which can be piled up on a glide plane of length L; but at least for 
large N there is a linear connection between L and N, and in the following we shall 
use Heilbronn’s approximate formula and set N=L7/2A. 

(ii) The dislocations in each slip band are considered to have originated from 
one source and, therefore, are of one type and one sign only, for instance, positive 
edge dislocations. The bending stress set up by these dislocations will be partly 
relieved if the neighbouring slip bands contain like dislocations, but not if—in our 
example—they contain negative edge dislocations. Therefore there is a high” 
probability that a slip band which develops between two adjacent slip bands with 
like dislocations will contain the same sort of dislocations as its neighbours. ‘Thus 
we expect whole packets of slip bands with one sort of dislocation only, as was 
observed by Nye (1949). 

(iii) Cottrell (1949) has shown that when m dislocations are piled up behind an 
obstacle under an applied shéar stress 7, their pressure “‘ produces a high local 
stress in the vicinity of the obstacle of m times the external stress”. In our model 
this stress will, in equilibrium, be /7?/2A. 


§3. THE MOVEMENT OF THE DISLOCATIONS 


The dislocations, which according to the assumptions of the first section are 
formed ona single atomic plane or on a few very closely spaced planes, will move 
one after the other along their glide planes until the first one from each course is 
held up by an obstacle. The subsequent dislocations will now distribute them- 
selves behind the first ones ina way determined by the resolved shear stress and their 
mutual repulsive interaction. If the shear stress or the number of piled-up 
dislocations is sufficiently large, the first dislocation will pass the obstacle and leave 
the part of the glide plane under consideration, either moving on to the next 
obstacle or, if there are no more obstacles, leaving the glide plane altogether 
through the surface of the metal. After the first dislocation has surmounted the 
obstacle, the following ones will move up in the ‘ queue’, and as a consequence 
another dislocation will leave the source; the next dislocation will then pass the 
obstacle, and so on until either the second free region of the glide plane is filled up 
and the process stops, or the first dislocation passes the second obstacle and moves 
on to the third one. 

Obviously the number of obstacles in comparison with the number of sources 
in action, the mean distance between the obstacles, and their opacity will depend 
on the state of the material and can change during plastic deformation. These 
differences will give rise to the observed differences in the mechanical properties. 


§4. THE AG-EFFECT 


Before discussing the model in greater detail, we must deal with a special 
aspect of it. Masing (1944, 1948) has pointed out that every. model of plastic 
deformation which involves the presence of an appreciable number of dislocations 
in the material will lead to an apparent modulus of rigidity which is lower than that 
for an ideal crystal. In our case, under a given stress there will be a certain number 
of dislocations in each slip band, divided into groups by the obstacles. On 
relieving or reversing the stress there will be a rearrangement of the dislocations 
within each group, but there will be no appreciable interchange between the groups 
if the obstacles are symmetrical in the sense of § 1, unless the real resolved shear 
stress is nearly equal and opposite to that which was required to move them into 
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_ the glide plane. ‘The rearrangement of the dislocations, however, will cause a 


reversible deformation in addition to the elastic strain. 
If the stress is reduced from 7 to zero the dislocations will rearrange and on the 
average move by q/, where q is a pure number which must be less than}. The 


additional shear strain caused by this rearrangement will be Ay =(N/W)(qlA/L), 


where W is the average distance between neighbouring slip bands normal to the 
glide plane and the other symbols have the same meaning as in §2. Let the 


__ observed modulus of rigidity be denoted by G’, the elastic shear strain by y; we 


find ea G 
re? y+Ay — 1+ Ay/y" 
Thus Ay —G _ GNqlr 


y COS WG. 
and with the relation N= L7/2A of § 2, Ay/y =Iqn(1—p)/WS1-21/W. 

Changes of several per cent in the modulus of rigidity or Young’s modulus 
after cold work, or abnormally low values of these quantities by factors as large as or 
greater than two, have been observed repeatedly (see for example Kawai 1930, 
Laurent and Eudier 1949, Crussard 1945, Maddin, Mathewson and Hibbard 1948, 
Raffelsieper 1950), and there is no obvious disagreement between the above formula 
and experiment. On the other hand, detailed experimental investigation of the 
AG-effect could give very valuable information about the magnitude of / and its 
dependence on degree of cold work, grain size and so on. 

As a rough estimate of the maximum number of dislocations which can be 
piled up in each interval we may take the ratio of theoretical yield stress to the true 
breaking stress at low temperatures, provided failure does not occur by cleavage. 
Until more information is available it therefore seems reasonable to think of values 
for n from 1 to 10 or 20, but hardly higher. 


§5. FACTS WHICH MUST BE EXPLAINED BY A THEORY OF 
PLASTIC DEFORMATION 

A theory of plastic deformation should explain, or at least not contradict, the 
main facts connected with slip. ‘They are: (i) the behaviour of slip bands; (11) the 
stress-strain curve and its dependence on grain size and temperature; (ii1) the 
creep curve under constant stress; (iv) the great influence which thin surface 
layers can have on the creep properties of single crystals; (v) the lattice bendings, 
which often develop in the course of plastic deformation and which manifest 
themselves in asterism, deformation bands and kinks; (vi) Bauschinger effect, 
after-effects and recovery. 

It is thought that the two assumptions of § 1, together with their three conse- 
quences, which led to the model of the movement of the dislocations described in 
§ 3, are in qualitative or quantitative agreement with these facts. In the following 
sections we shall try to justify this. 


§6. SLIP BANDS . 
The existence of slip bands follows at once from the first assumption, but their 
fine structure, which was discovered by Heidenreich and Shockley (1947, 1948), 
‘apparently cannot be explained on the basis of this model alone. 
' The fact that slip bands often start and end within the crystal suggests that the 
sources need not necessarily lie at the surface of the specimen or the grain 
boundaries, but may lie within the crystals. 
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It is difficult to think of a way in which internal sources of dislocations of one 
type and sign only can act without producing fine cracks, but such cracks would 
probably be filled up to a certain extent by slip in the neighbouring parts of the 
metal. This difficulty would be removed if we assumed that the sources release 
pairs of dislocations or dislocation loops. The assumption, however, that only 
such sources exist, or that they are in the great majority, at present seems to lead 
to several other difficulties, and we shall therefore not consider it. We shall 
continue to consider sources of dislocations of one type and one sign only, and 
believe that this is no great restriction because the other type of sources can be 
regarded as pairs of sources of opposite sign next to each other. 

The number of active sources, and hence of slip bands, at a given stress in this 
model will, if it depends on the temperature at all, increase slowly with increasing 
temperature, because if there is an energy of activation involved in the starting of a 
slip band at all it is assumed to be small. ‘The energy of activation necessary to 
move the dislocations through the obstacles, on the other hand, is assumed to be 
larger than the former, and therefore the number of dislocations which leave an 
average source at a given stress within a given time will increase more rapidly with 
increasing temperature than the number of active sources. ‘The strain due to the 
dislocations of one slip band increases proportionally to their number and the 
length of the slip band. Thus while with increasing temperature the number 
of slip bands for a given stress will slowly increase, for a given strain it will. 
decrease. ‘This is in agreement with the result of Yamaguchi and Togino (1929) 
for aluminium single crystals. 

It is thought, however, that the spacing of the slip bands and the fact that they 
tend to be arranged regularly (Kuhlmann 1950) is not a consequence of the 
number and distribution of the sources alone, but that some sort of stress mechanism 
and/or an energy principle, as proposed by Orowan (1941) and Bragg (1942, 1948), 
is acting. ‘Thus, especially at stresses well beyond the yield stress, the number of 
sources whose threshold stress is smaller than, or equal to, the applied stress is 
presumably much higher than the number of slip bands. 

The following energy consideration may be useful in this connection : the 
elastic energy of a sheet of metal which is bent elastically to a given radius of 
curvature is higher than that of a similar sheet in which the same curvature is 
produced by a row of parallel edge dislocations lying in the neutral plane when its 
thickness exceeds a certain value. ‘This may be seen as follows: Let the thick- 
ness of the slab be 2d, its length in the glide direction L, and its width one atomic 
distance (Figure 2). ‘The elastic energy of the slab is (Figure 2 (a)) U,, = Ldo,,2b/3E 
if +o,, 1s the stress at its upper and lower surface. The number of dislocations 
which will produce the same radius of curvature (Figure 2 (b)) is n’ =2Lo,,/EA and 
their energy U,=n'P, where P, the energy of one dislocation per atomic plane, is 
a few electron volts (Koehler 1941). For a given stress o,, the energy of the slab 
with the dislocations will be smaller than that of the elastically bent slab if U;< U 


el» 
pau at Lho,2b/3E>2L0,,P/EX or  b>6P/Noyy. 


For o,,=3 kg/mm? and P=3 ey. 6 is then greater than about 10-4cm., while for 
aluminium single crystals at a resolved shear stress of 3 kg/mm? Yamaguchi (1928) 
found that the thickness of the lamellae between the slip bands was about 
2x10-*cm. Moreover, Yamaguchi found a linear increase of the density of slip 
bands with increasing shear stress, while the present calculation entails-,a 


On the Theory of Plastic Deformation 145 


proportionality between 1/b and o,,; and finally the conception that bending 
_ stresses play an important réle in plastic deformation seems to be sound because, 
__as has already been mentioned in § 2, and will be discussed in $9, the dislocations 
themselves set up bending stresses. 
Thus it is very tempting to pursue this idea, which may possibly lead to a theo- 
retical explanation of the sources of dislocations, but it needs some modification 
and, therefore, so far, is to be regarded as a suggestion only. 


§7. THE STRESS-STRAIN CURVE 
(1) Yteld Stress 


The magnitude of the yield stress for unstrained metals in this model is 
governed by the number of obstacles and their opacity. In agreement with 
generally adopted ideas on the structure of real crystals, we assume that either 
or both increase sharply with decreasing grain size and crystal perfection. 

Let us follow the ideas of Orowan (1934a) and take the yield stress, 7,, as that 
stress under which the dislocations start to pass the first obstacles at such a rate that 
the velocity of deformation equals the testing velocity. The formula derived by 


"Becker (1925) may be applied here : 
W=K exp{—V(r_—7)2/2GRT} 


It gives the probability that within the volume V the applied stress 7 is increased to 
Tp by thermal activation. K isa constant, R and T have the usual meanings, and G 
is, as before, the modulus of rigidity. Thus if it is necessary for the stress 7p to act 
within the volume V for a dislocation to pass through an average obstacle, where in 
most cases Tp will be much smaller than the theoretical yield stress, the probability 
that this will happen under the applied stress 7 is given by Becker’s formula. 
Therefore the yield stress will be reached when W assumes a certain value, say 


W,, and we get UE Wee Gers)? 

BT fea GR 
In §2 it has been pointed out that the stress in the vicinity of an obstacle behind 
which dislocations are piled up is v7, and if / is taken not only as a measure for the 


average distance between two successive obstacles but also for that between the 
sources and the first obstacles, we have nr =/7?/2A. This leads to the equation 


2 
dl ( ) , so that 7,2 = 2A (rn— (2GRTp/V)%}, 


P= 7GRT\"® 2A 
or, if we set (2Arp/l)¥?=7, Vrp?/2G=U), andU,/kRp=T,, we have 
Ty=T {1 —(T/Ts)"?}!2. This formula is in qualitative agreement with experi- 
ment, because with decreasing grain size and crystal perfection / will decrease and rx, 
increase and, therefore, both 7) and 7's will increase with decreasing grain size. 
Ts is the temperature at which, in our equation, the yield stress becomes zero, 
but it is not the melting point, although it will be near to it. This is obvious 
from the fact that 7, depends on p and therefore on the testing velocity. With 
_ decreasing testing velocity p decreases slowly and hence 7’, increases. 

As we have taken all constants to be independent of temperature, which 
certainly isincorrect for G, and possibly for V and rx, the above formula will only be 
a first approximation. Another complication arises at high temperatures because, 
as Orowan (1935) has pointed out, Becker’s formula is not valid for very low stresses. 
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There a second term, of the form Kexp{—V(tg+7)?/2GRT}, which takes 
account of the stress fluctuations in the unfavourable direction, and which 
cannot be neglected. On this basis Kochendérfer (1941) explained the 
striking result of Schmid and his co-workers (1926, 1929, 1930) that near the — 
melting point the yield stress of single crystals of cadmium, bismuth, zinc and mag- 
nesium remained constant, or even increased, with increasing temperature. In 
view of the strong influence which surface layers may have on the plastic properties 
of single crystals, however, we think that this effect was probably due to oxide 
layers, especially as at that time the effect which surface layers may have was not — 
known, and, therefore, apparently no precautions against oxidation were taken. 
In Figure 3 the above-mentioned results of Schmid and his co-workers have 


Figure 4. Plan and cross section of 


Critical Shear Stress, ty ( gm/mm’) 


Figure 3. Experimental points for bismuth (x), specimen as used in the creep 
zinc (@) and cadmium (+) are from Schmid experiment by Andrade (1949). 
and co-workers (1926, 1929, 1930). The If the thickness of the annulus is . 
curves represent ry =79{1—(T/Ts)"2)} "2, where d at a distance r from the axis, 
Tg in all cases is taken as the melting point. then d/dy=r 2/7. 


been used to test our equation for 7,(7). The curves are calculated with 7; equal 
to the melting temperature, although a better fit can be obtained if Ty is chosen 
freely. 

For very soft materials and low speeds of deformation the critical shear stress is 
often taken as the stress which produces the first detectable deformation. Here 
another consideration may be useful. Until now we have calculated as though the 
first dislocations would move right up to the obstacles even under very low 
stresses, but in fact some types of obstacles—as for instance an oxide layer—will 
give rise to an image force which in the most extreme case could have the same 
effect as if the dislocations which approach them were repelled by their mirror 
image with respect to the oxide layer. In §2 we mentioned that two like disloca- 
tions on the same glide plane repel each other with a force Ar =AA/d; hence no 
dislocations at all can leave the sources at stresses below 7=A/2/. With the 
experimental value r,~5 gm/mm? for the softest single crystals we find from this 
formulal> GA/4n(1 —4)7y ~3 x 10-3cm. This is a reasonable value, but because 
of our ignorance of the nature of the sources and the image force of the obstacles 
the above consideration may only be regarded as a qualitative suggestion. 
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(ii) Work Hardening 


For single crystals of cadmium at low strain rates Andrade and Roscoe (1937) 
observed a linear law of work hardening and found that the number of slip bands was 
independent of strain. ‘Therefore they suggested that the law of work hardening 
for each individual slip band is linear, and pointed out that this, together with the 
result of Yamaguchi (1928) that in aluminium single crystals the number of slip 
‘bands increases linearly with the shear stress, would also explain the parabolic law 
of work hardening in aluminium. 

The mutual interaction of the dislocations in the slip bands leads to such a 
linear law of work hardening because, as was stated in § 2, the equilibrium number of 
‘dislocations in each individual slip band is N=LZ7/2A. The strain produced by 
them is y=}NA/W=LA)r/4AW, 1.e. increases linearly with stress, as soon as the 
length of the slip bands L has reached its final value. The hardening coefficient 
is thenh =d0r/dy =4AW/AL, where W, as before, is the average distance between two 
‘successive slip bands normal to the glide plane and is assumed to be constant. 
However, this cannot in general be the only cause of work hardening, though it 
may be so in favourable circumstances. The fact that the stress in single crystals 
‘can increase to five times the yield stress or more, and the marked temperature 
_ dependence of the work hardening coefficient at low temperatures, especially in 
hexagonal single crystals, show this. Therefore we must look for other reasons for 
work hardening and, apart from the trivial redistribution of stress during deform- 
ation (stress hardening = Spannungsverfestigung), in our model these can only be 
an increase in number and height of the obstacles which we shall call ‘real’ work 
hardening. Among the kinds of obstacles mentioned in § 1 only slip bands of other 
glide systems, internal stresses and mosaic, polygonization and grain boun- 
daries can increase in number or opacity, or both, during plastic deformation, and 
we cannot believe that this increase is in general very strongly dependent on the 
temperature. But there is another obstacle we have not yet mentioned which 
was first introduced by G. I. Taylor (1934): the mutual interaction of disloca- 
tions in neighbouring slip bands of the same glide system, especiaiiy when they are 
unlike and, therefore, move in different directions. It is thought that this inter- 
action may play an important réle at low temperatures, but the second consequence 
in §2 should be borne in mind. In the present paper we shall not attempt to 
investigate the matter in detail nowever, and shall only discuss it briefly later. 


§8. THE NORMAL CREEP CURVE AND THE INFLUENCE OF 
SURFACE LAYERS ON THE CREEP RATE 

The normal creep curve under constant stress consists of three parts—the 
sudden extension, the transient creep and the steady creep (Andrade 1910, 1914). 

A sudden extension due to the movement of the dislocations from the sources to 
the first obstacles on each active glide plane is to be expected, but this explanation 
can only be valid for sudden extensions below a shear of at most one-tenth of one 
per cent. However, the sudden extension at higher stresses may rise to much 
greater values because some rapid creep in the most efficient slip bands, which 
is completed almost immediately after the application of the load, will be 
indistinguishable from the sudden extension as defined above. 

The transient creep is due to the gradual decrease of the creep rate as the slip 
bands are being filled up with dislocations and as ‘real’ work hardening goes on. 
Where mosaic boundaries are the dominant obstacles the real work hardening will 
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arise as follows: Before the deformation these boundaries will have positions such. 
that the blocks fit between their neighbours as well as possible; the passage of 
dislocations through a boundary will shift the part above and below the glide plane 

relative to the lattice into a necessarily less favourable position, and consequently 

the energy of activation for the passage of dislocations will rise, at least for the first 

few. 

In a single crystal the glide plane extends through the whole specimen and the 
dislocations will finally leave the specimen through the surface. For poly- 
crystalline material the situation is rather more complicated, as the glide pianes end 
at the crystal boundaries. The dislocations, when leaving the grains, create steps. 
in their surfaces, to which the neighbouring grains set up an increasing resistance, 
and either the creep comes to an end for this reason or, if the applied stress is. 
sufficient, equilibrium may be reached by the initiation of a suitable complex slip 
in the neighbouring grains. As before, the dislocations can ultimately leave the 
specimen only through the surface. Thus steady creep means that a dynamic 
equilibrium is reached in which as many dislocations leave the specimen through 
its surface as are formed. Naturally the rate of steady creep will be determined by 
the number and opacity of the dominant obstacles. 

If the external stress is too small to force the dislocations through the highest 
obstacles, or if for geometric or other reasons the dislocations cannot leave the 
specimen, no steady creep will occur. The latter case is illustrated in a recent 
experiment by Andrade and his co-workers (1949)* : 

‘“‘ A polycrystalline metal wire loaded by a device that keeps the longitudinal 
stress constant shows, in a suitable temperature range, a sudden extension, followed 
by a creep that can be divided into a transient component and a permanent or 
quasi-viscous component. Experiments have been carried out on creep under a 
different system of stressing. ‘The metal is in the form ofa circular disc in which is. 
cut a deep coaxial groove, so that a thin annulus of metal separates an inner and an 
outer part of the thick disc. In the stressing apparatus the outer part of the disc 
is firmly held, while a constant torque is applied to the inner part of the disc,. 
which is gripped between metal plates. The cut face of the annulus is subjected. 
toauniform shear. [See Figure 4 of the present paper. ] 

‘“Under these conditions, when the surfaces through which the stress is. 
applied are not free, the law of flow with the two metals so far tested, lead and 
cadmium, are modified. ‘There is a transient component obeying the Andrade. 
t/3 Jaw, but the permanent or quasi-viscous flow does not take place.” 

The explanation of this result on the basis of our model is straightforward. 
Under the conditions described the dislocations will on the average travel 
everywhere parallel to the surfaces and, therefore, cannot leave the specimen, 
and hence no steady creep can occur. 

Andrade and Randall (1948) discovered that the creep rate of acadmium crystal 
suddenly changes if the crystal isimmersed in certain electrolytes. Andrade (1949), 
Phillips (1949) and Hall (private communication) independently proved that the 
effect is due to a surface layer, probably of oxide, which dissolves or thickens in the: 
electrolyte. If the surface layer is removed by electrolytic polishing the creep rate 
can no longer be increased by immersing in an electrolyte. The layers concerned. 
in considerably lowering the creep rate are much too thin for it to be possible to 
explain the effect by their mechanical strength. According to the model of creep 


* T am indebted to Professor E. N. da C. Andrade for this description of his experiment. 
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presented here, the explanation is that in these very soft single crystals the surface 
layer is the dominant obstacle, so that when it is removed the creep rate increases 
because now the less efficient obstacles become dominant. In agreement with this 
explanation the effect is absent in polycrystalline material. There the grain 
boundaries are the obstacles which govern the creep rate. 


§9. ASTERISM, DEFORMATION BANDS AND KINKS 


As was stated in §2, it is to be expected that—contrary to the assumption of 
G. I. Taylor (1934)—packets of slip bands with like dislocations will form in the 
course of plastic deformation. Now any block of material which contains a 
number of parallel slip bands with, say, N positive edge dislocations each, is bent 
‘see Figure 5). If Wis the distance between neighbouring slip bands and L their 


7 (Shear Strain) 
Resolved Shear Stress, t (kg/mm?) 


Total Angle of Asterism, 2 (degrees) 


Figure 6. Values of shear stress (x) and shear strain (e) 


Figure 5. Schematic illustra- corresponding to the same experiments are plotted for 
tion of the bending set up varying values of the angle of asterism. Experimental 
by the presence of parallel data from Yamaguchi (1929) for aluminium single 
rows of dislocations. crystals. 


length, the angle of rotation is p = NA/W, where, as before, Ais the atomic distance in 
the glide plane and in the glide direction and, if equilibrium is reached, N =/7/2A. 

This bend will manifest itself as asterism in the Laue pattern, and its axis of 
rotation is the same as that found experimentally by ‘Taylor (1928), Yamaguchi 
(1929) and Burgers and Louwerse (1931), namely, parallel to the glide plane and 
perpendicular to the glide direction. From the foregoing it follows that asterism 
will occur even when great care is taken to avoid inhomogeneous deformation, 
because it is merely a consequence of the presence of regions with a large surplus 
of dislocations of one type and sign, which is a characteristic feature of our model. 
It seems more difficult to find reasons why in some materials, for instance «- and 
B-brass (Elam 1936, Maddin, Mathewson and Hibbard 1948, 1949, ‘T'reuting and 
Brick 1942), the asterism is practically absent. 

As the density of dislocations within the material depends on the stress, but not 
directly on the strain, we expect a close connection between asterism and stress but 
not between asterism and strain. This seems to agree with experiment. Figure 6 
shows the values of p for aluminium single crystals taken from the paper by 
Yamaguchi (1928), plotted against stress and against strain. While the crosses 
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all lie within a narrow region the dots are scattered. Taking for 1/W the 

experimental result of Yamaguchi (1928) for the same material, namely , 
1/W=C(r—7,), where 7,20-8kg/mm? and C=200(mm.kg/mm*)* 

we find p= La(1—p)GC(7? — 17) = 0-1L(72—77,) 


Ot igi 


for L in millimetres and 7 in kg/mm*. The curve plotted in Figure 6 is — 


p=0-027 (7?—0-87), so that L=4mm. 

L is only an approximate measure for the average width of the bent regions, 
for we have neglected the interaction between neighbouring glide bands as well 
as glide packets, and assumed that the equilibrium number of dislocations is. 
always reached. 

Deformation bands, in particular regions in which the lattice curvature 
changes sign, will occur at the boundaries where a region containing dislocations 
of one sign meets a region containing dislocations of the opposite sign. We 
have to expect a greater number of deformation bands at low temperatures 
than at high temperatures, for in §6 we showed that at a given stress the 
number of active sources does not vary very much with temperature, but 
that the average length of the slip band increases markedly, and according 
to §2 each slip band sets up in its neighbourhood a region in which only 
slip bands with the same type of dislocation are likely to be formed. Therefore 
at high temperatures the relatively few sources which act first, and are thus 
almost independent of each other, will be the origin of long slip bands, and 
hence will influence comparatively large regions, while at low temperatures their 
region of influence will be much smaller because there will be a great number of 
sources acting before the slip bands reach an appreciable length and their zones 
come into contact with each other. From then on the number of deformation 
bands will to a first approximation be independent of stress and strain until 
secondary slip starts. As the energy of a number of like dislocations on different 
glide planes, according to elasticity theory, is lowest when they are arranged 
in a wall perpendicular to the glide plane, there will be a tendency for the 
deformation bands also to form planes perpendicular to the glide direction and 
to the glide plane. 

The deformation bands will act as strong obstacles against the movement of 
dislocations because the unlike dislocations lock each other. They thus will be 
a source of real work hardening of the type mentioned in § 7. 

If single crystals of zinc or cadmium with the hexagonal axis perpendicular 
to the specimen axis are axially compressed, a type of deformation called 
‘ kinking ’ occurs (Orowan 1942). It consists of sharp lattice bendings or ‘ kinks’. 
Hess and Barrett (1949) in an experimental investigation found that kink bends 
are a special type of deformation band. They explained kinks as produced by a 
great number of pairs of edge dislocations which under the applied stress 
originate in a section of the specimen and separate, thus producing a great surplus 
of like dislocations within restricted regions. 

The close relation between the ideas of Hess and Barrett and those put forward 
here is obvious. The only difference between kinks and the lattice bendings 
which produce deformation bands and asterism is that the latter arise 
unintentionally even in ‘homogeneous’ deformation, while those which cause 
kinking are of a larger scale and are formed under the influence of macroscopic 
bending stresses. 
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§10. RECOVERY, BAUSCHINGER EFFECT AND AFTER-EFFECT 
. (i) Recovery 

According to §8, steady creep goes on without further work hardening, i.e. 
the process has become stationary. If now a single crystal which has crept at a 
_ steady rate for a while is suddenly unloaded, no more dislocations are created, 
but under the influence of their mutual repulsive forces a movement of the 
dislocations will go on, though with a smaller speed. If the obstacles are of the 
symmetrical type as many dislocations will move back to their sources and be 
annihilated there as will be annihilated at the surfaces of the specimen, i.e. leave 
the specimen through the surface. While this elimination of dislocations goes on, 
the number and opacity of the obstacles will remain unchanged and the net 


__ plastic flow will be zero. But the properties of the material which depend on the 


number of dislocations present (for instance hardness, yield stress, conductivity, 
modulus of rigidity, thermoelectric power and so on) will change: the material 
will recover. 

In an earlier paper (Kuhlmann, Masing and Raffelsieper 1949) such a recovery 
has been considered, and with some simplifications the following law of recovery 
has been derived: x=B—A*ln(t+¢%)). Here x is the number of dislocations 
present in the material and ¢ is the time of recovery, A* is proportional to the 
absolute temperature at which the recovery has taken place, B= A*In(A/K,)+ Kg, 
where K, and K, are constants, and t) = A*/vy where vy is the velocity of recovery, 
dx/dt, att=0. Although the agreement between this law and the experiments 
seems to be good (Kuhlmann, Masing and Raffelsieper 1949, Cottrell and Aytekin 
1950) the author is aware that the calculation is too much simplified and needs 
improvement, and hopes to do this in a later paper. 


(ii) After-effect and Bauschinger Effect 

If, in the course of a deformation, macroscopic parts of the material (i.e. grains 
or regions which are comparable with the dimensions of the specimen) behave 
elastically, while others undergo plastic deformation, large-scale internal stresses 
will be set up (Wartenberg 1918) such that after unloading some part of the 
deformation will slowly disappear. This reversed creep is called after-effect, and 
its time and temperature law has been derived in a recent paper (Kuhlmann 1947). 
It is the same as for recovery. 

Similarly, but without the help of large-scale internal stresses, an after-effect 
will occur when a specimen which has been strained but has not yet reached the 
stationary state of steady creep is unloaded. In this case the dislocations will 
be preferentially eliminated at their sources, at least in the initial stages of the 
process, and thus a small plastic deformation will be connected with the recovery. 
As before, its direction will be opposite to that of the original deformation. We 
therefore expect an after-effect (or reversed creep) in all cases where steady creep 
has not been reached. The after-effect independent of large-scale internal stresses 
can be calculated in very much the same way as that which occurs when they 
are present. | 

The Bauschinger effect is closely connected with the after-effect because the 
tendency of the material to creep back spontaneously and partially to annihilate 
the preceding deformation will naturally lead to a lowered critical yield stress of 
the material when the direction of the stress is reversed. Masing (1923, 1926, 
Masing and Mauksch 1925) has treated the case where large-scale internal 
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stresses are involved. ‘The corresponding case without large-scale stresses, 
investigated for instance by Sachs and Shoji (1927) for «-brass single crystals, 
can be understood in the light of the above considerations on after-effect. 

Apparently no after-effect has yet been observed in single crystals, although, 
according to this model of plastic deformation, it is to be expected in certain cases. 
The reason may be that it takes place too slowly. 


§11. EXPLANATION OF A SPECIAL EFFECT CONNECTED WITH SUDDEN 
ALTERATIONS OF THE STRESS IN CREEP EXPERIMENTS 


The best way to test the present model of plastic deformation seems to be to 
employ it for the explanation of creep experiments, preferably in the region of 
steady creep. Here, as discussed in § 8, we do not expect changes of the number 
and opacity of the obstacles or the number of the active sources. In the steady 
creep range we, therefore, on the one hand avoid the great complications 
connected with asterism, deformation bands, true work hardening and changes 

in the number of slip bands, and on the other hand cannot employ these in a 
theoretical explanation of the observations, so that an agreement in this region will 
be of greater significance. 

Fortunately there are a few experimental investigations which are suitable for 
testing the model, namely the creep experiments employing small changes of 
load during steady creep carried out by Carreker, Leschen and Lubahn (1949), 
by Phillips (1949) and by Cottrell and Aytekin (1950). The last investigation is 
the most complete of these, and we shall discuss it. 

When the stress on single crystals of zinc was lowered by different small 
amounts Ar after the steady creep rate had been reached, it was found that for a 
certain period of time the creep rate was reduced to a very small value and then 
increased to a steady value again. The time 5 which elapsed before the creep 
rate reached its new steady value was in approximate agreement with that expected 
from the rate of recovery 7 calculated from the relation r=vh (Orowan 1934b), 
where v is the rate of steady creep and h is the hardening coefficient, i.e. 8 =Az/vh. 

Both the creep rate v and the rate of recovery 7 (determined by strain hardening 
and annealing experiments) were proportional to exp{—(U)—«r)/RT}, in which 
%isaconstant, 7 is, as usual, the resolved shear stress, and U, an energy of activation. 
The magnitudes of U, and « were found to be the same in creep and in recovery, 
and the coefficients of proportionality were essentially the same when allowance 
was made for the presence of the hardening coefficient. 

The idea adopted by Cottrell and Aytekin was that while creep had been going 
on at the higher stress, the yield stress remained constant because of a balance 
between strain hardening and recovery, and that when the applied stress was 
suddenly lowered the creep rate would be small until the yield stress had fallen 
by the same amount, due to recovery. 

Obviously the experiments of Cottrell and Aytekin are in agreement with our 
model, in so far as, according to § 10, the same energy of activation U, is involved 
in creep and in recovery, namely that necessary for a dislocation to pass an 
average obstacle, and, moreover, in the cases where recovery takes place under 
stress, as in the experiments of Cottrell and Aytekin, the energy supplied by the 
acting stress, «7, is the same for recovery as for creep. 

To understand these experiments quantitatively let us make a simplified 
model of the process and consider a slip band of length LZ with a sequence of 
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v obstacles, to a first approximation of equal opacity and regularly arranged at 
intervals, 7=L/v. The equilibrium number of dislocations in such a slip band is 
(s2) N= Lnfl=L7/2At. 

The dislocations will leave the slip band at a rate proportional to 
exp{—(U)—«7)/RT} and will be replaced by new dislocations from the sources 
in the manner discussed in § 3. 

Let the stress now be reduced by Ar =7/n to 7’ so that (n—1)L/I dislocations. 
are the equilibrtum number with the new stress. There will now be a surplus 
of v dislocations in the slip band. The first one will leave the glide plane through 
the surface of the specimen more rapidly than would be normal for the new 
stress, but this just brings the density of the dislocations to its equilibrium value 
in the last segment, and no dislocation from the preceding one can move up. 
As the strain produced by moving dislocations is proportional to the total distance 
traversed by them, and as the first dislocation only causes a slight readjustment 
of the dislocations between the last two obstacles, it will lead only to a negligible 
strain. The next dislocation will leave the surface after a life-time corresponding 
to the steady creep rate at the new stress and will cause all dislocations in the last 
section to move up and enable one dislocation to pass over from the penultimate 
to the last segment, but for the same reason as above it will not be replaced from 
the preceding section. The slip produced by the second dislocation to leave 
the surface after the stress was lowered will thus to a first approximation be A/v, 
where, as before, A is the interatomic distance in the glide plane in the glide 
direction. This process will go on, the rate at which the surplus dislocations 
leave the surface being Kexp{—(U)—«7’)/RT} while the slip produced is 
(k—1)A/v for the kth dislocation. When the vth dislocation has left the glide 
plane, the new equilibrium number of dislocations in the slip band is reached 
and thereafter each dislocation which leaves the slip band will, as before, be 
replaced by one from the source and produce the slip A. Thus this simplified 
model leads to a creep velocity which rises linearly with time until the new steady 
rate is reached, and the strain increment A during that time will be half as large 
as it would be if the final velocity v had acted all the time. 

Though it cannot be expected that the creep rate will necessarily increase 
in a truly linear manner from zero to its final velocity, as this simple model 
predicts, the second prediction, that the strain increment during the transition 
time will be A~4v6, should be fulfilled. ‘The reason is that we start with a certain 
number of surplus dislocations, of which practically all leave the specimen at the 
rate appropriate to the new stress, but on the average they travel only over half the 
length of the slip bands. The reverse effect to that discussed, namely that a sudden 
increase of the stress gives rise to a transient rapid creep, can probably be explained 
on the same basis. There, in addition to the normal creep, the dislocation density 
is brought up to its new equilibrium value. ‘The present calculation is equally 
applicable to polycrystalline material and to single crystals, but to single crystals 
with an oxide layer only as long as other obstacles besides the layer are present. 

Expressing the result in symbols, we find that the timed which elapses before the 
last surplus dislocation has left the average slip bandis6 =(v/K) exp {((U)—«7’)/RT}, 


+ As n=I7/2A dislocations are piled up behind each obstacle there will be a stress magnification. 
factor n, but the difference between U=U,—«7, the experimental result from Cottrell and Aytekin,, 
and U=U*—7?, which follows from the present model, would hardly have been noticeable because: 
the stress range in which measurements of U for single crystals are practicable is rather narrow, so 
that, taking into account the normal scatter of experimental values, the curve U(r)=U*—Br™® 
cannot be distinguished from its tangent U=U,)—«r. ‘The same applies to any other reasonable: 
function U(r)= U*—f(z). 
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or in general, since N=L7/2A, 5=(LA7/2AK)exp {(U)—«7’)/RT}. The new 
steady rate is v=(A/W)Kexp {(U,—a7’)/RT}. For the present case there is no 
change in the number and opacity of the obstacles during plastic deformation. 
We also assume that W is constant during deformation; according to Andrade 
and Roscoe (1937) this is true for single crystals of cadmium at very slow strain 
rates. ‘The same is true, according to Hoyt (1927), for single crystals of zinc 
provided a certain minimum deformation has been exceeded. For such a case 
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Figure 7. Experimental points by Cottrell and Aytekin (1950) for the creep of a single crystal of zinc. 
After the crystal had crept for some time with a constant rate at 82° c. and tr=27-7 gm/mm? 
the stress was reduced to 26:3 gm/mm? at t=0. (The curve is that predicted in § 11.) 


we found in §7 the relation h=4AW/AL. From this we get 8=2Ar7/hv. This 
formula is in good agreement with the result of Cottrell and Aytekin, 8=Az/hv. 
Considering the time law, we find that the shear during the transition time is 
given by y=(v/28)t?. In Figure 7 the readings for one of Cottrell and Aytekin’s 
creep experiments are plotted together with that curve. It will be seen that the 
agreement is very satisfactory. 


ACKNOWLEDGMENTS 


The author wishes to express her thanks to Professor N. F. Mott for extending 
to her the facilities of the H. H. Wills Physical Laboratory and for his kind 
interest and encouragement. ‘Thanks are also due to Messrs. F. R. N. Nabarro 
and F. C. Frank, and especially to Mr. G. Wyllie of that laboratory, for helpful 
discussions, and to Professor E. N. da C. Andrade and Professor A. H. Cottrell 
for kindly allowing the use of some unpublished results. 


REFERENCES 


AnpraDE, E. N. pa C., 1910, Proc. Roy. Soc. A, 84, 1; 1914, Ibid, 90, 329; 1949, Lecture on 
Surface Effects in Plastic Deformation; Conference on Plastic Deformation and 
Recrystallization, Cambridge, September 1949; see also The Strength of Solids 
(London : Physical Society), 1948, p. 20. 

ANDRADE, E. N. pa C., and RANDALL, R. F. Y., 1948, Nature, Lend., 162, 890. 

ANDRADE, E. N. pa C., and Roscog, R., 1937, Proc. Phys. Soc., 49, 152. 

Becker, R., 1925, Phys. Z. Sowjet, 26, 919. 


On the Theory of Plastic Deformation 155 


Brace, L., 1942, Nature, Lond., 149, 511; 1948, Proc. Camb. Phil. Soc., 45, 125; 1948, 
The Strength of Solids (London : Physical Society), p. 26. 

Burcers, W. G., and Louwersg, P. C., 1931, Z. Phys., 67, 605. 

‘CARREKER, R. P., LESCHEN, J. G., and LuBAHN, J. D., 1949, Trans. Amer. Inst. Min. Metall. 
Engrs., 180, 139. 

CoTTRELL, A. H., 1949, Progress in Metal Physics, 1, 77 (London : Butterworth). 

CorTTreLL, A. H., and AYTEKIN, V., 1950, #. Inst. Met., in the press. 

CRUSSARD, C., 1945, Rev. Metall, 42, 286, 321. 

Exam, C. F., 1936, Proc. Roy. Soc. A, 153, 273. 


_ Frank, F. C., 1948, The Strength of Solids (London : Physical Society), p. 46. 


HEIDENREICH, R. D., and SHOCKLEY, W., 1947, ¥. Appl. Phys., 18, 1029; 1948, The Strength 
of Solids (London : Physical Society), p. 57. 

Hess, J. B., and Barrett, C.S., 1949, Metals Trans. (A.I.M.M.E.), 185, 599. 

Hoyt, S. L., 1927, Trans. Amer. Inst. Min. Metall. Engrs., 74, 116. 

Kawal, T., 1930, Sci. Rep. Tohoku Univ., 19, 209. 

KKOCHENDORFER, A., 1941, Plastische Eigenschaften von Kristallen und metallischen Werk- 
stoffen (Berlin: Springer), p. 81. 

KOEHLER, J. S., 1941, Phys. Rev., 60, 397. 

KUHLMANN, D., 1947, Z. Phys., 124, 468 ; 1950, Z. Metallkunde, 41, 129. 

KUHLMANN, D., and MasInec, G., 1948, Z. Metallkunde, 39, 361. 

KUHLMANN, D., MasInc, G., and RAFFELSIEPER, J., 1949, Z. Metallkunde, 40, 241. 

LAURENT, P., and Eupirr, M., 1949, C. R. Acad. Sci., Paris, 228, 225. 

Mappin, R., Matuewson, C. H., and Hisparp, W. R., 1948, Trans. Amer. Inst. Min. 

Metall. Engrs., 175, 86; 1949, Metals Trans. (A.I.M.M.E.), 185, 527. 

Masinc, G., 1923, Wiss. Veréff. Siemens Konzern, 3, 231; 1926, Ibid., 5, 135; 1944, 

Nachr. Wiss. Gottingen, p. 81; 1948, Z. Phys., 124, 586. 

MasInc, G., and Mauxscu, W., 1925, Wiss. Veréff. Siemens Konzern, 4,74. 

Mott, N. F., 1946, 7. Inst. Met., 72, 367. 

Mort, N. F., and Nasarro, F. R. N., 1940, Proc. Phys. Soc., 52, 86; 1948, The Strength of 

Solids, (London : Physical Society), p. 1. 

Nasarro, F. R. N., 1946, Proc. Phys. Soc., 58, 669; 1949, Z. Metallkunde, 40, 81; 1950, to 
be published. 

Nye, J. F., 1949, Rev. Métall., 46, 371. 

Orowan, E., 1934a, Z. Phys., 89, 605; 1934b, Ibid., 89, 614; 1935, Ibid., 98, 382; 1941, 
Nature, Lond., 147, 452; 1942, Ibid., 149, 643. 

Puiuires, D. J., 1949, Thesis, Bristol. 

RAFFELSIEPER, J., 1950, Z. Metallkunde, in the press. 

Sacus, H., and SHoji, H., 1927, Z. Phys., 45, 776. 

Scump, E., and Boas, W., 1929, Z. Phys., 57, 575; 1930, Ibid., 61, 767; also GrorciErFF, M., 
and ScuMID, E., 1926, Ibid., 36, 759; FAHRENHORS?T, W., and ScuMID, E., 1930, Jbid., 
64, 845. 

‘Taytor, G. I., 1928, Trans. Faraday Soc., 24, 121; 1934, Proc. Roy. Soc. A, 145, 362, 388. 

TREUTING, R. G.,and Brick, R. M., 1942, Trans. Amer. Inst. Min. Metall. Engrs., 147, 128. 

VON WARTENBERG, H., 1918, Verh. dtsch. phys. Ges., 20, 113. 

YaMacucutl, K., 1928, Sci. Pap. Inst. Phys. Chem. Res. Tokio, 8, 289; 1929, Ibid., 11, 151. 

Yamacucul, K., and Tocino, S., 1929, Sci. Pap. Inst. Phys. Chem. Res. Tokio, 9, 277. 


tI-2 


Factors Controlling Dislocation Widths 


By A. J. FOREMAN, M. A. JASWON anp J. K. WOOD 


Department of Mathematics, Imperial College, London 
Communicated by H. Jones ; MS. received 21st Fuly 1950 


ABSTRACT. The work of Peierls and Nabarro is extended to a family of edge 
dislocations of greater widths. ‘The weaker the shearing forces between adjacent atoms 
in the slip plane, the wider is the dislocation. The external shear stress required to move 
the dislocations is extremely sensitive to the width, becoming vanishingly small at widths 
of the order of three atomic spacings. The theory is applied to bubble raft dislocations, 
and satisfactory agreement is found with experiment. 


§1. INTRODUCTION 
HE early theories of the edge dislocation (Taylor 1934) were based on the © 
| classical elastic dislocations of Volterra, in which singularities of stress and 
strain occur at the centre. A more direct approach was made by Peierls 
(1940) and Nabarro (1947), who took into account the atomic structure along the 
slip plane whilst still making use of classical elasticity theory elsewhere, and in this 
way determined the atomic displacements throughout the dislocation. Their 
calculations, however, were based on an assumed form for the interaction forces 
between adjacent atoms in the slip plane. It is the purpose of the present paper 
to extend this work by showing how the width of a dislocation and the external 
shear stress required to move it depend on the nature of the assumed interatomic 
forces. It is concluded that dislocations in single pure annealed crystals have a 
greater width than that obtained by Peierls and Nabarro, and that the external 
shear stress required to move them is less than the observed yield stress. Finally 
an interesting illustration of the theory is provided by the bubble raft dislocations 
constructed by Bragg and his collaborators. 


§2. THE PEIERLS-NABARRO DISLOCATION 

A diagram of the cross section, one atom thick, of a single positive edge 
dislocation is shown in Figure 1, the problem being two-dimensional. Following 
Peierls, it is convenient to regard the dislocation as constructed by slicing the 
material along the slip direction x’ Ox, thus forming two semi-infinite planes. 
The upper plane is cut along Oy and an extra row of atoms inserted. The two 
planes are then ‘stitched’ together along x’ Ox, the upper one being compressed, 
the lower extended, so that there is perfect alignment of corresponding atoms at 
infinity, as illustrated in Figure 1. Each half-plane is now separately in equili- 
brium under the action of the forces applied to the boundaries A’ A, B’ B respec- 
tively, it being assumed, for purposes of calculation, that the planes can be treated 
as isotropic elastic continua. 

The forces acting along A’A which maintain the upper half-plane in its de- 
formed state, are taken to be a system of shear stresses directed towards the origin, 
these being equal and opposite to the shear stresses in B’B. The displacement 
U(x) along the slip direction, caused by these stresses, of an atom at the point . 
x in A’A will be equal and opposite to the displacement U(x) of the corresponding 
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atom in B’B, provided that U(«) varies slowly from atom to atom. It follows 


therefore that the relative displacement r(x) of the two atoms is given by 
r(x) = U(x) — U(x) +4d=2U(x)+4d  (x>0) 
r(«) = U(x) — U(x) —4d=2U(x)-3d — (x<0), 


the initial relative displacement in the undeformed state being +4d. Since 
there is perfect alignment at infinity, 7(+ 0) =0, so that 


oe co) dite oS es (2) 


It has been shown by Nabarro and by Leibfried and Liicke (1949), on the basis 
of classical elasticity theory, that the shear stress p(x) and the displacement U(x) 


at x are related by the equation 
a(1—v) pet ad U(x) ax ) 
ple) =| i ord (3) 


where G is the shear modulus and » is Poisson’s ratio. 


Figure 1. Cross section of a positive edge dislocation. The arrows indicate directions of shear 
stress and displacement for a pair of corresponding atoms. 


The shear stress distribution along A’A is provided by the interaction forces 
between the atoms in A’A, B’B, so that p(x) is a function of r(x). This function 
p(r) has not yet been calculated, but it must be periodic with period d, and 
alternately positive and negative in each half-period. Furthermore, it has the 
limiting form p=Gy7/d for small relative displacements, in accordance with 
Hooke’s law. A simple function of this type, viz. 

ves ce ae sth eee eg re (4) 
was put forward by Peierls (1940) as the most convenient assumption as to 
the law of force. 

It is clear that each assumed function p(r) implies.a corresponding displacement 
function U(«) and a corresponding distribution of shear stress p(x). ‘Thus, 
substituting for 7 in (4) from (1) gives 

| G . 4nU 


ee ie eee 0 5 
p=— 5-sin—G| (5) 
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which, when put into (3), leads to the integral equation for U: 
(=p). 470 [" LUG de 
= 9} sin we = ah ae Pe eae Ho Shoe (6) 
As may be verified by substitution, a solution of (6) which satisfies the boundary 
conditions (2), and which represents a single dislocation, is given by 


U(t) = — £ tant, p=), fy 3 (7) 


From (5) and (7) the distribution afi shear stress which maintains these displace- 
ments is Cas 


is lace prgz ve ee he = (8) 


These displacements and stresses are both odd functions as required, and their 
general shape is entirely satisfactory, but it is unfortunate that they vary rather 
rapidly from atom to atom at the centre. As a result of this, the equation (1) 
for relative displacements is slightly in error at the first pair of atoms. Further- 
more, the strain along the slip direction, 1.e. dU/dx, reaches a maximum of the order 
of 1/5, which is above the limit usually accepted in classical elasticity. It is, 
however, considered unlikely that the errors in the theory are serious on this 
account, particularly as the elastic problem is essentially one of simple tension 
and compression. 

The dislocation described above is such that the atoms in the slip plane become 
fairly well aligned at a short distance from the centre, so that the major part of 
the slip caused by the dislocation takes place over a region of only a few atoms. 
This is classed as a narrow dislocation, the strains near the centre being large. 
It is found convenient to define the width of a dislocation as the range of x for 
which U(x) is numerically less than d/8, i.e. half its limiting value. The 
dislocation represented by equation (7) is then of width 

w=d/(1—v), 
i.e. about 1:5 atomic spacings if »~0-3, which is very narrow. It is probable 
that in actual crystals the dislocations are wider than this. The extreme narrow- 
ness here is due to the assumed law of force (4) between atoms, which is very 
unlikely to hold in real crystals. In consequence, a flatter type of periodic 
function of stress against relative displacement has been suggested. This would 
almost certainly lead to a wider dislocation, for the tendency of adjacent atoms to 
align themselves correctly on either side of the dislocation will not be so marked 
if they are attracted by weaker forces. It has been found possible to investigate 
this question quantitatively, and hence to determine how the width of a dislocation 
varies with the assumed law of force. In order to obviate the difficulty of having 
to solve integral equations of the type mentioned, the method of approach adopted 
here is to choose some suitable type of function U(x), find the corresponding stress 


distribution from (3) and, by eliminating x between them, find p as a function 
of U and hence of 7. 


§3. MODIFICATION OF THE PEIERLS-NABARRO DISLOCATION 

It is required to find displacement functions U(«) such that they (i) represent 
wider dislocations than that obtained by Peierls and Nabarro, and (ii) give curves 
of stress against relative displacement having the same initial slope as the 


previously assumed sinusoidal curve (4), thus ensuring that Hooke’s law holds 
for small values of 7/d. 


i) is given by 
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One type of displacement curve which implies a greater dislocation width 


d C 
— ath — 
U 7, tan ee (9) 
where ¢ =2(1 —v)x/d and a is a parameter greater than unity. When. a=1, this 
becomes the solution (7). The width w of the dislocation (9) is easily: seen to be 
adj/(1—v), and so increases linearly with a. On substituting from (9) into the 
relation (3), the corresponding stress distribution is found to be 


Ge 


Pp om e: e+ jae Fawr 5s OP) ei erererers (10) 
from which the function P(r) is given by 
Gas 
p= Jpg Sin 2arld. SOO (11) 


This is a sinusoidal curve as before, but has an incorrect initial slope, except 
whena=1. This example shows at once that the main difficulty of the problem 
is to find some type of function U(x) which will ensure that the initial slope of 
the corresponding function f(r) is given by Hooke’s law, i.e. p = Gr/d for r/d small. 

A displacement function which satisfies the conditions mentioned, and which 
is of a suitable form, is 


d 
U=- x frat! +(a-Narah, Rievereters (12) 


where €=2(1—v)x/d, and the parameter a determines the width. When a=1, 
U reduces to the function (7), whilst for larger values of a the dislocation is wider. 
Displacement curves are shown in Figure 2 for the cases a=1, 2,4, 10. In these 
wider dislocations U(x) varies more slowly from atom to atom, so that the maximum 
elastic strains are smaller, and the theory accordingly more consistent. 

Setting the displacements (12) into (3) gives the corresponding stress 
distribution, the integration being facilitated by observing that (12) may be written 


as 
d 0 aac 
U= = x4 1-(a-1) 57h tan Be 
and that the operator can be brought outside the integral sign. The resulting 
stress distribution is then 


G 
p= {1=(a-1 she zp 


= - ae g + 2a(a— 


reg NS: 
a+? ) 


1) tml. 
(a? + G*)2 
As the dislocation widens, that is as a increases, it is found that the maximum 
stress decreases, its position moves outwards, and the sharp rise in stress which 

occurs when a=1 gives way to a smoother type of curve (Figure 2). 
It will be noted from (13) that as| £] > ©, the relation between p and £ becomes, 
for all values of the parameter a, 


Be We bat 14 
Leaeey ( ) 
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which is the shear stress along the slip plane in Koehler’s classical dislocation 
(1941). It will also be noted that the limiting form for U as | £|—> co is, from (12), 


a7 red 
U=- mi Ar Ont (¢>0), 
ec (15) 
joe (<0) 
4° Wnt [ 
Hence from (1), for | ¢| large and, therefore, for the relative displacement r small, 
aE Ae eid=ljnt 5 ee (16) 


It follows from (16) and (14) th s 7-0, the relation between p and r is p =Gr/d, 
which is Hooke’s law, as required. 


0 0:25 05 | 
2(1-v) x/d r/a 
Figure 2. The upper diagram illustrates the family Figure 3. Curves of stress against 
of displacement curves assumed. The lower relative displacement. 
one illustrates the corresponding shear stress 


distributions. 


From (12) and (13), using (1), the function p(7) corresponding to the assumed 
displacements is written most conveniently in the parametric form 


ovens {2(2a—1)sin 0 —(a—1)sin 20}, 

Gia oa ‘ 
2nr a—1 2+. °2 = oy eee ( ) 
il =0—- sin 6, 


where @ is given by €=acot6/2. The cases a=1, 2, 4, 10 are illustrated in 
Figure 3. As can be seen, the initial slopes of the p(r) curves are constant for all 
values of a, but the amplitudes decrease as a increases, i.e. as the width increases. 
It is interesting to note that the relation between the width w and the amplitude 
Pmax 18 given very closely by the equation 
Gd 1 
De Pla) oe ae (18) 
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as is illustrated by Table 1. This relation has been found to hold approximately 
for several families of dislocations of the same type as that rcs here> ‘This 


Table 1 
a al D 3 4 5 Om CO 
(1 —v)w/d 1-000 1:255 1-653 2:076 2°507 4-697 0:442a 
(277/G)pmax 1-000 0-787 0-594 0-471 0-389 0-207 2:202/a 


27(1—v)wpmax/Gd 1-000 0-988 0:982 0-978 OL975 0-973 0:973 


suggests that the form of a dislocation depends primarily on the amplitude of 
the corresponding p(r) function rather than on its precise shape. 


§4. EXTERNAL STRESS REQUIRED TO MOVE A DISLOCATION 


We now consider the minimum external shear stress T which will move the 
dislocation through the lattice. In the case of the dislocation discussed by 
Nabarro this was found to be 


ia Ao exp {—2n/(1—»)} a (19) 


which is very much smaller than the theoretical shear strength G/2z of the lattice. 
‘The value of T for the dislocation represented by the displacements (12) may be 
found following the methods of Nabarro (1947); the procedure is outlined below. 

The interaction energy of the atoms in A’A and B’B when the centre of the 
dislocation is displaced by ad has been given by Nabarro as 


Vom Got, > | H+eon {ean 2) (24 5), 


apart from a constant term; the series may be transformed into a more rapidly 
converging one by the relation (‘Titchmarsh 1937) 


i Kn) = 2 : hie VGrcosOata 7) hee (20) 


lhe eas energy for the dislocation (12) may be shown to be 


voiced «hein SE [tool taan—n 22) 


the subsequent analysis being simplified by the use of this operational form and 
by the use of the relation (20). Equating the constraining force F on the 
dislocation to —d-! or gives 


Fs {1- (a-1) 2 += aes 1) tise % sexp {— 2msa/(1—v)} sin 42s, 


in which the series is dominated by the first term. From this it follows, using 
the Mott-Nabarro relation F =rd, that the shear stress is given by 


Pate) re (21) 


where K =2n(a—1)/(1—v). 


162 A. F. Foreman, M. A. faswon and fF. K. Wood 


The shear stress T is extremely sensitive to the parameter a and hence to the 
width w of the dislocation, as can be seen from Table 2. It has been pointed 
out by Cottrell (1949) that the value of T for a=1 is from ten to fifty times larger 
than the observed yield stress for annealed pure single crystals. 


Table 2 
a 1 2 4 10 
p38 w 15d 1:9d 3-1d 71d 


Ip 2x 10-*G 5 <10-7G 2DXAOHG 4x 10-*G 


p=1/4 w 13d 1:7d 2°8d 6:3d 
if: 610 4G 3x10G bx10TVG 1K10T3GC. 


If the yield stress is determined by the force required to move a single 
dislocation which is straight and parallel to a simple crystallographic axis, then the 
above results suggest that the widths of dislocations in crystals are greater than 
1-5d. A dislocation of width 2d would provide a lower limiting value for w, as 
it would require a stress of the order of magnitude of the observed yield stress. 
to move it. 


§5. BUBBLE RAFT DISLOCATIONS 

A useful check on the theory is provided by the edge dislocations which have 
been obtained experimentally in bubble rafts by Bragg and Nye (1947), even 
though the lattice of bubbles differs from the simple cubical structure, and it is. 
a case of plane stress rather than of plane strain.* By shearing of adjacent rows 
past each other Bragg and Lomer (1949) have obtained what is effectively the 
function p(r) for bubbles of 1:3mm. diameter (Figure 3). ‘This function is of 
the same form as those described here except that the region in which Hooke’s. 
law holds is more clearly defined. ‘The amplitude of Bragg and Lomer’s curve 
is rather less than half that of the sine curve with the same initial slope, and is 
very nearly equal to that for the case a=4 described here. As the width of a 
dislocation depends primarily on the amplitude of the function p(r), see equation 
(18), the displacements corresponding to a=4 will provide a reasonable approxi- 
mation to those in the bubble raft dislocation. If a small correction is made for 
the slightly different spacing of the bubbles, and Poisson’s ratio taken as 0-5 for 
1-3 mm. bubbles as indicated by Bragg and Lomer (1949), it is found that the 
calculated displacements of the bubbles in the slip plane of the dislocation are in 
reasonable agreement with the measured values, as indicated in Table 3. 


Table 3 
x/d 1 2 3 4 5 6 7 8 9 10 
2u(x)/d meas. 0°18 0°32 0:39 0:43 0:45 0-47 0-48 0-48 0-48 0-48 
calc. 0:20 0:32 0:39 0:43 0-46 0-47 0:47 0:47 0:48 0-48 


It is interesting to note that the widths of these bubble raft dislocations are 
very dependent on the size of the bubbles used in the experiments, and hence it 
is possible to obtain wide or narrow dislocations by varying the size of the bubbles. 

* The theory for plane stress may be derived from that of plane strain by replacing 1—v by- 


1/(1+¥). ~ . 
+ Dr. W. M. Lomer has kindly supplied a photograph of a dislocation in 1:32 mm. bubbles. 
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| _ For example, the dislocations in 0-7 mm. bubbles correspond to about a=10 in 


the theory, whilst for 1-9 mm. bubbles they are very narrow, being of the type 
a=2. It is hoped to check the theory more fully when the experimental (p, 7) 
curves for bubbles of varying sizes become available. 
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Further Studies — Experimental and Theoretical — of a 
Wide-Angle @-Spectrometer which uses a Prolate 
Spheroidal Magnetic Field 


By T. H. BRAID ann H. O. W. RICHARDSON 
Department of Natural Philosophy, The University, Edinburgh 


Communicated by N. Feather ; MS. received 4th August 1950 


ABSTRACT. A f-spectrometer is described which focuses a conical sheaf of rays 
emitted at about 80° to the axis of a prolate spheroidal magnetic field. In agreement 
with a previous theoretical prediction, the trajectories cross over at two sites in the rotating 
meridional plane. The field behaves in some ways like a magnetic mirror and can be 
arranged to form an aberration-free image of a point on the axis of symmetry. It cannot, 
however, give a clear image of an area because the Abbe sine condition is. not fulfilled. 

The wide-angle focusing action is compared with that recently found by Slatis and 
Siegbahn, for which a different optical analogue is suggested, based on the isolation of a 
meridional focal line between two refractions. 

The adjustment of the spectrometer and measurements of its resolving power and 
solid angle of collection are described, an example being 1:4°% of 47 with a line-width 
of 1/113 at half-height. By working at less than unit magnification it is possible to focus 
rays emitted at less than 80° to the axis. 


§1.-INTRODUCTION 
HE mode of focusing in the present instrument (Richardson 1949) bears 
| some resemblance to the forming of an optical real image by a curved 
mirror. The resemblance appears when the paths of particles emitted 
from a point source on the axis of symmetry are plotted in the usual rotating 
meridional 7, z plane. The particles are repelled by the bowl-shaped meridional 
potential barrier and adjacent ‘ rays’ cross over near the turning point rather in 
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the manner of a specular reflection. ‘The turning back, however, is gradual so 
that it cannot be sharply distinguished from refraction. In spite of this it is 
convenient to describe a system in which a strong repulsion coincides with a 
cross-over as a magnetic mirror. The concept has proved helpful in predicting 
the behaviour of the rays. An alternative description of a similar system in terms 
of refraction is discussed in §5*. 

A remarkable feature of the focusing is that it seems possible to form an | 
image of a particular object point which is corrected for spherical aberration, in 
the sense that the point at which the rays return to the axis is stationary over a 
narrow range of angles of emission. This angular range is defined by stops 
which select a ‘ sheaf’ of rays emitted at a fairly large inclination to the axis. 

The magnetic field has lines of force which are, approximately, ellipses lying 
on confocal prolate spheroids, so that the equipotentials are the orthogonal 
family of confocal hyperboloids of revolution. ‘The field is maintained between 
two soft-iron poles of roughly hyperboloidal shape. 


- Figure 1. Wide-angle focusing in the spheroidal field. A sheaf of mono-energetic trajectories 
represented in the rotating meridional 7, z plane, showing the selecting ring-slits B1 B2 
and B3 B4 

An earlier computation of electron trajectories in this field had predicted 
promising focusing properties with a large solid angle of collection. ‘The theory 
indicated (Richardson 1949) that the focused conical sheaf of rays has the 
property of crossing over at two sites in the meridional plane, as shown in 
Figure 1. The first cross-over, which simulates reflection, takes place at the 
outer envelope or caustic formed by mono-energetic rays of different inclinations. 
The second cross-over can be arranged to be on or near the axis, at a site suited 
for entry into a small Geiger-Miiller tube. 

The site of this second cross-over is controlled by choosing the correct radius 
for the ring-slit B1B2 between the central stop B1 and the ring baflle B2. The 
radius given by the theory was 9cm. for the case of an object point O and an 
image I symmetrically placed 12cm. apart on the axis, in a field in which the 
distance PQ between the foci of the confocal spheroids was 20 cm. 3 

These dimensions were used in the design, the gap between the magnet poles 
being set at 16cm. The equivalent hyperboloids 16cm. apart at the vertices 
are shown in Figure 1 and the actual magnet profile in Figure 3. 


* The measurements made with the ‘ short ’ slit, with Figures 3, 4 and 5, were presented in a 
paper read at the Edinburgh Meeting of the Physical Society of 19th February, 1949. 
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The momentum-resolving action of the baffles can be seen by considering 
first a ‘ central’ ray, defined as a ray falling normally on the mid-plane, half-way 
between the pole-pieces. The central ray of the selected mono-energetic sheaf 
is inclined at 81° to the axis at O. After traversing the ring-slit B1B2 it reaches I 
by a path symmetrical about the mid-plane. The central rays of sheaves of other 
energy would be stopped either by B1 or by B2 which also stop all rays of lower 
energy. The non-central asymmetrical rays of higher energy can still enter B1B2 
obliquely but most of them will be stopped by the second set of baffles B3 and B4, 
which define the angular width of the accepted sheaf. 


§2. COMPARISON WITH A THIN MAGNETIC LENS 


Some insight into the nature of the correction of spherical aberration at large 
angles may be gained by considering the passage of the thin wide-angle sheaf of 
Figure 1 through a thin magnetic lens. The real image would then have a large 
aberration with two focal lines. The first, a focal ring round the axis, is the 
figure of revolution formed from the meridional intersection of the sheaf where 
it touches the caustic, while the second is a short straight line along the axis. 
The radius of the focal ring increases rapidly with the angle of emission « and the 
sign of the aberration is positive because rays with larger « return to the axis 
nearer to O. : 

Let the field now be transformed to the spheroidal form, with a strength 
chosen to focus the mono-energetic sheaf from O on the axis, as in Figure 1. 
The focal ring will have migrated to the mid-plane and, if the sheaf is thin, the 
axial focal line will have shortened to a point. The shortened second focal line 
is now both an axial and a second meridional intersection of the sheaf and so 
constitutes a three-dimensional point image of the particular object point O. 

The point image I will in general only be partially corrected. In the special 
case in which O and I are equidistant from the mid-plane, perfect correction 
is possible if certain conditions of symmetry and reversibility of rays between 
O and J are satisfied. Consider two rays emitted at angles «, and «,, above and 
below the central ray of the same mono-energetic sheaf. Let these two rays be 
incident on the mid-plane at angles 7, and 7, at the same point. Then, if 
1, = —15, the symmetry of the field ensures that the part. of one ray lying on one side 
of the mid-plane will be a mirror image in that plane of the part of the other ray 
lying on the other side. Hence, both rays pass through O and I. Several such 
symmetrical pairs of rays must form an envelope which is also symmetrical about 
the mid-plane and the extent of its symmetrical part will be a measure of the range 
A« over which correction holds good. A point image at I will certainly exist for 
a small A« if the envelope has a horizontal tangent where it crosses the mid-plane. 

A second symmetry property is shown by the central rays which, being incident 
normally, are necessarily mirrored in the mid-plane. If the baffles are removed, 
all central rays from O converge to a perfect image at I, which is achromatic 
since these rays all have different energies. 

In the spectrometer the radius of the ring-slit B1B2 is chosen to be a little 
larger than is appropriate for a point image. ‘The resulting astigmatism 
introduces a second focal ring which is adjusted to fall on the ring-slit 5 around 
the counter (Figure 4). The second meridional cross-over at this ring reverses 
the order of the rays at the axis so that the spherical aberration is now negative. 
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The use of a partially corrected wide-angle sheaf makes possible an increase 
in the solid angle of collection. Suppose that Ax, is the semi-angular width of a 
paraxial beam which can be focused without intolerable aberration and that Ax 
is the corresponding thickness of a wide-angle sheaf. The ratio of the two solid 
angles is Ax)? : 2Axsina. If Ax) and Ax are small and nearly equal it is clear 
that an appreciable gain is possible in passing from paraxial to wide-angle 
conditions of operation. 


§3. THE QUALITY OF THE IMAGE AND THE POSSIBILITY OF 
MAGNIFICATION 

The question arises of the applicability of this kind of magnetic reflection to 
forming sharp magnified images with an appreciable solid angle of collection. 
It is usually considered impossible to correct an electron beam for spherical 
aberration in a system consisting of any combination of electrostatic and magnetic 
lenses. This is because each component of the system always has a positive 
spherical aberration. However, Cosslett (1946) has remarked that an 
electrostatic electron mirror has negative spherical aberration so that a 
combination of mirror and lens is potentially capable of correction. This 
suggests that the correction of a particular zone in the present case may be due 
to some combination of lens and mirror action. 

In order to make a rough estimate of the quality of the image formed of a 
small disc source at O, with its normal along the axis, we adopt a simple 
aspherical mirror as the equivalent optical instrument and leave any lens action 
out of account. This would be fully justified if all the bending occurred suddenly, 
at the first cross-over, as in specular reflection. In fact, the bending has a broad 
maximum near the cross-over and then falls smoothly to zero at the axis,where 
the bottom of the meridional potential valley is flat. However, as the mirror 
gives values for the inclinations of conjugate rays at O and I in rough agreement 
with the computed trajectories, it should also give an approximate description 
of the conjugate rays through neighbouring points, which build up the image. 

The symmetry of the sheaf in Figure 1 with respect to the mid-plane means 
that the part of the aspherical mirror selected by the ring-slit B1B2 must have 
all its centres of curvature in that plane. The fact that this ring zone of the 
mirror is able to form an aberration-free image of O shows that these centres of 
curvature must coincide with those of the particular aplanatic surface for reflection 
from O to I which is tangential to the mirror round the ring. This surface is the 
well-known prolate spheroid with foci at O and I which passes through the slit 
B1B2 and satisfies Fermat’s principle. 

The mid-plane is not the only site of the ring-slit which gives aberration-free 
images of O. The slit may be moved along the axis and its radius adjusted to 
form again a sharp image I. The ring B1B2 now defines a zone of a spheroidal 
aplanatic surface which is no longer equatorial so that the magnification is not 
unity. By considering the image formation using different zones we can get 
what is an ‘ upper limit’ to its quality, because additional defects such as the 
equivalent of chromatic aberration are present in the electron-optical case. 

Proceeding in this way it is evident that the imaging of an area near O is very 
poor because the image of a point distant Ar from O has meridional and sagittal 
intersections which are separated by distances of order 2Ar. The image can 
easily be drawn for rays which intersect the axis, using successive zones of the 
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‘spheroid. It is erect and of unit magnification for the equatorial zone. As the 
reflecting zone is moved towards the image end of the spheroid, the image of a 
disc becomes a cone which finally turns into an inverted diminished disc. Most 
-of the intensity at the image, however, is due to the skew rays which do not 
intersect at the same points as the axial rays. The only case giving a common 
intersection is the final polar site of the zone, when the arrangement becomes an 
ordinary concave mirror giving a paraxial image. 

The image diffuseness is connected with a failure to satisfy the Abbe sine 
‘condition for aplanatic imaging with absence of coma, except in the paraxial 
polar zone near the long axis. Since this zone is obstructed by the magnet 
pole-piece, there seems little practical possibility of getting good images. This is 
in accord with the work of Schwarzschild (1905) and others who have shown that 
aplanatic wide-angle imaging by a single mirror is impossible but that it can be 
‘obtained if two reflections are used, arranged for mutual correction as in the 
reflecting microscope (Burch 1947). 

From the point of view of 6-spectroscopy it seems that the best image of a disc 
will be a diffuse sphere of similar radius 7 so that the resolution will not be 
improved by making the collecting slit S much narrower than r. By working 
with the ring-slit B1B2 displaced from the mid-plane towards the image, a small 
advantage may, however, arise from the magnification being less than unity. 


§4. THE USE OF OTHER SHAPES OF FIELD 


Wide-angle focusing is clearly possible for field shapes differing from the 
spheroidal form. The criterion is that the surfaces of equal vector potential must 
delineate the contours of a bowl-shaped meridional potential barrier with a 
curvature capable of reflecting a diverging beam of incident particles back to a 
single point on the axis. The appropriate sheaf of rays must be selected by 
ring-bafHles with radii chosen to bring the second cross-over on to the axis. 

A quite different shape of field has recently been shown by the independent 
experiments of Slatis and Siegbahn (1949) to give wide-angle focusing. ‘They 
keep the radius of the ring-slit B1B2 fixed and vary the shape of the field by 
adjusting the ratios of the currents flowing in five adjacent short coaxial coils. 
The contours of the meridional barrier may be plotted from the published diagram 
of the magnetic lines of force. 

Two features of the focusing distinguish this from the spheroidal case. First, 
the bending of the trajectories in the rotating meridional plane, instead of being 
greatest near the first cross-over, has two well separated maxima about 8 cm. from 
each of the poles (which are 61 cm. apart). Secondly, the contours in the regions 
of strong bending are convex towards the axis and so have curvatures opposite to 
those of the rays which cross them. This antithesis is characteristic of the short 
magnetic lens but not of the spheroidal field in which the contours are concave 
to the axis between the poles. 

The lines of equal vector potential are very similar to those of the field plotted 
in Figure 2 (a) due to two coils of one turn separated by a distance of twice their 
diameter, obtained by the graphical method of Stormer (1933). It is thus clear 
that wide-angle focusing can also be obtained using two thin magnetic lenses 
placed near the object and image, with an interposed ring-slit B1B2 to select 
the first cross-over in the mid-plane. 
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§5. AN EQUIVALENT REFRACTING OPTICAL INSTRUMENT 


For fields such as that of Slatis and Siegbahn, in which the first ray cross-over 
does not take place in a region of strong bending, the analogy with reflection at an 


aspherical mirror is not close. It is, however, possible with a very simple refracting 


device to reproduce crudely the feature of a cross-over at some distance from 


the axis, taking place between two refractions which bring a wide-angle sheaf 


from O back to an aberration-free point image I. This alternative to the mirror 
is the thick symmetrical lens of Figure 2 (5). 


3 


La) 


Figure 2 (a). The field of two equal circular currents of radius A placed coaxially 44 apart, 
showing the lines of equal magnetic vector potential and equal potential energy V in the 
rotating meridional plane. 


(b) 

Figure 2 (6). An optical analogue of wide-angle focusing in a field of the type shown in Figure 2 (a) 
in which a focal ring is isolated inside a thick lens. ‘The vertical scale is exaggerated, an 
actual set of values in equation (1) being n=2:24, p=12-6 cm., 5;=18-76 cm., 0=23:38°,. 
a= 39-43°. 


Define a ‘ central’ ray of the sheaf from O as a ray which is refracted parallel 
to the axis inside the lens. Then the large positive spherical aberration will 
cause a pencil from O which includes a central ray to intersect again in a focal 
line perpendicular to the meridional plane. The pencil, if rotated around the axis, 
sweeps out a wide-angle sheaf which now intersects in a focal ring. If we choose 
the thickness of the lens so that this ring lies in its mid-plane then, by symmetry, 
all the meridional pencils leaving points on the focal ring are brought to a second. 
intersection on the axis at I, which is both meridional and axial and so is an 
aberration-free point image of O. The focal ring and the width of the sheaf must 
be defined by ring-baffles B1B2, B3B4, as in Figure 1. As the refractive index 
determines the position of the ring, the instrument will select a narrow range of 
wavelength. 7 
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The half-width of the lens can be found easily from the two standard 
equations for astigmatism at a spherical refracting surface (Flint 1936). On 
inserting the condition that the axial intersection is at infinity, appropriate because 
the central ray is parallel to the axis, the equation for s,, the distance of the 
meridional intersection from the lens surface, reduces to 


P DS ne p cot?6 3 (1) 
a een nicos! 0A 72-sin20))2|7-- ie 


Here p is the radius of the surface, m the refractive index and s and @ are as 
indicated in Figure 2 (0). 
This method of selecting an intermediate focal ring and so neutralizing the 
_ aberration of a wide-angle sheaf, should be feasible with more general optical 
systems having many different magnetic analogues.* 

All systems in which the order of the rays is reversed by an intervening 
meridional cross-over resemble a mirror in being contra-current and in failing 
to satisfy the Abbe sine condition with a wide-angle sheaf, so that the images of 
off-axis points will be blurred. 

It is possible that systems in which the first cross-over occurs between two 
regions of sharper bending are less adaptable for other than unit magnification 
than are more mirror-like systems in which the sharpest bending is near the 
cross-over. 

§6. CHOICE OF POLE-PIECE SHAPE 

Unless large diameter coils and cores are used the outward leakage of flux 
round the edges of the pole-pieces distorts the inner part of the field from the 
spheroidal form. ‘This effect can be much reduced if peripheral ridges are used 
as shown in Figure 4. The pole profile was chosen after a study of the 
magnetostatic equipotentials by the electrolytic tank method of solving Laplace’s 
equation. 

A full-scale wooden model was made of a portion of the magnet, with 
conducting surfaces to represent the known equipotentials. It was assumed that 
the pole-faces were equipotentials and that the 9-in. long cylindrical iron cores 
inside the magnetizing coils had a uniform axial potential gradient. The rest 
of the iron circuit was given the potential of the mid-plane half-way between the 
poles. These assumptions were tested roughly by plotting the actual lines of 
force with a compass. 

In order to have good axial symmetry, the return part of the magnetic circuit 
consists of 12 flat soft-iron plates, each 12 in. wide, securely flanged on to 
dodecagonal end-plates of 48-in. diameter. 

The symmetry of the dodecagon makes it possible to reproduce the 
equipotentials by using a model of a 30° sector of the magnet in a tank tilted at 
this angle. 

The shape of the model pole-face was changed until it gave the equipotentials 
drawn in Figure 3 which show only small deviations from hyperbolae in the 
region traversed by the f-rays. 

To avoid changes in the shape of the field with changing intensity of 
magnetization it is necessary to operate far from saturation so that the 


* Bothe (1950 a, b) has recently treated the formation of a focal ring between two thin magnetic 
lenses by assuming that the deviation d of a ray occurs suddenly as it crosses the mid-plane of each 
lens and that d may be expressed as a power series in the radial distance (r—ro) of the point of 
crossing from a corrected ring-zone of radius 7). This leads to expressions for the effect of slit- 
width on the aberrations and on the resolving power. 
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permeability remains high at the pole-faces. With the available coils with 
7-in. diameter cores the upper limit of the spectrometer was thus fixed at about 
1-5 Mev. 


§7.. THE GEIGER-MULLER TUBE 


The Geiger-Miiller tube is designed to accept rays incident perpendicularly 
all round the axis over a range of about 290° of azimuth. The anode wire of 
0-1mm. diameter lies along the axis and the cathode cylinder of 1cm. inside 
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Figure 3. Magnetostatic equipotentials ( ) plotted using the electrolytic tank, compared 
with the hyperbolae (— — -— - - ) of the spheroidal field. 


SE 


FZ 


S 


y 


SESS NNO Pw 


BNA 


Cag? ZI LLL IL LL ELL LILA 


Figure 4. An axial section of the spectrometer showing the preliminary ‘ short’ slit attached to 
its driving rod R giving a small radial motion. 
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and 1-5cm. outside diameter receives the rays through a ring-window 2:-4mm. 
wide cut round its circumference and crossed by three bridge-pieces necessary 
for mechanical strength. The bridges obstruct about 70° of azimuth and one of 
them serves as a seating for a joint in a nylon film of 0-5 mg/cm? which is clamped 
round the cathode cylinder with a rubber gasket seal about 0:3 mm. thick. When 
the tube is filled with argon and alcohol vapour at 12cm. pressure, the plateau 
is about 200 volts long and the background counting rate is 6 per minute. The 
insulating bushes are of quartz and polystyrene and the tube is joined to a gas 
reservoir of about 150 ml. 

Although the cathode is irregular in shape, the plateau for focused f-rays is 
remarkably flat (about 0-01°% per volt), showing that the ions which they produce 
in their passage across the tube reach the middle part of the wire where the 
electric field is highest. In contrast the unfocused f-rays ejected by a y-ray 
source give a plateau with a slope as large as 0-25°% per volt. 

Encircling the ring-window is the narrow collecting ring-slit S of diameter 
2-2cm., on to which the f-rays are focused. 


§8. RESOLVING POWER TESTS 


It was first necessary to find the correct radius for the ring-slit B1B2 to give 
a sharp focal ring round the counter slit S. The value predicted by the theory 
was 9-1cm. but the actual field is not truly spheroidal. To allow the radius to 
be varied, a short portion of a ring-slit 1:-25cm. long, corresponding to 8° of 
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Figure 5. Measurements of the F line of thorium B with the short preliminary slit B1B2. The 
values of the three parameters, viz. line-width at half-height, A(Hp)/Hp,  slit-width 
B1B2 (mm.), and source diameter (mm.), in that order are as follows for the four lines : 
(A058 2:2 552:059 (0) M75; 2:2; 1-0; (ce) 1/143; 0°95: 1-0: (d) 1/200; 0-95; 0-34. 
A solid angle of about 2:°5% of 4m was found for line (6) when normalized to correspond to 
collection over 290° of azimuth instead of the 8° accepted by the short slit. 


azimuth, was machined to a radius of 9cm. and provided with a radial sliding 
motion and a driving rod, shown in Figure 4, passing through a rubber vacuum 
gland. Wing-like baffles b3 and b4 were attached to the slit and moved with it. 
They did the work of B3 and B4 in Figure 1 by defining the width A« of the 
accepted sheaf and excluding rays of higher momentum. 
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An advantage of a short slit is that the image it forms is not sensitive to small 
deviations from axial symmetry so that it should give an upper limit of attainable 
resolving power. 

The correct radius of B1B2 was found to be 1 or 2% larger than the predicted 
value. Some shapes of the F-line of thorium B at 0-148 Mev. are shown in 
Figure 5. Fairly high resolutions are possible with a disc source which is 
sufficiently small. 

The preliminary short slit has now been replaced by a nearly complete 
annular slit built up of 18 movable 20° sectors possessing a concentric radial 
motion of about 3mm. to permit adjustment. f-ray lines found with this 
annulus are shown in Figure 6. 
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Figure 6. -ray lines measured with the complete adjustable ring-slit B1B2 set at a radius of 9:3 cm. 
anda width 0:7 mm. Line (a) had this slit in the mid-plane with B3B4 of width 1 cm. and 
radius 7cm. In (}) and (c) B1B2 was displaced from the mid-plane 5 mm. towards the 
image. In (c) only 100° of azimuth was used and B3B4 was narrowed to 3 mm. at a radius 
of 4:5cm. In these lines the source diameter was 0:34 mm. and the y-ray background 
was negligible. The values of A(Hp)/Hp and the solid angle x 1/47 normalized to 290° 
of azimuth were: (a) 1/104; 0:8%; (6) 1/113; 1:36%; (c) 1/140; 0-75 %. 

Line (d). The H line (1925 Hp) of radium B, taken with a thin radon tube lying along 
the axis, of diameter 0:25 mm. and length 8 mm. 


The large solid angle of collection makes it possible to work with only a 
fraction of the 360° range of azimuth ¢ when broadening due to axial asymmetry 
is suspected. Figure 6(c) shows a resolution of 1 in 140 obtained using a range 
of ¢ of about 100° selected by closing part of the smaller ring-slit B3B4. 


§9. PROCEDURE FOR CENTRING AND ADJUSTMENT 


‘The alignment of the pole-pieces was verified by passing pointed test rods 
through }-in. test holes bored accurately along their axes. The centring of the 
Geiger tube and other parts was checked by these rods or by sighting through the 
4-in. holes in the pole-pieces. ‘The axial symmetry of the field was tested by 
swinging a search coil around the axis. 

A striking feature of the focusing is the narrowness of the beam which enters 
the receiving slit S. This can be shown by moving the source along the axis and 
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finding the variation in counting rate as in Figure 7(a). For small axial motions 
Az of the source we can expect an equal but opposite displacement — Az of the 


_ image, if the ring slit B1B2 is in the mid-plane. This behaviour is a consequence 


of reflection which makes the system ‘ contra-current’. The curves show that 
the counting rate from a disc source is appreciable over a region Az of length 
about 0:‘5mm. If we subtract from this the width 0-35 mm. of the collecting 
slit S we see that the focal ring at S must be very narrow. The curves were 
obtained with slit widths of 0-7mm. at B1B2 and 6mm. at B3B4, using a flat | 
part of a continuous f-spectrum to ensure that the number of incident rays 
entering the ring B1B2 did not vary with source position. The beam at S 
broadens markedly if the centring is poor. 

It is clearly necessary to adjust the axial position of the source with care. This 
may be done by moving it longitudinally by a screw, passing through a vacuum 
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Figure 7. Longitudinal adjustment curves, showing the thinness of the sheaf near the 0:35 mm. 
wide collecting-slit S. (a) Axial displacement of the source alone, reckoned positive towards 
the counter. (b) Displacement en bloc of source, counter and slits. 

The tail at the right of curve (6) is duc to a small ‘ stray’ radioactivity on the sides of 
the mm. platinum wire, whose flat polished end received the active deposit from thoron. 
In (a) the deposition occuired through a masking hole of $ mm. diameter, on to a larger disc. 


gland, until the peak counting rate is found. A simple alternative is to move 
the whole vacuum tank axially between the magnet poles. A curve found in this 
way is shown in Figure 7 (6) and is narrower than curve (a), presumably because 
the previously mentioned ‘ contra-current ’ motion of the image is now combined 
with an opposite ‘ concurrent’ motion of the counter slit. 

The width and triangular shape of these curves suggest that most of the 
incident beam is spread evenly over a width Az nearly equal to that of the 
collecting slit S. This width of 0:35 mm. is nearly equal to the diameter of the 
disc source, and so is in rough agreement with the estimated image diffuseness. 


§10. RESOLVING POWER WITH AN AXIAL LINE SOURCE 


The remarkable result that rays are only counted if they come from a region of 
the axis within about 0-6 mm. of the object point O allows fairly good resolution 
to be obtained from a long line source along the axis. ‘The source is effectively 
shortened by the selecting action of the slits. The shape of the 0-264 Mev. H-line 
of radium B is shown in Figure 6(d).- This was obtained with a glass radon tube 
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of about 3-5mg/cm? wall thickness. This thickness of glass should give a 
standard straggling curve with a width at half-height of about Hp/90, according 
to White and Millington (1928). This accounts for only 10°% of the width of the 
line, to which it contributes as a statistically independent source of variance. 


§11. FOCUSING WITH LESS THAN UNIT MAGNIFICATION 
The lines shown in Figures 6(4), (c) and (d) were obtained with the ring-slit 
B1B2 displaced from the mid-plane about 5 mm. towards the counter. In sucha 
case the angle of emission «, at O will be less than the angle «, at entry into the 
counting slit S and the focused sheaf will no longer fall normally on the mid-plane. 
Such an unsymmetrical arrangement makes it possible to study rays coming from 
the source at smaller values of «, than 80°. 


§12. TESTS OF THE SOLID ANGLE OF COLLECTION 
Experimental values for the solid angle are attached to some of the curves 
in Figures 5 and 6. Each value assumes that the rays are collected over a range of 
azimuth of 290°, the full limit of the Geiger tube. For the ‘short’ slit in 
Figure 5 the absolute source strength was found by counting the «-rays with a 
spark counter kindly made available by Dr. R. D. Connor. In the other cases 

they were counted with a proportional counter insensitive to f-rays. 
The values given are probably lower limits because the slit S was rather 
narrower than was justified by the aberrations and errors in centring. In sucha 


case widening the slits would at first increase the solid angle without much loss 
in resolution. 
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ABSTRACT. The exposure was made at the Jungfraujoch, altitude 3,457 m., using 
Ilford type C2 nuclear research emulsions. The frequency of stars under thicknesses of 
carbon up to 220 gm/cm? was measured, the experimental value of the absorption length of 
the radiation causing them being 166+8 gm/cm?. The proportion of the stars due to fast 
m-mesons has been evaluated and, on subtracting this proportion, the absorption length 
of the radiation responsible for the remainder is found to be 143-+-10 gm/cm?. The energy 
distribution of the stars having from 3 to 20 tracks is found to be of the form dE/E?**, with 
no apparent discontinuity. ‘The frequencies of slow 7- and u-mesons at various depths is 
reported. ‘The z-mesons exhibit a transition effect which is interpreted in terms of their 
production in stars occurring in the carbon absorber. ‘The ratio of negative to positive 
7-mesons is found to be 3-1 to 1. 


$1. INTRODUCTION 

N a recent paper (George and Jason 1949a) an account was given of 

experiments with Nuclear Research photographic emulsions exposed under 

lead absorbers to the cosmic radiation at the Jungfraujoch at a height of 3,457 m. 
In the present paper we describe similar observations on emulsions exposed 
under carbon absorbers, also at the Jungfraujoch. For brevity we shall continue 
to refer to the strongly interacting component of the cosmic radiation causing the 
disintegration ‘stars’ as the ‘ N-radiation ’. 

A comparison of the absorption of the N-radiation in lead and in air had led us 
to the conclusion that the N-radiation consisted of nucleons. From the cloud 
chamber observations of Blackett, W. M. Powell and Hazen, it was thought that 
protons formed only a small part of the N-radiation, and that most of the stars 
were therefore produced by cosmic-ray neutrons. A similar conclusion was also 
published by Bernardini, Cortini and Manfredini (1949). 

Perkins (1947) and Bielovitski et al. (1949), on the other hand, have suggested 
that the stronger absorption of the N-radiation in air than in lead might be due to 
the fact that the particles of the N-radiation were unstable. ‘This would then give 
an absorption anomaly similar to that reported by Rossi, Hilberry and Hoag (1940) 
for the u-mesons. 

The present investigation was undertaken in order to decide between these 
viewpoints by comparing the absorption of the N-radiation in air and carbon. 
The preliminary results, consistent with the nucleon hypothesis, were communi- 
cated to the Cosmic-Ray Conference at Bristol. A more detailed account of the 
observations is given here, together with a discussion of the transition effects in 


carbon observed for the z- and the y-mesons. 
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§2. EXPERIMENTAL 


The absorber was constructed of bricks of compressed coal dust, built up in 
three cylindrical sections in the manner indicated in Figure 1. ‘The plates were 
exposed with the plane of the emulsion vertical and were placed at various depths 
in the absorber at intervals of about 50 gm/cm?, asshown. An additional packet of 
plates at the bottom of the absorber was screened with a further 75cm. Pb. 
The shape of the pile was chosen so as to give, as nearly as possible, a constant path 
length in the absorber for N-radiation incident at angles less than 45° from the 
vertical. The errors introduced by the absorber geometry at angles greater than 
45° are discussed below. The pile was covered with roofing felt to prevent water 
from melting snow dripping into the cracks between the bricks. Ilford Nuclear 
Research plates, type C2, 100 microns thick were used. Each packet was supplied 
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1. Observations: Z,=173 gm/cm? 


2. After subtraction of m1-meson 
contribution Z,=149 gm/cm? 


Frequency of Occurrence of Stars(No/cm?/day) 


pane sy W/), °b 0 100 200 300 
0 50. 100m. Thickness of Absorber (gm/cm?) 
Figure 1. Arrangement of carbon absorber. Figure 2. Frequency of ‘stars’ plotted 


against absorber thickness. 


with an electrical heater, and, using a distant reading thermometer, the temperature 
of the plates was maintained between the limits of 10° and 25° c. In order to 
minimize the effect of the background events accumulated while they were not in 
the absorber, the plates were taken to Switzerland on the day of manufacture, and 
were processed at the Jungfraujoch immediately on removal from the absorber. 
Corrections for the background, amounting to a few per cent, have been applied 
where appropriate. 

The absorber was built on the edge of the main roof of the Scientific Station 
in order to minimize the screening effect of local rocks. The roof beneath 
consisted of granite and concrete and was about 60cm. thick. Throughout the 
experiment the absorber was kept clear of snow, which fell frequently. 
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A sample of the coal was analysed and the results are shown in Table 1. 


Table 1. Analysis of Coal Absorber 


Constituent H C N O S Moisture Ash 
% by weight 3°9 80-9 1:9 0-2 2:0 1:3 9:8 


Although impurities constitute 19° of the total weight, fortunately they consist 
mostly of substances with atomic weights close to that of carbon. The absorption 
length L, in gm/cm? of the N-radiation in the coal would be 4% greater than the 
value for pure carbon, assuming L, oc A"?, Our final result, therefore, will be 
reduced by 4% in order to give the value of Ly for carbon. 


§3. EXPERIMENTAL RESULTS CONCERNING THE NUCLEAR 
DISINTEGRATIONS 


The plates were examined microscopically, using 4mm. objectives. Nuclear 
disintegrations giving rise to ‘stars’ of three or more tracks, of which at least 
one was longer than 60 microns, were recorded. In this section we give the 
experimental results concerning (i) the frequency of occurrence of stars at various 
depths in the absorber, and (ii) the size-frequency distribution of the stars. 


(1) Frequency of Occurrence of Stars as Function of Depth 


The frequencies of occurrence of stars at various depths in the absorber are 
given in Table 2. ‘The stars produced by 7-particles at the ends of their ranges 
are not included in this Table. The upper set of points in Figure 2 shows the 
logarithm of the star frequency plotted against thickness of absorber. 


Table 2 
Absorber thickness 220'E 
Ts (cexiom?) 0 55 110 165 220 +85 Pb 
Number of stars 1837 417 652 258 415 533 
Frequency per se? Mor 8:3 6°53 4°35 Beo7 
cm3/day +0-35 +0:56 +0:33 +0-41 +0:21 +0:15 


In plotting the last point, we have assumed the A’? law connecting absorption 
length with atomic weight. On this basis 85gm/cm*Pb is equivalent to 
33 gm/cm?C. From Figure 2 we deduce the absorption length of the N-radiation 
in the material of the absorber to be L, =173+8gm/cm?. This must be reduced 
by 4% to give the absorption length for carbon, with the result 


Ly, =166 + 8gm/cm2. ee Hh) 


(ii) The Size-Frequency Distribution of the Stars 


In previous communications (George and Jason 1949 a, b) we have plotted the 
logarithm of the frequency of occurrence of stars with n or more tracks, /( >”), 
against m the number of tracks. The results appeared to lie on the two straight 


lines 
I( >n) =I, exp (—0:48n), 3<n<8 ] 


and I( >n) =0-28 Ip exp (—0-33n), n>8. 


A similar result has been published by Page (1950). ‘Taking the stars at all depths 
together, the present results may also be expressed in the same form. 


ea) 
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The reason for the result (2) was never very clear, and we decided to plot the 
frequency of occurrence of stars against the energy released EF, rather than against 
the number of tracks”. In order to do this we obviously need a relation connecting 
E with n. From observations on stars produced in ‘ electron sensitive ’ plates, 
a formula has been given by Brown et al. (1949b) of the form 


E( Mev.) 3/0247" 0). eae (3) 


However, before using a formula of the type (3) to convert from an m scale to an EF 
scale, two modifications have to be introduced for the following reasons: (i) due to 


Number of Stars per MevV. Interval 


2 3 
R (units of 140 gm/cm? Carbon) 

ed ied ee & i ; 
i00 a 500 1000 2000 5000 10000 Figure 4. ‘The dependence of [N(>R)]-/? on R. 


Energy of Star (mev.) N(>R) is the integral range spectrum of 

m-mesons emitied from stars. The experi- 

Figure 3. Differential energy distribution of mental points are those of Camerini et al. 
stars. The full line is a plot of 14+ R/b, with b=2-12. 


the limited sensitivity of the C2 plates, the value of observed in our plates will be 
smaller than that in the plates of Brown et al.; (ii) as stated by Brown ez al., the 
energy total of the ‘shower’ particles has not been included by them. 

Taking into account these two effects, we find thatthe relation between Eandn 
in our plates may be represented approximately by 


E(Mev.) =26n + 12-8n?. Roce ool a 
The details of the argument leading to (4) are given in the Appendix. Figure 3 
shows the differential energy distribution of the stars, plotted on a logarithmic 
scale. ‘he energy has been obtained from the star size, using the relation 
expressed approximately by equation (4). It is seen that over a wide range of 
energy the differential energy distribution may be represented by the relation 

I(E)dE « dE/E”, wa taten(o) 

where y =2:6+ 01. 
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§4. DISCUSSION OF RESULTS CONCERNING NUCLEAR 
DISINTEGRATIONS 


(i) The Absorption of the N-Radiation 
From the absorption length of the N-radiation. in air, L, =150+4gm/cm? 
{George and Jason 1949 b) and, assuming L, « A1/3, we deduce that the value of Ly 


for carbon should be L, =141 gm/cm?. The result of the present investigation, 
166 + 8 gm/cm®, appears to be significantly greater than the expected value. ‘Thus 


_ we may conclude either thatthe assumption L, « A’ is not valid in going from air 


to carbon, or that there is a small absorption anomaly due to unstable particles that 
produce stars in the carbon absorber but decay before doing so in the air. The 
extrapolation of the expected value of L, is hardly likely to be in error over such a 


/ small range of atomic number and it therefore appears that the presence of unstable 


particles must be assumed. 

It seems reasonable to attribute this excess of stars in carbon to the fast 7-mesons 
emitted in the more energetic nuclear disintegrations (Brown et al.1949a). These 
are known to produce stars with an efficiency comparable with that of protons. 
(Camerini et al. 1950), and hence must make some contribution to the star 
population observed in the carbon absorber. We now propose to evaluate this 
contribution approximately. 

To do this we shall need the range spectrum of the emitted mesons. Sands 
(1950) gives this for mesons produced in the atmosphere in the form 

b’ ; 
N(>R)= (+Ry with y=1-9. ere (6) 

‘Camerini et al. (1950) have shown that there is no appreciable difference in the 
energy spectra of the mesons produced in stars with different numbers of heavy 
branches, and that therefore the energy spectrum is, within the limits of measure- 
ment, independent of the nature of the nucleus in which the mesons are created. 
Also, they find good agreement with the spectrum given by Sands for the mesons 
producedinair. We may therefore take it that the energy spectrum of the mesons 
emitted in stars produced in the carbon absorber is the same as that found for stars 
produced in nuclear emulsion. 

For the purpose in hand, i.e. the evaluation of a correction, we have assumed for 
simplicity y=2. Then we have attempted to fit an expression of the form (6) with 
y =2, to the data concerning the energy spectrum of 7-mesons emitted from stars 
kindly supplied to us by Camerini et al. (1950), now published by them. For such 
a relation we should have, simply, 


[M(>R)P?2=14+ as ae (7) 


In Figure 4 we have plotted the empirical values of [M(>R)]|-1? against R and it is 
seen that they may be reasonably well represented by an expression of the form (7). 
For later convenience we have chosen to take 140 gm/cm* as the unit of length in 
the carbon pile. Expressing R in terms of this unit, b has the value 


6=2-12. Seaeate(S) 


We shall also need the angular distribution of the emitted mesons, but as yet no 
complete data are available and we have to proceed by assumption. One would 
expect the low energy mesons to be emitted almost isotropically, and those of high 
energy to be emitted at small angles to the incident particles, themselves strongly 
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collimated in the vertical direction. This would lead to a gradual change of the 
angular distribution with increasing energy of emitted mesons from the isotropic 
to the vertically collimated distribution. 

In the next section we show that from the observations on the frequency of 
slow 7-mesons stopping in the emulsions, approximately one-half of the emitted 
mesons fall into the low energy isotropic class. In order to display the main 
features of the argument we therefore make the assumption that mesons emitted 
with energy below the median value of their energy spectrum (270 Mev.) are emitted 
isotropically, while those with energy greater than 270 Mev. are collimated in the 
vertical direction. This is obviously a crude model, but it will serve to give us a 
reasonable approximation to the effects due to 7-mesons on the star absorption 
curve. 

Camerini et al. (1950) have shown that the interaction length of 7-mesons of 
energy less than 1,500 Mev. is close to that of protons of energy less than 600 Mev. 
The mean free path of neutrons will presumably be the same as that of protons. 
Hence, in what follows, we assume that the mean free path of the 7-mesons in the 
carbon absorber is the same as that of the stable N-radiation. A more refined — 
calculation in which the difference in path length of mesons and nucleons is 
included is hardly justified at present. 

We now calculate the contribution to the star population caused by 7-mesons 
in the carbon absorber, neglecting the loss due to decay, which is negligible 
under the conditions of our experiment. The median range from (7) is given by 
Ry =0-4146=0:88=123 gm/cm*. In view of the unit of length chosen 
(140 gm/cm?), the frequency of stars at depth x produced by the stable component 
of the N-radiation may be taken as Ae~“ and for the purpose of this calculation we 
may put A=1. Let a be the number of mesons emitted per star. Those 
emitted with ionization range less than Ry are assumed to be isotropic, and the 
number of stars they produce at depth x will be given by 


Sx) = 2 ann Sasi Reost+ R)}2nsind| (5-72) =0-5 | aRao, 


For x> Rp, this becomes 


S,(x) =0-5 ae-# ir (Ge) (eer) -05 | dR 


or S1(%) =0°-27 ae-* (a> Re)ask sone gees (9) 


For x«<R, the lower limit in the integration over @ is modified and the integral for 
S\(x) has been evaluated by numerical methods. We may express the result in 
general as 

9 (*) Ge ay = eee (10) 


and values of F(x) will be given after we have discussed the star contribution of 
the remainder of the emitted mesons. 

For those mesons emitted with range greater than Ry, we are assuming that they 
are collimated in the vertical eae Se Their integral range spectrum may 
therefore be written 


f(R) =0°5 for R<R, | 
. seco (te 


b 2 
x (R) for R>R,. | 


M 


i 


——— 
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‘The star frequency at x due to all of these collimated mesons produced above x is 


So(x) =ae-* | “ f(R)dR 
0 
or So(x) =0-5axe-* PEI | 


Sc(x) =0-5ae-# | Ro+ aierad x>Ry. 


Let us express this result as 
eer Yee el se ie (13) 
In order to be able to derive numerical values from equations (9), (10) and (13) 


__we shall need to know the value of a for carbon. The value for nuclear emulsion 


may be derivedas follows: Table III of Brown et al. (1949 b) shows that the number 
of ‘shower’ particles per star is 0-25, while, according to Fowler (1950), at least 90% 
of the ‘shower’ particles are 7-mesons. Hence there are at least 0-22 7-mesons 
per star emitted with energies greater than 70 Mev. ‘To this we must add 0-05, 
representing the number of mesons emitted with energies less than 70 Mev. (Fowler 
1950), giving a total of 0:27 7-mesons per star. ‘This will err on the low side, but 
only by a small amount, since we have ignored those mesons emitted in stars with 
0, 1 or 2 heavy branches. Thus we may take a =0-27 for nuclear emulsion at the 


Table 3 
ip S(x) 
FF F ——<—<—<—<$<<$<$<$<$< 
cfs eee EB) EONS rata ake} 
0) (@) 0:033 (0) 14:7 +0-35 
55 0-393 0:073 0:052 10:1 +0-49 
110 0:786 - 0-076 0-106 7:0 +0:28 
165 1:18 0-080 0-156 5-28+0:°32 
220 1°57 0:080 0:195 3:41+0:16 
253 1°81 0-080 0215 2:60+0-11 
(6) 0-080 0:52 


Jungfraujoch, and we believe that the same value may be taken for carbon for the 
following reasons. From the number of slow z-mesons stopping in our emulsions 
we are able to deduce that there are 0-15 mesons per star in carbon emitted with a 
range less than 120 gm/cm? (cf. Harding and Perkins 1949 and §6 (1) below). 
The centre packet of plates is screened by this amount of carbon in all directions. 
From the fact that the energy spectrum of the mesons produced is substantially the 
same in nuclear emulsion and carbon, and that 120 gm/cm? is the median value of 
the range spectrum (7), it follows therefore that the total number of emitted mesons 
is approximately twice the number deduced from those stopping, that is, the total 
number of emitted mesons in carbon must be close to 0:3 per star. We have chosen 
to keep to the measured value of a equal to 0:27. 

In order to compare the results under carbon absorbers with those found for 
absorption in air, we must obviously eliminate the 7-meson contribution. This is 
done by dividing the observed star frequencies by 1+a[F,(x)+ F,(x)] with 
a=0:27. The numerical values are given in Table 3. 

The results corrected in this way are given in the last column of Table 3, 
and are plotted in Figure 2. The absorption length of the stars after elimination 
of the meson contribution is found to be 149 + 10 gm/cm? coal, corresponding to 
143 + 10 gm/cm? carbon. 
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It therefore appears that if one subtracts the stars attributable to 7-mesons: 
generated in the carbon absorber, the absorption length of the radiation responsible 
for the remainder is very close to the value expected from the absorption of the 
N-radiation in air. This radiation may be taken to consist of nucleons. In 
dense absorbers there is an additional component of the N-radiation, the 7-mesons. 
In absorbers of thickness up to several hundred gm/cm? at the Jungfraujoch, 


these 7-mesons may contribute up to 30°, of the observed star frequencies, and — 


this contribution does not vary exponentially with depth. The resulting curve 
of the frequency of stars due to both nucleons and 7-mesons together plotted 
against absorber thickness does not however depart very strongly from an 
exponential curve, but the absorption length deduced from this quasi-exponential. 
will of course be greater than that for nucleons alone. In the atmosphere the 
m-mesons decay, in general, before causing further stars. 

The above deals with the major contributions to the star population i in dense 
absorbers of moderate thickness, but is by no means the whole story, since we 
know that the p-mesons also make a small contribution which under sufficiently 
thick absorbers may be quite important (George and Evans 1950). We estimate 
that the p-meson contribution is approximately 0-1 stars/cm®/day at the 
Jungfraujoch, which means that, because of their small absorption, they make the 
dominant contribution under carbon absorbers more than 700 gm/cm? thick. Atthis 
thickness of absorber, therefore, the departure from the exponential absorption of 
the N-radiation should be quite marked, the star frequency thereafter decreasing 
only in the slow manner characteristic of the hard component. The results of 
Barford et al. (private communication) on the star frequency under large thicknesses 
of ice show a departure from the exponential curve at this depth which might be 
explained by this -meson contribution. 

The possibility must be considered that stars may be caused by particles of the 
cosmic radiation other than the nucleons, 7-mesons and x.-mesons discussed above, 


for instance, photons, heavier mesons, etc. It would appear, however, that the — 


contributions of these additional particles must be small, since the main features of 
the star absorption curve at all depths up to 6,000 gm/cm? can be understood 
without invoking them. 

It has been found by Cortini and Manfredini (1949) and by Malaspina et al. 
(1950) that a maximum occurs in the transition curve for stars of three and four 
heavy branches at approximately 1cm. Pb and approximately 2cm. Al. No 
maximum is found for larger stars. ‘This result has been interpreted by Dallaporta, 
Merlin and Puppi (1950) in terms of secondary neutrons generated in evaporation 
processes in the condensed material surrounding the plates. It is clear that the 
m-mesons discussed above must make some contribution to the frequency of the 
stars under these thicknesses of absorber, but on evaluating it from equations (9) 
and (13), the contribution is found to be much smaller than the magnitude of these 
transition maxima found experimentally. ‘Thus the transition effect due to 
m-mesons discussed above is different from, and should not be confused with, the 
transition effect observed in small stars under small thicknesses of absorber. 


(it) The Size—Frequency Distribution of the Stars 


It has been seen that when we convert the size distribution to an energy 
distribution, the resulting differential energy distribution of the stars may be 
represented by a power law of the form N(E)dE a dE/E” with y=2-6+0-1, for 


a. 
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| j values of EF in the range 300 to 6,000 Mev., corresponding to stars with from 


three to twenty' branches in C2 emulsions. - We note that a similar value for the 
exponent in the energy distribution of the 4-mesons has been found by Wilson 
(1950). The use of equation (4) may have caused some fortuitous smoothing 
of the energy spectrum due to fluctuations, but this is hardly likely to account for 
the linearity of Figure 3 over such a wide range. 

A simple interpretation of the observed energy spectrum of stars is that it 
is a reflection of the energy spectrum of the N-radiation, which in its turn is deter- 
mined by the energy spectrum of the primary radiation. These would all be power 
spectra with the same exponent if the energy loss £, in a nuclear collision of an 
incident particle of energy Ej, were a homogeneous function of E/E). This 
latter assumption is one very frequently made in the literature (Clementel 1950, 
Heitler and Janossy 1949 a,b). Granted its validity, then a primary energy spectrum 
of the form (5) would give rise to a similar distribution in the star energies and 
(apart from effects due to decay) in the energy spectra of the 7- and y-mesons. 
This appears to be in agreement with the energy spectrum of the protons ejected 
from stars measured by Camerini et al. (1950), who have shown that the differential 
spectrum follows a power law of the form dE/E?® for values of the kinetic energy of 
the order of 1,000 Mev. 

The fact that the size distribution in C2 plates may be represented as the sum of 
two exponentials must be regarded as fortuitous. The interpretation in terms of 
an energy spectrum would appear to be far more logical, and in any case easier to 
understand. Other interpretations are not excluded, but we have none to offer. 


§5. EXPERIMENTAL RESULTS CONCERNING SLOW MESONS 


Mesons stopping in the emulsion were identified by means of their characteristic 
scattering and variation of grain density with range. hey were classified according 
to their behaviour, following the scheme initiated by Lattes, Occhialini and Powell, 
ritO 


(i) p-mesons: mesons stopping in the emulsion and not accompanied by 
tracks of secondary particles. 

(ii) c-mesons: mesons stopping in the emulsion, accompanied by nuclear 
disintegrations at the ends of their ranges. 

(iii) w-mesons: mesons stopping in the emulsions, accompanied by the track 
of a secondary meson (u-meson) commencing at the end of the range of 
the first particle. 

(iv) p-mesons: the secondary particles arising from 7-mesons above. 


It is known that only 71°% of negative 7-mesons stopping in photographic 
nuclear emulsions give rise to recognizable ‘stars’; the remainder are presumably 
captured by atomic nuclei with the resultant emission of neutrons only. Hence 
some 77~-mesons will be included in the category of p-mesons. 


(1) p-Mesons 


The p-mesons stopping in the emulsion and possessing a projected range greater 
than 60 microns were divided into two groups according to whether they entered 
the emulsion from above or below the horizontal. Then, using the geometric 
factor given by Lattes, Occhialini and Powell (1947), it follows that only 0-56 of the 
p-mesons stopping in the emulsion would be recorded. The observed numbers 
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were accordingly divided by 0-56 in order to give the total numbers of p-mesons 
stopping in the emulsion. The results at various depths in the absorber obtained 
in this way are given in Table 4. 

The frequency of downward p-mesons at T=0, 4:15 + 0-25/cm?/day, is in good 
agreement with the value 4-1/cm3/day reported by Camerini, Muirhead, Powell and 
Ritson (1948), also obtained using nuclear emulsions exposed at the Jungfraujoch. 


Table 4 
Thickness, 7’, of absorber x ! 220 C 
(emjenhe) 0 55 110 165 220 85 Pb 
‘Total p-mesons 4-15 4-18 3-95 4-80 3-60 3-86 
downward/cm?/day +0:25 +0-43 +0:-28 +0°47 +0:25 +0-21 
‘Total p-mesons O54: 0-46 0-33 0-46 0:34 0-26 
upward/cm3/day +0-09 +0:-14 +0-07 +0-14 +0-07 +0-05 


These results may be compared with the number expected from the observations of 
various investigators using systems of Geiger counters. ‘Thus ifz, is the number 
of mesons stopping per gramme of air per sterad persecond in the vertical direction, 
under very thin absorbers, then since the slow meson intensity varies as cos*?6 
(Kraushaar 1949), and taking 4-0 as the density of C2 nuclear emulsion, we obtain 
the following expression for the expected number of u-mesons stopping per cm? 
per day in the emulsion: 


gee i 670 LOS at a3 201 Ue ees (14) 
(the 1-5 in the denominator is the relative stopping power of air to emulsion). The 
earlier measurements of Rossi, Sands and Sard (1947) give 

7, =1-8 x.10-> gm/sec/sterad 
at a height corresponding to the Jungfraujoch, while the more extensive 
measurements of Kraushaar give 7, =1-3 x 10-° gm/sec/sterad. 

Using (14) these results give 5-8 and 4-2 respectively for the expected number of 
4-mesons stopping per cm? per day in the emulsion and these are seen to be in good 
agreement with the observed number of p-mesons. 

The differential range spectrum of the downward stream of p-mesons is seen 
from Table 4 to be constant in the range 0 to 265 gm/cm? (for ionization loss 
85 gm/cm? Pb =45 gm/cm?C). By way of illustration, we assumed the frequency 
of mesons stopping to be the same at all thicknesses of the absorber, and then 
applied the x” test to the observations to see whether they were consistent with 
this assumption. The best fit was obtained for an assumed constant value 
of 3-9 p-mesons/cm*/day, and the corresponding value of x? is 6. According 
to Fisher’s tables, this value of y? would be expected in 40% of the cases as a result 
of random fluctuations. 

There are no measurements using nuclear emulsions with which our p-meson 
results under absorbers may be compared. Using systems of counters, the 
differential range spectrum has been shown to be flat up to 100 gm/cm? air 
equivalent at Echo Lake by Kraushaar (1949), and a similar result at several points 
in the atmosphere was obtained by Rossi, Sands and Sard. We may note that if 
the differential range spectrum at the Jungfraujoch is at all similar to that at sea 
level, then it should remain fairly flat up to about 1,000 gm/cm?, after which it 


z 
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should fall rapidly (Rossi 1948). Our maximum absorber was considerably less 
than 1,000 gm/cm?. 

The upward stream of p-mesons has been attributed by Camerini, Muirhead, 
Powell and Ritson (1948) to the decay in flight of the upward stream of 7-mesons. 
We may note here that the frequency of upward p-mesons is of an order of 
magnitude comparable with that of the z-mesons (Table 5). | 

We will leave the discussion of the upward p-meson results until we have 
considered the results for 7-mesons. 


(it) a-Mesons 

The frequencies of o-mesons giving rise to stars, N(c), and those showing 
p-decay, N(z, ~), were noted, as well as their directions of motion on entering the 
emulsion. In order to determine the total numbers of a+- and z-particles 
stopping in the emulsions, the observed frequencies were corrected for the loss due 
to the geometry of the emulsion. A simple extension of the method of Lattes, 
Occhialini and Powell (1947) shows that 18-69% of the positive and 11:3% of 
the negative mesons would not be recorded for reasons of geometry. Furthermore, 
Adelman and Jones (1949) have shown that only 71°% of 7-mesons give rise to 
nuclear disintegrations. Hence the actual frequencies of 7-- and 7+-particles 
were obtained as follows: 


N(a*) =1-23 N(x, 1), 
N(a-) =(1-13/0-71) N(o) =1-6 N(o). 
The total frequencies of 7-mesons stopping in the emulsion at various depths in 
the absorber, obtained with the aid of equation (15), are given in Table 5. 
In Table 5, z-mesons of both signs of charge have been added in order to 


give the total frequencies entering the emulsion in the downward and upward 
directions. 


Table 5 
‘Thickness, T, of absorber é 220KS 
Goin) 0) 55 110 165 220 +85 Pb 
a-mesons down/cm3/day 0-36 0-72 0:77 0-73 0:56 0-36 
+0:-05 +0714 +0:09 +0-13 == O20 Seem 0-05 
a-mesons up/cm?/day 0-86 0:80 0-54 0-55 0-47 0-23 
+0:08 -+0-15 +0-08 +0-:12 +0°07 =+0-04 


The results are plotted on a logarithmic scale in Figure 5. "The downward stream 

of 7-mesons appears to exhibit a transition effect, passing through a flat maximum 
in the vicinity of 100 gm/cm? (cf. Harding and Perkins 1949). ‘The results for the 
upward stream of 7-mesons may be fitted to an exponential curve, Figure 5 (4) ; 
the line drawn has the same slope as the line in Figure 2, giving the variation of 
nuclear disintegrations with depth. We show below that the frequency of upward 
a-mesons would be expected to be proportional to the frequency of nuclear 
disintegrations, and it is seen that, bearing in mind the limited accuracy of the 
results, the observed points fit the line reasonably well. 

The total frequency of z-mesons for T=0 is 1-22 + 0-1 and is consistent with 
the figure of 1-42 + 0-16 given by Camerini, Muirhead, Powell and Ritson (1948) 
for unscreened plates exposed at the Jungfraujoch. 
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Another feature of the results that may be obtained is the ratio of negative to 
positive 7-mesons. In Table 6, we give the results obtained by disregarding the 
direction of motion of the 7-mesons, and taking simply the ratio of N(w~) to N(7*). 

It will be seen that the ratio of 7-- to 7+-mesons remains substantially constant 
with depth in the absorber, and the mean value, averaged over all depths, is found 
to be 

N(x) 
N(a*) 


The only result with which this may be compared is that of Camerini, Muirhead, 
Powell and Ritson (1948), who report that the values of N(c) and N(z, 4) are 


Har 25 oe Oe (16) 
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Figure 5. The frequency of occurrence of slow 


m-mesons. ‘The calculated vluaes are 
those given in ‘lable 5. 


Table 6 

Thickness, T, of absorber 220 C 
Cac C) 0) 55 110 165 220 85 Pb 
N(7-)/cm3/day 0°93 1:14 1:0 6:97 0:74 0-44 
+0-09 +0-18 +0:11 --0:16 +0:09 -+-0:05 

N(a*)/cm3/day 0-30 0-38 0-28 O30 0:29 0-15 

N(a-)/N(ar*) $81 3-0 3°6 Sl 2:6 3-0 

+0°6 +0-9 +0°8 a= 1-0 +0°6 +0:7 


0:74+0-11 and 0-68 +0-23/cm%/day respectively. This gives for the ratio of 
negative to positive mesons 


Nr) _ (0:74 + 0-11) 
N(x*) — 0-71(0-68 + 0-23) 


which does not seem consistent with the value (16). Examination of Table 2 of 
these authors reveals that they would appear to have used a much larger factor 
(~4) to correct for the loss of positive mesons than for the negative ones. In the 
light of present knowledge concerning 7-mesons, we believe that the values (15) are 
to be preferred. 


=1-440-5, 
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§6. DISCUSSION OF MESON RESULTS 
In the previous section it was seen that the results for the downward stream of 
p-mesons were consistent with the observations of other workers using nuclear 
emulsions and systems of Geiger counters. In this section we discuss the frequency 
of the downward and upward stream of 7-mesons, the upward stream of p-mesons 
and the negative/positive ratio of 7-mesons. 


(1) The Frequency of m-Mesons as a Function of Depth in the Absorber 


Examination of Table5 shows that the numbers of 7-mesons inthe upward and 
downward streams are approximately of equal magnitude, at least in the lower half of 
the absorber when the transition effects of the downward stream are in equilibrium. 
This forms the experimental justification of the assumption made in §4 that the 
lower energy 7-mesons are emitted isotropically from ‘stars’ (cf. Franzinetti 1950). 
The number upward is a little lower than the number downward, but this is to be 
expected since the frequency of stars below a given point in the absorber is lower 
than that occurring above. By carrying through a calculation similar to that of 
Harding and Perkins (1949), but including the upward as well as the downward 
stream of 7-mesons, it is found that the number of 7-mesons emitted isotropically, 
and with energies less than 270 Mev., and which come to rest in the emulsion, is 
0-15 per star. This figure is to be compared with the total number of mesons 
emitted, viz. 0:27 per star, mentioned in §4(i). The difference accounts for the 
m-mesons producing stars already discussed in §4(i1). We now show that the 
same assumptions and method of calculation which we used in §4(i) will also 
account for the observed frequencies of mesons stopping in the emulsions. 


(a) Contribution from collimated mesons. 
The differential range spectrum n(R) corresponding to equation (11) is for this 
group 


n(R)dR=0 forR <R, 
2 
: at 5 for R> Ry. | pee 17) 
The frequency of 7-mesons from this group stopping in dx at x is given by 
a—Ry 252 
N((7) = ae-*dx i, Oia_e) dx 
=ae~*dx[0-5—{b/(b+x)}"] x>R, | 
Ey fer i ere (18) 
The number of stars starting in dx is e~*dx and the ratio is 
GIO OO al ny sub ds Pukey (19) 


Hence to obtain the frequency relating to our conditions we multiply the observed 
frequency of stars at x, S(x), by (19) and by a factor D expressing the probability of 
surviving decay 

N¢(7) =DS(x)a[0-5 — {b/(6 + «)}"] percm’perday. ...... (20) 
These mesons will all be observed in a downward stream. Numerical values will’ 
be inserted after we have evaluated the other contributions. 


E3=2 
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(5) Contribution from isotropic mesons. — 


The frequency of mesons stopping per cm® per day arising from the isotropic 


group is : 
Down: 
bi —R b> ve Ros a3 
N,(mp) = DaS(«) | pepe | eben’ sin 0 db 
Ry b2 m/2 ; > 
+ DaS(x) | ope | e* 08? sin Od0 fot & a Re poner (20) 
+2 yp 
{cos$=x/R,) and 
rR b2 742 
N,()=Das(q) | par | elt cos) go forw Ry Oe ee (20) 
0 0 
Up: 
N,(my) =DaS(x) [e® "aR |" oP 4sino dd 
1(7y) = Da: («)| é (+ Ry oe sin 


= DaS(x) 0197 = ee (20c) 


The integration over 6 is straightforward and the remaining integration over R has 
been performed by numerical methods to give the results in Table 7. 


Table 7 
emfem® c) Nelma)lem* {day Ny(orp) —Ny(ore) Nom») +. NuCorp) 
0 0 0 0-86 0) 
14 0 0:47 0:77 0:47 
28 0 0-61 0-72 0-61 
55 0 0:75 0-62 0-75 
110 0 0-68 0-45 0-68 
165 0-16 0-49 0-33 0-65 
220 0-20 0-33 0-24 0-53 
220 C+85 Pb 0-28 0:22 0:21 0:50 


‘The value of D has been determined from the mean measured lifetime of the 
‘m-mesons, T=1:5 x 10-’sec. (Richardson 1948, Martinelli and Panofsky 1950), 
and a mean value of D =0-9 was used in evaluating the results in Table 7. 

The results in the last two columns of Table 7 are plotted as the dotted lines in 
Figures 5(a) and (0), and it is seen that the agreement between the calculated and 
observed values is reasonable. 

We must assume that on account of the low density of the atmosphere, and 
consequently the small value of D, the frequency of downward 7-mesons arriving 
from the atmosphere must be negligibly small. ‘Thus the expected frequency 
of downward mesons is zero at T=0, whereas the frequency observed was 
0:36 +0-05/cm3/day. The increase in downward mesons with thickness T is 
however very rapid, and at T=14gm/cm? the expected number is 0-47. The 
plates used were 3 in. x 4in., standing on their longer sides, and in actual fact 
the absorber thickness at the nominal value JT =0 should be taken as approximately 
10 gm/cm?, made up of the light covering over the plates, and of the glass and 
emulsion of the plates themselves. The observed frequency at the nominal 
value of T=0 can therefore be understood as being due to mesons generated 
in the local material, and not necessarily due to mesons arriving from the air. 

The calculated value of downward z-mesons is a little too high and that 
of the upward z-mesons is a little too low. This is probably due to the curved 
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@ shape of the top of the carbon absorber, and better agreement could be obtained 


_ by taking into account the geometry of the pile, but we have not thought it worth 
while to push our calculations to this degree of refinement. An approximate 
estimate shows that the remaining discrepancy can be removed by these 
considerations. 

Thus our crude model concerning the angular distribution of the 7-mesons is 
perhaps not too far fetched, in so far as, at one and the same time, it shows the 


__ expected number of slow z-mesons to be in agreement with the observed upward ° 


_ and downward streams, and serves to explain the small difference in the observed 
absorption lengths for the N-radiation in air and in dense material of comparable 
atomic weight. 

(11) The Upward Stream of p-Mesons 


There are four possible sources contributing to the upward stream of p-mesons: 
(a) 29°, of the upward stream of 7-particles do not give rise to stars, and hence 
will be recorded as p-mesons; (b) about 10° of the upward stream of 7*- and 
7 -particles will decay in flight, giving u-mesons ; (c) both downward and upward 
7~-particles coming to rest in the glass backing of the emulsion will give rise to 
p-mesons; these pz-mesons are emitted isotropically and hence 50% will coine 
to rest in the emulsion travelling in an upward direction ; (d) back-scattered 
p-mesons originally moving downward. 

The contribution from (d) is uncertain, but it must be small, since the sum of the 
contributions (a), (b) and (c) is sufficient to account for the greater part of the 
observed frequency of upward p-mesons, as will be seen from Table 8 where the 
individual contributions are tabulated. 


Table 8 

Observed 
T (gm/cm?) C (a) (5) (c) (a) +(6)+(c) frequency of 
0-29 N(wz). 0-1 N(ayg) = 05 N(ar*) up p-mesons 

O 0-19 0-09 0-15 0-43 0:53+0-09 

55 0-18 0-08 0-20 0-46 0:46+0-13 

110 0-12 0-05 0-13 0:30 0-32+0-07 

165 0-12 0-05 0-15 0:32 0:46+0-14 

220 0-09 0-04 0-15 0-28 0-34-+0-06 

220 C—85 Pb 0-05 0-02 0-08 0-15 0-26 +0-06 


Values are in number/cm?/day. 


(iii) The Ratio of Negative to Positive 7-Mesons 


It was found that there was a negative excess in the 7-mesons stopping in the 
emulsion at all depths in the absorber, the mean ratio of negative to positive mesons 
being 3-1+0-25. 

A negative excess in the mesons which are ejected from stars and come to the 
end of their range in the same piece of nuclear emulsion has been noted by C. F. 
Powell (1949). The mean energy of this group of mesons is 1-6 Mev., and this 
negative excess has been attributed to the effect of the nuclear Coulomb field. 
An explanation along similar lines would hardly hold in the case of our mesons 
since their mean energy is much higher, of the order of 100 Mev. 

A more likely explanation seems to be the following: The ratio of negative to 
positive 50 Mev. mesons produced by 345 Mev. protons in various elements has 
been studied at Berkeley by Bradner and Jones (1950) and by Richman and 
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Wilcox (1950). For carbon the ratio is 4:8 + 0:5 positive mesons io one negative. 
For incident neutrons, one would expect from the charge symmetry of nuclear 
forces that the situation would be completely reversed, since in the carbon nucleus 
there are equal numbers of protons and neutrons, i.e. one would expect 4:8 

‘ negative mesons for 1 positive meson. Thus if « is the proportion of cosmic-ray 
“ stars’ produced by protons, the expected ratio of negative to positive 7-mesons 
stopping in the emulsion would be 


N(w-) _ 48—3-8« 
Mat) 143-80 ° 


‘The observed ratio of 3-1 + 0-25 corresponds to a value for « of 0-11+0-04. This 
is similar to the proportion of charged particles in the N-radiation at the 
Jungfraujoch. Page (1950) gives « =0-17, and Barton (1950) finds 0-16. ‘These 
experimental figures are upper limits, since not all the charged particles in the 
N-radiation are protons: 7-mesons are probably also included in these observed 
values of «. In any case they are not inconsistent with the value « =0-11 + 0-04 
required above. 

One consequence of the above interpretation is that as the depth is increased 
the negative/positive ratio should decrease steadily from the value 3-1. This is 
due to the increasing proportion contributed by the more energetic collimated 
mesons to the total number stopping (cf. Table7). ‘The value of « for the energetic 
events is close to 0-5 (Brown et al. 1949), which from equation (21) is seen to result 
in a lower value of the charge asymmetry of the 7-mesons produced in carbon. 


ee 


§7. CONCLUSIONS 
The conclusions at which we have arrived, based on measurements at the 
Jungfraujoch, are as follows : 

(i) The measured value of the absorption length in carbon for the 
radiation producing ‘ stars’ is L, =166+8gm/cm?, using absorbers up 
to 250 gm/cm? in thickness. 

(11) In dense absorbers, the N-radiation is composed partly of 7z-mesons, 
whereas in air the 7-mesons are absent, due to the high probability of 7—u 
decay. 

(ii1) The presence of these 7-mesons, generated in the absorber, gives rise to an 
increase in the apparent absorption length of the N-radiation of about 14% 
in carbon, for the absorber thickness used in our experiments. 

(iv) The absorption length of the N-radiation in carbon, when corrected for the 
presence of these secondary stars due to z-mesons, is 143 + 10 gm/cm?, 
which is close to the value expected from the known absorption of the 
N-radiation in air. 

(v) From the size distribution of ‘ stars’ we conclude that their distribution in 
energy is of the form dE/E” with y=2-6+0-1, in the range of E=300 to 
6,000 Mev. 

(vi) The differential range spectrum of the downward stream of -mesons at 
the Jungfraujoch is constant up to 265 gm/cm?. 

(vii) The mean ratio of negative to positive 7-mesons with energies between 


0 and 450 Mev. produced in carbon by the cosmic radiation at the Jungfrau- 
joch is 3-1+0-25. 
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| _ (viii) About one-third of the upward stream of p- eee at the Jungfraujoch is 


composed of 7 -mesons that do not give rise to ‘stars’. Most of the 
remainder is due to the decay of 7+-mesons in the elas backing of the 
emulsions. 
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APPENDIX 
The Relation between Energy and Size of ‘ Stars’ observed in C2 Emulsions 


The following relation has been a by Brown et al. (1949b) for ‘electron 
sensitive’ nuclear emulsions: 


PMC VN STH TAN go we ie (22) 


E is the energy released among the nucleons, and n is the number of branches in 
‘stars ’ with a grain density greater than 1-54,,;,. We require the corresponding 
relation for C2 emulsions for the total energy released. ‘This may be done by using 
the data in Table IV, given by Brown et al., and taking into account (1) the limited 
sensitivity of the C2 emulsions, (ii) the energy released in the form of fast particles 
with grain density less than 1-5 9,,.... 


(i) Limited sensitivity of C2 emulsions. 
Brown et al. divide the heavy tracks into three groups : (a) short black tracks 
<70p, (6) long black tracks >70p, (c) grey tracks with grain density >1:5 Znin 


and <5 Zin 


In C2 emulsions all the tracks in groups (a) and (6) would be detected, but a 
fraction of group (c) would not. Taking the upper energy limit for protons of the 
C2 emulsion as 100 Mev. (Page 1950), we may determine the fraction of the (c) 
tracks lost from Figure 14 of Brown et al. as a function of the size of ‘ star’ as 
observed in ‘ electron sensitive ’ emulsions, and hence obtain the mean number of 
branches that the ‘ stars’ considered by Brown et al. would have had in C2 
emulsions. 


Table 9 
Number of heavy branches in 4 
“ electron sensitive ’ emulsions 4,5 T8 12-15 16-24 
Fraction of grey tracks recorded 
in C2 emulsions 0-86 0-73 0-62 0-63 


Mean number of branches ob- 
served in C2 emulsions 4.03 71 11-9 16:4 
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(ii) Shower particles. 


From Table III of Brown et al., we may determine the mean number of shower 
particles, 7,, as a function of star size (‘Table 10). 


Table 10 
Number of heavy branches 4,5 Thy «s: 12-15 16-24 
Nig 0-12 0-31 0-9 De 


Most of the ‘ shower ’ particles are 7-mesons, and from the results of Camerini 
et al., we find the mean value of the total energy of these shower particles is 640 Mev. 
According to Hooper (private communication), the energy released in the form of 
neutral particles other than neutrons is approximately one-half of that released in 
the form of shower particles. Hence the mean energy per shower particle may be 
taken as : 


F=0600uex. re (23) 


From Tables 9 and 10 and equation (23), we may now construct the following 
Table expressing the relation between the number of tracks and total energy 
released in nuclear disintegrations observed in C2 emulsions: 


Table 11 
Mean number of branches n 4-03 il 11-9 16-4 
Energy in emitted nucleons (Mev.) 2a 507 1183 2035 
Energy in shower particles (MevV.) 2 295 864 1980 
Total energy released FE (Mev.) a28 802 2047 4015 


The relation between the values of F and » given in Table 11 may be represented 
approximately by E(Mev.) =26n + 12-8n?. 
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ABSTRACT. The‘ prong distribution’ for the disintegrations produced by the o-mesons 
of the cosmic radiation in photographic emulsions exposed below ground is found to be 
different from that observed in emulsions exposed above ground. The underground 
distribution contains relatively more ‘stars’ consisting of one prong. The observations 
are consistent with the hypothesis that most of the o-mesons below ground are in fact 
negative fi-mesons and not negative 7-mesons. It is concluded that 8-7+1:7% of the 
negative “-mesons captured in the emulsion lead to disintegrations accompanied by the 
emission of charged particles. The average number of prongs per disintegration produced 
by the capture of negative “-mesons is 0°10+0:02. The mean number of prongs and the 
associated prong distribution suggest that a nuclear excitation of about 15 Mev. is produced 
by the capture of a w-meson. 


§1. INTRODUCTION 

HE disintegration of nuclei following the capture of slow mesons in 

photographic emulsions exposed to the cosmic radiation was first reported 

by Perkins (1947). It was shown by observations carried out at Bristol 
that most of the disintegrations of this type are produced by negative 7-mesons. 
Conversi, Pancini and Piccioni (1945, 1947) showed that the interaction between 
p-mesons and atomic nuclei is very weak, and this was confirmed by Camerini 
etal. (1948), who found that few, possibly none (9:5+9%), of the negative 
-mesons stopping in silver bromide produced observable disintegrations. 

On the other hand, Sard and collaborators (1948, 1950) have observed the 
production of neutrons by charged particles arrested in lead absorbers, and it is 
likely that several of the stopped charged particles were negative u-mesons. 
Further, Groetzinger and McClure (1948, 1949) concentrated negative mesons 
by means of a magnetic focusing arrangement into an absorber, selecting those 
which stopped in it. The long path length in the iron cores discriminated in 
favour of y-mesons, and they confirmed the production of neutrons by meson 
capture. From these two counter investigations it may be concluded that 
something like three neutrons are emitted when a negative x-meson is absorbed 
by a lead nucleus. 
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A question of some interest is whether the capture of 1-mesons may sometimes 
lead to nuclear disintegrations with the emission of charged particles. Several 
attempts to investigate this have been made at sea level or mountain altitudes, 
using either cloud chambers (Wang and Jones 1948, Chang 1949, Cool et al. 1949) 
or nuclear emulsions (Perkins 1949), but the results have not been conclusive. 
These investigations are all faced with the following difficulty: the slow mesons 
of the cosmic radiation which are arrested in absorbers or emulsions consist of a 
mixture of z- and p-mesons. The negative 7-mesons are very efficient in 
producing nuclear disintegrations containing charged particles, and even at 
sea level most of the o-meson ‘stars’ must be ascribed to 7-mesons. Hence, in 
those few cloud-chamber pictures which show heavy tracks associated with 
stopping mesons it is not clear which type of meson is involved. 

We have been able to overcome this difficulty by working underground at 
various depths up to a maximum of the equivalent of 60 metres of water. Under 
these conditions, the ratio of 7- to y-mesons is reduced by a factor which may 
be as great as 20, depending on the depth, and we have been able to isolate the 
disintegrations produced by the capture of slow p-mesons from the small 
background of stars contributed by z-meson capture. 

In order to avoid confusion we retain the phenomenological description 
“o-meson’ for those slow mesons which, when captured by a nucleus, produce 
nuclear disintegration ‘stars’, since we believe that only a small proportion of 
the o-mesons observed below ground are in fact 7-mesons. 


§ 2. EXPERIMENTAL ARRANGEMENTS 


The observations were made with Ilford Nuclear Research emulsions, 
Type G5, 200 microns in thickness. ‘They were coated, stored and developed 
underground, using stations on the underground railway in London at depths 
equivalent to 20, 34 and 60 metres of water, in a manner described previously 
(George and Evans 1950). The plates were only brought to the surface for 
examination after fixing. The possibility of confusion due to pre- or 
post-exposures of the plates was thus completely eliminated. 


§3. EXPERIMENTAL RESULTS 


In the examination of 250 cubic centimetres of emulsion, 30 o-mesons, 
5 positive 7-mesons showing 7— decay, and 1,014 «-mesons have been observed 
to come to rest in the emulsion. ‘The features of the disintegrations due to 
a-mesons will first be discussed in detail. 


(1) The Distribution in ‘ Prong’ Number of the o-Stars produced Underground 


The distribution in prong number of all the underground o-stars is given in 
Table 1. In this Table the expected distribution for an equal number of 
o-mesons of the cosmic radiation at mountain and balloon altitudes in the 
atmosphere is given for comparison (Menon et al. 1950). 


Fable 
No. of prongs 1 2 3 4 5 
Underground o-mesons 23 5 2 0 


o-mesons observed at altitude 10 8 i 4 1 
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The o-mesons at altitude are known to consist almost exclusively of negative 
a-mesons, and hence the discrepancy between the prong distributions given in 
‘Table 1 strongly suggests that the majority of the underground o-mesons are 
not 7-mesons. Since they are obviously mesons, and since ju-mesons are observed 
in large numbers, we make the tentative suggestion that most of these underground 
o-mesons are, in fact, negative u-mesons, and that they produce nuclear 
disintegrations substantially different in character from those due to negative 
m-mesons. 


(11) The Directions of Motion of the Underground Mesons 


The mesons were divided into two classes, according to whether they entered 
the emulsion from above or below the horizontal. The results are given in 
‘Table 2. 

The p-mesons show a marked collimation in the vertical direction, as also 
do the o-mesons, while such 7+-mesons as have been observed appear to be 
distributed nearly isotropically. The directions of motion of the slow 7-mesons, 
of both signs, observed under thick absorbers at the Jungfraujoch have been 


‘Fable: 2 
‘Total no. of mesons observed : 
Down Up Up (% of whole) 
[é-mesons 988 26 2-6--053 
o-mesons OF 3 10-46 
am*-mesons (showing 7—y decay) 2 3 ~60 
Table 3 
Depth below ground Frequency per cm? per day 
(metres water) o-mesons ju-mesons Stars 
20 0-0032 +0-0007 0-098 +-0-005 00-0062 +0-0012 
34 0-0011 +0-0005 0:058+0-004 0-0056-+0-001 
60 0-0005 +0-0003 0-017 +0-002 0-0047 +0-0009 


shown by Barton et al. (1951) to be distributed almost isotropically, and the 
m-mesons produced by the y-rays from the 350 Mev. electron synchrotron at 
Berkeley have been shown to be distributed isotropically also (McMillan et al. 
1949). The low energy z-mesons observed in the present investigation are 
probably produced following electromagnetic excitation of nuclei by the fast 
p-mesons of the penetrating component (George and Evans 1950), and we should 
therefore expect the underground 7-mesons to be isotropic also. It is seen from 
the Table that this is consistent with the observations for the positive 7-mesons, 
but is certainly not so for theo-mesons. The distribution of the latter resembles 
more closely that of the y-mesons, and appears to be consistent with the 
hypothesis that many of the c-mesons are pz-mesons. 


(iii) The Frequency of Occurrence of o-Mesons at different Depths below Ground 
The results are given in the second column of Table 3, and the figures for stars 
and slow u-mesons are also given for comparison. 
Under thick absorbers, the frequency of 7-mesons is known to be proportional 
to that of the stars (Harding and Perkins 1949, Barton et al. 1951). Between the 
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extreme levels of observation the frequency of stars, and hence, presumably, of 
m-mesons, only varies by a factor of at the most 1-5, whereas the frequency of 
o-mesons is seen to vary by a factor of about 6. This last figure is in agreement 
with the figure for slow u-mesons, which again supports the idea that at least the: 
greater part of the underground o-mesons are p- and not 7-mesons. 


§4. DISCUSSION 

The distribution in prong number, the angular distribution of the c-mesons, 
and the depth variation of the o-stars all point to the suggested interpretation of 
the o-mesons. ‘This is further supported by the negative/positive ratio the 
m-mesons would have if the o-mesons were 7-mesons. ‘Thus, assume that the 
o-stars are all produced by negative 7-mesons. Then from the observed number,. 
30, we may deduce that 41 negative 7-mesons came to rest in our plates, since: 
28°, do not produce a visible star (Adelman and Jones 1949). In the same 
volume, five positive 7-mesons were observed, and this would therefore correspond 
to a negative/positive ratio of 8:1. On account of the probable electromagnetic 
origin of the underground 7-mesons (George and Evans 1950) the value of the 
negative/positive ratio would be expected to be about 1:7:1, found for the 
production of 7-mesons by energetic photons. In fact, assuming 1:7 negative 
7-mesons for every positive one, and bearing in mind the proportion of ‘star-less’ 
negative 7-mesons, we are led to expect that the number of negative 7-mesons. 
giving rise to o-stars would be approximately equal to the number of positive 
m-mesons. Since we have seen five examples of 7—~ decay we may expect that. 
something like five of the 30 o-mesons are, in fact, negative 7-mesons. 

We have examined the o-stars to see whether, using the criteria laid down by 
Menon et al. (1950), we could identify stars characteristic of 7-mesons. We find 
that the two stars with 3 prongs and one of the 2-pronged stars are characteristic of 
m-meson disintegrations in light elements, and are therefore unlikely to be produced 
by »-mesons on account of the small probability of capture of the «.-mesons by the 
atoms of the gelatine. ‘Two examples of the stars ascribed to 7-mesons are shown 
in Figure 1 (Plate*). In both cases the mesons were travelling upwards, and this 
further supports their classification as z-mesons. In addition we have classified 
one other 2-pronged star as a 7-meson disintegration. From the known prong 
distribution of the 7-meson stars we should also ascribe at least one of the 1-pronged 
stars to the 7-meson group, though in this case we do not know, of course, which one. 

By the above analysis we have selected from the 30 o-mesons five cases which we 
ascribe to 7-meson disintegrations, the prong distribution being one 1-pronged,. 
two 2-pronged and two 3-pronged stars, a result in reasonable accord with the 
detailed studies made in other laboratories on negative 7-mesons (Menon et al. 
1950). ‘The distribution in the directions of motion of these particles is two down 
and three up, which again appears to be in reasonable accord with the assumption 
of a strictly local production of z-mesons. The plates were standing on the floor 
of a large platform tunnel, and assuming an original isotropic distribution in 
m-mesons, the upward direction would be favoured for those seen in the plates, on 
account of the opportunity for decay in transit for the downward travelling: 
mesons in the air-space over the plates. 

This leaves three 2-pronged stars and twenty-two 1-pronged stars which are. 
ascribed to disintegrations produced by z-mesons. The three 2-pronged stars are: 


*For Plates see end of issue. 
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‘shown in Figure 2, and are seen to be very untypical of disintegrations produced by 
-m-mesons, on account of the extremely low energy, of the order of 1 to 2 mev., of 
the emitted charged particles. 

In nineteen out of the twenty-two l-pronged stars the single heavy track 
-escaped, while in three cases it ended in the emulsion. Photo-micrographs of one 
of these three examples are shown in Figure 3. We may conclude that the mean 
range of the charged particles is greater than 300 microns, corresponding to an 
-energy of 7 Mev. for protons. 

From the measured value, 1:27, of the ratio of positive to negative particles 
with momentum 101° ev/c, at sea level determined by Owen and Wilson (1949) we 
may assume that 458 out of the total of 1,039 u-mesons were negative. Of these, 
54% will be captured in Ag and Br (Menon et al. 1950) and 20% of the remainder 
will be captured in gelatine. This gives a total of 290 negative jz-mesons which 
experienced nuclear absorption. This now leads to the distribution in the number 
of prongs for the disintegrations produced by negative u-mesons, given in the 
second row of Table 4. 


Table 4 
No. of prongs ° 0) 1 2 Sy 
No. of cases observed 265 22 3 0 
Calculated, assuming Poisson distribution 262 26 ies 0:5 


From the observations of Sard et al. (1948, 1950) and Groetzinger and McClure 
(1948, 1949) we may take it that the neutral particles emitted in the first category 
are neutrons, and further, that on the average two to three neutrons are emitted. 
The average number of prongs per disintegration is 0-10 + 0-02, and this corres- 
ponds to an average excitation of about 15 Mev. (Fujimoto and Yamaguchi 1949). 
The calculated prong distribution assuming the numbers of charged particles to 
follow Poisson’s law is given in the third row of Table 4, and is seen to be in 
reasonable agreement with observation. 

The excitation of nuclei by «-meson capture has been discussed by Tiomno 
and Wheeler (1949), asssuming a reaction of the form x” +P+N + po, where py 
is a neutral particle, possibly a neutrino. The nuclear excitation caused by the 
recoil momentum of the neutron is considered under various models. For the 
limiting cases in which the nucleus is assumed to follow the Hartree model or 
the Fermi gas model, the calculated values of the mean energy of nuclear excitation 
are 18-5 mev. and 12-6 Mev. respectively, corresponding to mean prong numbers 
of 0:14 and 0:09. Our observed values lie between these two limits. 

With such an excitation, in the case of the 2-pronged stars most of the energy is 
lost as binding energy, leaving only 1 or 2 mev. available for kinetic energy. The 
2-pronged stars would therefore be expected to have very short tracks, as is in 
fact found (Figure 2). Further, at this excitation, the nuclear potential barrier 
favours neutron emission, the theoretical number of neutrons corresponding to 
our observed prong distribution is about two, which agrees with results of the 
counter investigations discussed earlier. 


§5. CONCLUSIONS 
From the observations reported in this paper, we have drawn the following 
conclusions: 
(i) Under absorbers of thickness greater than 2,000 gm/cm? most of the o-stars 
are produced by negative u-mesons. 
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(ii) 2:5 + 0-5% of all the slow y-mesons of the cosmic radiation, or 8-7 + 1-7°% of 
all the negative .-mesons captured in the emulsion, give visible stars (cf. Camerini 
et al. 1948). 

(iii) 90°% of the stars produced by u-mesons have one prong and 10% two. 
prongs. 

(iv) The nuclear excitation corresponding to the observed prong distribution of 
the o-stars due to u-mesons is about 15 Mev., in agreement with calculations based. 
on the reaction ~~ + P+N + po. 
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Some Gamma-Rays from Light Elements under Proton 
Bombardment 


By J. H. CARVER anp D. H. WILKINSON, 
Cavendish Laboratory, Cambridge 


MS. recetved 29th September 1950 


ABSTRACT. The energy of the photo-protons from deuterium, measured in an 
ionization chamber, has been used to determine the energy of gamma-rays produced by the 
proton bombardment of several light elements. The results are expressed in the following 


Table : 


Target Proton Gamma-ray 

thickness energy energy 

Isotope . (kev.) (kev.) (Mev.) 
*Be 2 1000 e390 
uB 200 950 12:50+0:21 
13 i 8 -14--0-08 
Cc 20 570 L. 5-81-£0-25 
at thick 740 12-09 +0:-28 

AG thick 500 14-4 +0-4 

15N) unknown 900 4-45 +0-04 


The absence of certain well-known lines from the above table merely means that their 
energy was not determined in this investigation. 


section of the deuteron we have determined the energy of several gamma-rays. 

Photo-protons spent their energy in deuterium or a mixture of deuterium 
andargon. It was assumed, for the gases individually, that the ionization produced 
was proportional to the energy of the proton. Such assumptions are always 
questionable, but in this case may be well justified by three remarks. Firstly, 
the smallest proton energy involved was 1:1 Mev. and the dependence of W, the 
energy loss per ion pair, on particle energy is always weak at such particle 
velocities. Secondly, hydrogen and argon are both well known for their constancy 
of W,-even at low energies. Thirdly, Bakker* has reported values of W for 
protons of 340 mev. which differ by but 3 or 4% from those obtaining at low 
energies. 

Our measurements consist essentially in the comparison of the energies of 
the photo-protons produced by the gamma-ray in question with those of the 
photo-protons produced by the gamma-ray resulting from the proton 
bombardment of fluorine. These latter have been extensively investigated 
(see Chao, Tollestrup, Fowler and Lauritsen 1950 for references), and it is known 
that a thick target bombarded by protons of about 500 kev. results in the almost 
pure radiation of 6-14 + 0-03 mev. which has been taken as standard in this work. 
For the reduction of photo-proton to gamma-ray energy the value of 2-230 Mev. 
given by Bell and Elliott (1950) for the deuteron binding energy was used. 

The photo-proton pulses were amplified and analysed by the ninety-nine 
channel kick-sorter described by one of us (Wilkinson 1950). Owing to the 


* Report to the 1950 Harwell Conference. 


I; the course of a series of measurements of the photo-disintegration cross 
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momentum of the gamma-ray there exists a natural spread in the energy of the 
photo-protons. The magnitude of this spread is often found surprising: it has a 
minimum value of 0-196 in terms of the mean photo-proton energy. ‘This sets an 
unfortunate limit to the precision of this method of measuring gamma-ray 
energies. The method possesses, however, some notable advantages of which ~ 
the principal is the fact that its efficiency increases with decreasing gamma-ray 
energy—down to 4:5 Mey. at any rate—whereas that of the pair spectrometer 
decreases very rapidly with decreasing gamma-ray energy. The relatively constant 
efficiency of the method from 3 to 15 Mev. or so makes it well worthy of 
consideration in those studies where the relative intensities of gamma-rays, as 
well as their energies, must be known accurately. The efficiency need not be so 
low as is commonly supposed: with the deuterium-filled chamber used in this. 
work, a counting rate of about 4,000 per minute is given by 1 millicurie of 
6-Mev. gamma-rays at a distance of 10cm. The greatest disadvantage of the 
method is that it is very sensitive to fast neutrons whose flux must be less than 
0-001 of that of the gamma-rays for clean results to be obtained. 

The gamma-rays on which we report all derive from the proton bombardment 
of light isotopes and are due to radiative capture with the exception of those 
from ©N. The bombarded isotopes are : 

°*Be: <A thin target and a proton energy of 1,000kev. were used. The 
measured gamma-ray energy was 7:°39+0-15mMev. Measurements by Fowler, 
Lauritsen and Lauritsen (1948) had given 7-4+0-2 Mev.,and by Walker (1950), 
7:37 + 0-07 Mev. 

4B; Protons of 950 kev. were used with a natural boron target of thickness 
200 kev. Gamma-rays of 16-6, 11-8 and 4:3 Mev, are known from this reaction 
for a thick target and a proton energy of 900kev. (Fowler, Gaerttner and 
Lauritsen 1938). Only the middle gamma-ray was measured here and found to 
be of 12-50+0-21 Mev. ‘The most accurate measurements are those of Walker 
(1950), who finds 11-76 + 0-18 and 12-12 + 0-12 Mey. for protons of energy 510 and 
1,200 kev. and thick targets. 

13C: Protons of 570 kev. were used with a target of the separated isotope 
of stopping power about 20kev. Lauritsen, Lauritsen and Fowler (1941) have 
reported 8-1 and 5:8Meyv. gamma-rays of equal intensity from this resonance. 
We found a strong gamma-ray of 8-:14+0-08Mey., but only a weak radiation of 
energy 5:81 + 0-25 Mev. whose intensity was estimated as about 0-07 of that of the 
stronger line. It is possible, but not, we believe, likely, that this weaker line is 
due to fluorine contamination (6-1 Mev.). The measurements were carried out 
at 0° to the proton beam. 

19F: The radiative transition from the 669kev. resonance discovered by 
Devons and Hereward (1948) was measured and found to be of 12-09 + 0-28 ev. 
and to possess an intensity 1-72 +0-25% of that of the combined 6- and 7- Mev. 
lines from this resonance. ‘This energy is in accord with the 12:0+0-2 Mev. 
reported by Rae, Rutherglen and Smith (1950) and confirms their remark that 
the transition takes place not to the ground state of ?°Ne but rather to the excited 
state at about 1:5mev. known from other measurements. These two 
determinations of the gamma-ray energy locate this state at 1:50+0-2mev. _ 

7Li: Protons of 500kev. bombarded a thick target of metallic lithium. 
Gamma-rays of 14:8+0:3 and 17:-6+0-2Mmev. at the 440-kev. resonance are 
known from the work of Walker and McDaniel (1948). We measured the softer 
gamma-ray only, finding an energy of 14-4 + 0-4-Mey. 
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UN; The gamma-ray in this case comes from.!2C which results from 
alpha-particle emission following proton capture at about 900kev. An excited 
state of *C, formed in the °Be(«n)!2C reaction, at 4:40+0-05 Mev. has been 
reported by Pringle, Roulston and Standil (1950), and at 4-45 mev. by Bradford 
and Bennett (1950), Walker (1950) finds 4:43+0-15 from the reaction 
UB(py)’C. 'T. Lauritsen * has reported gamma-rays of energy 4-47 Mev. from the: 
reaction here studied; we find an energy of 4-45+0-04Mev. The gamma-rays. 
are markedly anisotropic at 900kev., having a distribution of the order 
1+0-3 cos? @, though interference terms are probably present and the distribution 
depends somewhat on the proton energy. A target of ®N on molybdenum 
was used. 
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CLEVIERS: TO THE- EDITOR 


Fall and Re-increase of the C(y,3*He) Cross Section 


Some preliminary results of an investigation of the photo-disintegration “’C(y, 3 *He), 
using nuclear emulsions irradiated by gamma-rays from the A.E.R.E. synchrotron, have 
already been reported (Goward, Telegdi and Wilkins 1950). This investigation has been 
continued, with peak gamma-ray energies of up to 25-5 Mev., and analysis of about 
600 stars confirms and extends some interesting trends shown by the earlier results. 

Figure 1 shows a histogram of the total energy release Ey of all the stars observed whose 
measurements indicated satisfactory momentum balance as carbon disintegrations 
(|A| <0°8 in the notation of Goward and Wilkins (1950)). The superposed cross section 
curve was calculated by assuming a theoretical gamma-ray spectrum (Rossi and Greisen 
1941), and is corrected for the appreciable influence of random errors in Ey. These 
errors are not closely known, but a reasonable assumption is a gaussian curve of total 
width 1-2 Mev. at half amplitude. The cross-section curve then shows a considerable 
and effectively vertical drop at Ep ~12 mev., or E,, (gamma-ray energy) =19-3 Mev. The 
subsequent re-increase in cross section is now firmly established by the observation of some 
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‘80 stars with Ey>14 Mev. The absolute cross-section values are considerably higher 
than the earlier estimate, mainly because of more experienced and thorough emulsion 
searching and more adequate correction for escape factors and overlooked stars. A 
reasonable estimate of the maximum cross section is 2:4 x 10-?8 but this may easily still 
be an underestimate by about 50%. 

In considering possible explanations of the cross-section curve shape, it must be 
remembered that most of the Figure 1 histogram was obtained by using the range-energy 
‘data of Lattes, Fowler and Cuer (1947); the more recent data of Rotblat, corrected to 
normal atmospheric conditions, would give Ep (and E,) values some 0:5 Mev. lower. ‘Thus 
the main cross section fall may well occur near E,=18-7 Mev., at which energy McElhinney 
et al. (1949) found evidence for a rapid rise in the 2C(y, n\UC cross section from zero to 
a value of the order of 10-2? cm?. 
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The re-increase of cross section at higher gamma-ray energies might suggest some 
alternative mode of disintegration to that already established for stars with Eq< 14 mev., 
where the predominant mode of disintegration was via the 3 Mev. excited level of 8Be. 
Analysis does not give significant support to this suggestion, however, as is shown by 
comparison of the experimental E* histogram and calculated curve in Figure 2 (a). 

A significantly increased proportion of stars is now found to involve ground state ’Be 
nuclei; the true overall proportion with correction for searching losses appears to be in 
excess of 10°% in the present experiments. Such stars are indicated by shading in Figure 1, 
the proportion apparently increasing with energy. ‘The considerable certainty with which 
these stars are identifiable is shown by Figure 2 (6). The fitted curve indicates the 
histogram shape expected if 43 stars involve ground state *Be, with a disintegration energy 
of 0:1 Mev. and negligible real level-width. 


Atomic Energy Research Establishment, J. J. WILKINS. 
Harwell, Didcot, Berks. F. K. Gowarp. 
17th October 1950. 
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The Reaction **Na(d, «)”1Ne 


t 

Some recent measurements on the reaction 2°Ne(d, p)?4Ne have pointed to the existence 
of levels in ?'Ne with excitation energies of about 0-35, 1-8 and 2-8 mev. (Elder, Motz and 
Davison 1947, Zucker and Watson 1950, Ambrosen and Bisgaard 1950). The reaction 
*8Na(d, «)?4Ne might be expected to lead to the same states, but earlier measurements 
by Murrell and Smith (1939) had revealed only the transitions to the ground state and to 
the 1:8 Mey. level of 74Ne. A new study of this process therefore seemed to be of interest. 

The energy distribution of alpha-particles from the reaction 2*Na(d, «) was studied 
with the aid of a differential ionization chamber and an air absorption cell. This technique 
has been described in an earlier letter (French and Treacy 1950). Thin targets of NasSO, 
or NaCl were used, with deuteron bombarding energies between 0°83 and 0-93 mMev.. 
In addition to the alpha-particles from the process *C(d,«) in a carbon contamination. 
layer on the target surface, four groups of alpha-particles from the disintegration of sodium 
were detected. With the Na,SO, targets, protons from the reaction !O(d, p)!7O* were 
also observed. The Figure contains the results of typical measurements made at 0-83 Mev. 
deuteron energy with a target of NaCl, about 0-4 mm. air equivalent evaporated on brass. 
In curve (a) only three of the sodium groups can be seen, the fourth being masked by the- 
strong peak from carbon. The chief object of this particular run was to demonstrate: 
the doubletstructure at the high energy endof the spectrum; in other runs, with fresh targets: 
of about twice the thickness, the fourth group was clearly resolved from the carbon peak 
(see curve (0) inset). 


Counting Rate (arbitrary units) 


0 10 20 30 40 50 60 70 
Absorber Pressure (cm.Hg) 
Alpha-particle groups (A, B, C, D) from deuteron bombardment of sodium at 0°83 Mev. Group D 


masked by the carbon group E in the main figure, appears clearly (inset) when a fresh target 
is used. 


Five separate determinations were made of the energy of each of the four groups from 
sodium. For calibration of the range scale, the following were taken as standards: (a) alpha- 
particles from ThC, mean range 4:70 cm.; (6) alpha-particles from *C(d,a), a Q value of 
5-20-+0-05 Mev. was assumed, after consideration of all the available data (Hornyak et al. 
1950); (c) alpha-particles from 'O(p,«). Q=3:97-40-05 Mev. (Freeman 1950); 
(d) protons from 1*O(d, p)7O*. Q=1-049 Mev. (Buechner et al. 1949). 

The new range-energy curves of Bethe (1950) were used in the analysis, although we 
have been forced to conclude from the accompanying text that his curve for protons refers 
to standard air, not to air at N.T.P. as stated on the graph. Tt 

In the Table below are given the Q values with their probable errors for the different 
groups and the corresponding excitation energies in *4Ne. ‘The values suggested by: 
Hornyak et al. (1950) on the basis of the ?°Ne(d, p) experiments are given for comparison. 

Our value of 0:38-0:02 mev. for the excitation of the first level in 4Ne agrees well 
with the figure. of 0-39 Mev. found by Ambrosen and Bisgaard (1950). The. smaller 
splitting (0-31 mev.) found by Elder, Motz and Davison (1947) and by Zucker and Watson. 


+ Professor Bethe has since confirmed our conclusion in a private communication. 
14-2 
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This work Hornyak et al. 
Group Excitation Q Excitation 
(Mev.) (Mev.) (Mev.) (Mev.) 
A 6:84-0:05 0 6°86 0 
B 6-46+0-04 0:38 -+0-02 6°52 0:34 
(© 5-09-+0-04 1:75-0:04 5-06 1-80 
D 4-08 +0-06 2:76+0:05 4-02 2-84 


(1950) may be due to the less complete resolution of the two groups in their experiments. 
For the higher levels our results show satisfactory agreement with the values chosen by 
Hornyak et al. (1950). 

Our figure of 6-84-+-0:05 Mev. for the energy release of the ground state transition in 
23Na(d, «)*4Ne agrees well with that computed from the values 2:35+0-04 Mev. in the 
process ?°Na(p, «)?°Ne (Freeman 1950) and 4:51-0:07 Mev. in ?°Ne(d, p)*'Ne (Hornyak 
et al. 1950). 

We are grateful to Dr. W. E. Burcham for pointing out that a new study of the reaction 
23Na(d, «) was desirable. One of us (D. M. T.) is indebted to the Department of Scientific 
and Industrial Research for a maintenance grant. 


Cavendish Laboratory, A. P. FRENCH. 


Cambridge. - D. M. THomson. 
28th October 1950. 
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Nuclear Spin of Erbium-167 


The paramagnetic resonance spectrum of a single crystal of erbium ethyl sulphate, 
diluted with lanthanum ethyl sulphate, has been observed at 1:2 cm. wavelength. For all 
directions of the external magnetic field a single strong line is found; centred on this is a 
set of eight almost equally spaced weak lines. The relative intensity of these lines, and 
their variation with angle, is consistent with their assignment to the hyperfine structure 
of the single odd isotope of mass 167 (natural abundance=22-8% (Leland 1950)). It 
follows that the nuclear spin of this isotope is 7/2. ‘This is the same as that of neodymium 143 
and 145 (Bleaney and Scovil 1950), and these three isotopes are the first odd-neutron nuclei, 
to have been discovered with spin 7/2. 
ie Lhe ee has been analysed in terms of the usual Hamiltonian, 


= B{g,H,S,+8\(H,S,+ H,Sy)}+ 48,14 BS 2. at Syly)+ OU — 41(I-+1)}, 


and the sees. values are found for the parameters (the effective value of S is taken as 3, 
since only one electronic transition is observed) : 


g| (parallel to axis)=1:47-++0-03; g, (perpendicular to axis) =8-85+0-2; 
A=0-0059 cm~*+-0:0002; B=0:0307 cm-!+0-0005; Q<0-:005 cm“ 


The anisotropy in the g-values is by. far the largest yet discovered. Susceptibility 
measurements of van den Handel (1940) reveal very little anisotropy at temperatures down 
to 15° K.; this result can only be reconciled with ours if there are other very low-lying 
levels besides those giving rise to the transition observed by us. This is also suggested’ 
by the very short spin-lattice relaxation time (see below). 
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The ratio A:B is nearly the same as the ratio &| :&1, as was also the case with 
neodymium. As a result, the width of the hyperfine structure in gauss as observed at 
constant frequency is almost the same in all directions. The theoretical work of Elliott and 
Stevens (1951) is in accordance with this, and shows also that no effect due to a quadrupole 
moment should be detectable. Our upper limit for Q is based on the absence of observable 
deviations from the normal spacing of the hyperfine structure, which would be of the order 
of Q?/B. 

The line-width showed an interesting effect, being dependent on the direction of the 
external magnetic field, and varying inversely as the g-value. This behaviour is to be 
expected if the dominant broadening is due to spin-lattice relaxation, and the relaxation 
time 7 is independent of the orientation of the magnetic field. We have 

eS, i h 1:14 10-7 

*“— Yacht IngBAH gAH -’ 

where Av, AH refer to the half-width at half intensity in wave numbers and gauss 
respectively. At >=0-845 cm@! the values of AH parallel and perpendicular to the axis 
respectively were 81+5 and 14+1 gauss, giving r=0-:95+0:05x10-® second at the 
lowest temperature, estimated to be 131°. from the vapour pressure of the liquid 
hydrogen. At 20° k. the lines were about ten times broader, showing that 7 varies roughly 
as the inverse seventh power of the temperature. It should be noted that as the orientation 
of the crystal is changed, the magnetic field changes, but the frequency of the transition 
remains constant. "The measurements show that 7 also remains constant under these 
conditions. ‘The theory of Van Vleck (1940) for caesium titanium alum shows that in strong 
fields the relaxation time should be isotropic and independent of the external field, if the 
“Raman processes’ predominate in the relaxation mechanism, as is certainly the case in 
our experiment. ; 

No absorption lines could be obtained in the ethyl sulphates of cerium, samarium, 
dysprosium and ytterbium at 1:2cm. wavelength, 20°K. and fields up to 12 ke. 
Measurements are now in* progress on the bromates of these elements. Work is also 
planned to observe the hyperfine lines at zero magnetic field, by which means a value of 
QO of 1/100 of that given above should be detectable. 

Our thanks are due to Mr. J. K. Marsh for the generous gift of the samples of erbium 
and other rare earths. 


Clarendon Laboratory, B. BLEANEY. 
Oxford. H. E; D. Scovir. 
14th November 1950. 
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A Preliminary Survey of the Paramagnetic Resonance 
Phenomena observed in Rare Earth Ethyl Sulphates 


In recent letters Bleaney and Scovil (1950, 1951) have reported the results of their 
experimental work on the rare earth ethyl sulphates of the form M(C,H;SO,)3;.9H,O. 
Here we report the results of some preliminary theoretical investigations using a crystalline 
electric field hypothesis to account for the magnetic properties at low temperatures. In 
most of our work so far we have used only first-order perturbation theory, with the result 
that we have only sought approximate agreement with experiment. 

From the crystal structure of cerium ethyl sulphate (Ketelaar 1937) it seems reasonable 
to assume that the ion is embedded in a crystal field possessing D3, symmetry. Using 
potential terms up to the sixth degree and neglecting terms of odd parity (which give no 
contribution in first order) our potential takes the form 


V=A,+A,(32?—r?) + A,(3524—30r?2?-+ 374) + Ag(2312°—31 5r224-+ 105r422— 578) 
+ A 8(x®—15x4ty2+ 15x?y4—y*). 
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The coefficients A can be estimated for cerium, and we take slightly different values for 
the other ions. 

Assuming that /< spin-orbit coupling, we can treat it as a perturbation on the lowest J 
manifold, and our results for the ions with an odd number of electrons are as follows : 

Cet+++ 4f! 2F,).. The sixth-order terms have no matrix elements inside J =5/2 so the 
‘states are eigenstates of J,. We find the doublet spanned by J,=-+5/2 lies lowest, which 
accounts for the absence of a resonance. 

Ndt++ 4f% 4Ig.. Neglecting A,®, we find Jz=+7/2 is lowest with Jz=+5/2 near. 
‘The effect of A,® will be to mix these states and change the lowest states to 


ali) +b|—%) and al—{)+6|2). 


Taking |a|=0-9 and |b|=0-4 gives g,=3-54 and g,=2-0, whereas experiment gives 
gy= 3°61 and g,=2-05. The only other states for which transitions are allowed are 
J,=+1/2, which yield g,=8/11 and g,—40/11. 

Smt+++ 4f5 °H.).. The sixth-order terms have no matrix elements and V gives 
J,=+1/2 lowest. The corresponding g-values are g,=2/7 and gt=6/7. No transitions 
have been observed as these values lie outside the range of observation. 

Gdtt+ 4f7 8S, Since this ion is in an orbital S-state, to our approximation the 
crystalline field gives no splitting, and other effects become important (Pryce 1950). 
g is isotropic and closely equal to 2. 

Dy*t** 4° ®H,;;5. Our potential field gives the ground states as approximately 
J,=+15/2, to which are admixed J,=+3/2, and there will be no transitions. If our 
estimates of the crystalline fields are slightly inaccurate (A4,,7” varies from ion to ion) it 
may happen that J,=-+13/2 is lowest, and to this are admixed small amounts of J,=+1/2 
and J,=+11/2. This would lead to small g-values which would not have been observed 
and which would appear to -be in disagreement with susceptibility measurements of 
-van den Handel (1940). 

Ert++ 4f11 41,5... The crystalline field without A,® gives J,=+5/2 and J,=+7/2 
close together at the bottom. When these are admixed we get the states a\i) +b|—8 
and a|—%)+6|8) with |a|=|b|=1/2. The corresponding g-values are g,;=1:2 and 
g,=8-9. On the other hand small changes in the crystalline field may give J,=+1/2 
lowest, to which are admixed small amounts of J,=-+11/2 and +13/2. The corresponding 
g’s are approximately g,y=1:2 and g,;=—9-6, and with the experimental g’s of gy=1-47 
and g,;=8-85 we cannot, at present, decide which levels are actually lowest. 

Ybt*+* 4f18 2F./.. Here we find J,=-+3/2 is the lowest level, and there will be no 
transitions in agreement with experiment. 

We now turn to the consideration of the nuclear hyperfine structure. The absence of 
hyperfine structure in gadolinium points to the fact that the admixture of states containing 
unpaired s-electrons (Abragam 1950) is small, and this will probably be true for the other 
ions also. With our approximate electronic states the hyperfine structure arising from 
magnetic interactions should show the same anisotropy as the g’s, as is approximately 
correct. The deviations are probably due to admixtures of states in other J manifolds 
{particularly J’=J+1) to the ground states, and show that second-order effects are not 
negligible. Using the values given by Bleaney and Scovil for the constants in their 
Hamiltonian, we can estimate the magnitudes of the nuclear magnetic moments and set 
upper limits to the nuclear electric quadrupole moments. These values are liable to 


considerable error as it is necessary to estimate 1/r? (Sternheimer 1950). The 


accompanying ‘T'able gives our results. In particular the upper limit set to the quadrupole 
moments is found to be large. 


A a) wd/rP) = Q/*) Ar’) “ Q’ 
143d 0-040 <0-002 47 <870 34 1°4 —<<o)iy 
145d 0-025 <0-002 29 <870 34 0°85 DD 
167 5 C << 740 ° <10 

Er 0:0059 <0:005 45 ee 75 0-6 . <17 


This work is continuing and will be reported more fully in due course. 
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A and Q are as defined by Bleaney and Scovil and are in cm=! ; 1//* is given in units 
of 10% cm~*; pw is given in nuclear magnetons; Q’ is the quadrupole moment in units 
of ex 10-*4 cm** and is defined (Mack 1950) by Q’=Jp,(3z*—r*)dr, where p, denotes the 
charge density for my=I. The two upper limits of Q for erbium arise from the uncertainty 
in the ground state. The smaller value is due to |+4). 


Clarendon Laboratory, R. J. Evviort: 
Oxford. K. W. H. STEVENs. 


14th November 1950. 
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Scattering of Neutrons by Ferromagnetic Crystals 


Consider neutrons with wavelength greater than the Bragg reflection wavelength scattered 
by a ferromagnetic crystal. Then inelastic magnetic scattering in which the neutrons lose 
energy by disorienting atomic magnetic moments within a domain would be expected to be 
small from energetic reasons. J. M. Cassels has suggested to the author that at temperatures 
not too far below the Curie point where there are an appreciable number of disoriented 
atomic moments within a domain, some of the disoriented electron spins may be reoriented 
during scattering into the direction of magnetization. This would result in loss of exchange 
energy by the scatterer, and consequent gain in energy by the neutrons. 

The purpose of this note is to report the result of an investigation of this effect using 
Bloch spin wave functions to provide a simple model of a ferromagnetic, with just one 
magnetic electron per atom. The approximate theory of Bloch spin waves (Seitz 1940) 
in which the waves are assumed to be non-interacting, is probably good up to room tempera- 
tures in iron. Using this theory the total cross section for inelastic magnetic scattering of 
monochromatic neutrons, with wavelength greater than the Bragg reflection wavelength, by a 
single crystal can be written c=Yo(n). (27/a)n is a vector in the reciprocal lattice, a being 


n 
the constant for a cubic lattice. We take the lattice constant and the Curie temperature and 
the body-centred cubic structure for iron, but just one magnetic electron per atom. 

In the scattering process discussed here, the neutron receives two contributions to its 
momentum. One is from the spin wave which in practice is always small, due to the energy 
relation. ‘he other is from the periodic field and corresponds to the momentum, though 
not, of course, the energy, relation appropriate to a Bragg scattering. Each o(n) is then the 
cross section corresponding to such a momentum relation. The first term, with n=0, is 
small angle scattering, where the energy transfer is small. For n not equal to 0, o(n) 
represents a scattering into a small solid angle of the order 4 steradian round the direction of 
the vector p—(2z/a)n where p is the wave number of the incoming neutrons. (In fact 
the momentum equation is p’=p—(2m/a)n+k, where p’ is the wave number of the 
scattered neutron, and k is the wave number of the annihilated spin wave and is very small.) 
Since |p|<4(2z/a)|n|, scattering occurs in directions near the directions of the reciprocal 
lattice vectors. It is found that only the terms for n?=2 are important, due largely to the 
decrease of the atomic form factor for large |n|. 

As an example, for a crystal temperature of 400° x. for n=(1,1,0), p=(27/a)(2,0,0) and 
for an unmagnetized crystal and unpolarized neutrons, it is found that o(n) is about 4. This 
corresponds to a differential cross section of 4 barn per steradian. 

If this effect were to be observed it would provide a test of the Bloch model. The 
possibility of observing it depends on the large differential cross section for wavelengths 
at which the nuclear scattering is isotropic, and of the order of magnitude of about 8 barns. 
Moreover, the energy of the scattered neutrons will be greater than that of the i oming 
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neutrons, the final wave number corresponding to o(m) being given approximately by 
p’|=|p—(2z/a)n|. As the direction of scattering corresponding to a(n) does not depend on 
pl, for p in the direction of n, such neutrons need not be monochromatic, provided the 
wavelength is greater than the Bragg reflection wavelength. 

In real iron, instead of just one spin per atom there are approximately two. The likely 
effect of this would be to broaden somewhat the angular range of the scattering corresponding 
to each o(m).. Such a change would probably not, however, invalidate the qualitative 
conclusions of this work. 


Department of Mathematical Physics, R. G. Mooruovuse. 
University of Birmingham. 
4th November 1950. 


SriTz, F., 1940, The Modern Theory of Solids (London and New York: McGraw-Hill), pp. 617 ff. 


A Note on the Radially Symmetrical Phase Growth 
Controlled by Heat Conduction 
In a recent paper Frank (1950) gave a function for the space-time dependence of the 


temperature in radially symmetrical heat conduction problems, when a sphere of new phase 
is growing. His function is 


U(s)=U ~+ Al[s—exp(— $s) — $2792 {1 —erf(4s)}]. =... (1) 


The temperature is dependent upon the only variable s, which is a combination of the 
space- and time-coordinates r and ¢: 


S=SfaTtie eo Vy ee eee (2) 


The temperature of the boundary is always Us(S)=const. (equilibrium temperature between 
the two phases) and its velocity of propagation is (S marks the boundary) 


vdrdt=}Sat? = ies =e Ge eee (3) 


There is another means of describing the temperature around a growing sphere of new 
phase. ‘The source function for a point heat source is, in 7-dimensional space 


U(r, P= QO a cyatit*) texpl-4eb. 6 ae (4) 


‘The heat is generated on the surface of the growing sphere and the corresponding temperature 
distribution is given by an integral over the surface having the value 


_ OQ 1S pe ipr \1-t” . 2pr 
U(r, j= War oy ayn *P (- 4a°t 4a°t } lr (4n— 1)J yn—s (2 an 40°) DR EOC GD) 


where Q is the quantity of heat, c the specific heat, y the density, a? the thermal conductivity 
and J the Bessel function of the first kind. One has also to integrate over all the time that 
heat is generated and this gives (for n=3) 


Qop(7)v(7)dr [r—pi7)]? [r+ p(z)]}* 

t es — =>] = Aronc a 
eas If 7=0 27 ®eyar(t—7)* ||P Gare—ay JP “4aXQt—a) Jf? (©) 
Qy is the quantity of heat for unit volume. p= Jv(n)dn and v are the position and velocity of 
boundary betwecn the two phases. 


In the above-mentioned case p and vw are determined by equations (2) and (3). Then 
equation (6) becomes 


Ow (2 “ZA—1 
Scent pe 0 (Z?+-1)? 


where Z?=(t'/?+ 71/?)/(t/?— 71/2), Integration with known formulae (see, for example, 
Magnus and Oberhettinger 1948) gives the function of equation (1) as it must be. 


> exp (— fs? Z?)dZ suckanatered (Ay) 
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: The sediconlions of the integral (6) lies in the possibility of choosing any function for v 
’ (and p) which ay be desirable = any problem suggested by theoretical consideration or by 


F2 A Sith haecessineedh cese of nee growth, in which the velocity v is assumed proportional 
” to the difference between the equilibrium temperature Ug and the temperature of the 
_ boundary U, s=C (Us—U), was treated by Geist und Dehlinger (1949). The temperature 
of the boundary U lies beneath Us on account of supercooling. To use the one-dimensional 
analogue to equation (6), (6’), in this case, obviously one must first determine v as a function 
_ of t from an integral equation resulting from the velocity law inserted into equation (6’). 


Nimberg, : D. GEIST. 


| Franx, F. C., 1950, Proc. Roy. Soc. oe 201, 586. 

_ Geist, D _ and DEHLINGER, U., 1949, Z _ Naturforsch., 4a, 415. 

é atscxcs, W., and OBERHETTINGER, F., 1948, Formein und Sétze fiir die speziellen Funktionen der 
mathematischen Physik (Berlin : "Springer-Verlag. ) p. 126. 


Calculation of Relative Transition Probabilities for First 
Negative Bands of Nz 


The relative transition probabilities between the lower vibrational levels of the two 
electronic states corresponding to the transition x*2; — BX; in Nj have been calculated here 
| _ using the simple harmonic wave functions. ‘Wace edn wave functions are well 
| known and take the form 


; S,=d7?(2’v ! o/z)-17H,(x/d) exp {—3(x/d)}, pe Seeet) 
where x=r— r, (7, being the equilibrium internuclear distance) and wv is the vibrational. 
quantum number. Further d?= =h/47 cue, where p is the reduced mass of the system and 


@, is the usual quantity occurring in the expression for the vibrational energy. The 
relative transition probabilities depend on the expression 


P(n', v3 n’ =| [- S(n’, v' |r) S(n”, v” |r) ae 


‘The p values for various (z’, v«”) numbers were calculated using the wave functions given in_ 

_ equation (1), the integration over the internuclear distance 7 being carried out graphically. 

: The values of the constants for the two electronic states, tegether with the reduced mass- 
of the system, are given in Table 1 (see Herzberg 1950). 


es Table 1 
a p=1-162 x 10-** gm 
p r, (10-8 cm.) w, (em). d(10-2° cm.) 
4 xEz 1-116 2207 4-67 
. Es 1-075 2420 4-46 
: Table 2 gives the final p values, together with the values (in brackets) obtained by Bates- 
. (1949) using the Morse potential function. The experimental values obtained by 
a Herzberg (1928) are given in square brackets. 
: Table 2 
ee \o" 0 1 2 
a oS 0-69 (0-65) [0-69] 0-34 (0-26) [0-26] — 
1 0-29 (0-30) [0-23] 0-26 (0-22) [0-27] = 
2 0-02 (0-05) [0-07] = = 
3 0-00 (0-00) [0-00] —_ 0-29 (0-41) [0-36 
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It appears that good agreement with experiment is obtained by using simple harmonic 
oscillator wave functions and the use of complicated Morse wave functions hardly seems 
justified in the view of the accuracy. For higher (v’, v”) numbers, however, where increased 
cancellation of the wave functions occurs the simple harmonic wave functions will certainly 
be useless and, although the Morse functions should give fairly accurate results, the calcu- 
lations become extremely lengthy. 

It would be of interest to apply the method of distorted wave functions (as adopted by 
Pillow 1949, 1950) to the transitions considered above. 


Department of Physics, G. STEPHENSON. 
Imperial College, London, S.W.7. 
23rd November 1950. 
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C, (Swan) Bands in Krypton 


Swan bands are known to be excited, among other methods, by inert gases. Among the 
gases that have so far been used are helium, neon and argon. While studying the selective 
efficiency of the production of these bands at various pressures of the above three gases in 
a discharge tube having carbon electrodes, some connection was noticed between the 
efficiency of the bands and the excitation potentials of the gases. ‘Thus the bands were 
much more efficient and purer in argon than in neon, and better in neon than in helium. 
The order of efficiency was observed to be generally in inverse proportion to the excitation 
potentials of the gases; krypton, having less excitation potential than argon, was expected 
to act as a still better exciter of Swan bands. This suggestion has led us to investigate these 
sbands in krypton. 

The discharge tube was of the H-type fitted with carbon electrodes. It was filled with 
mercury-free krypton gas, at successive pressures in the range 5 mm. to 5C mm. Hg, and the 
-discharge was excited by a neon sign transformer, with its primary fed from 230 Vv. A.c. 
mains through a variac. The best conditions of gas pressure and current density were 
-established to obtain bands as intense and pure as possible. As expected, the bands came 
out with comparative ease, as in argon, and probably a little better. The relative intensities, 
from visual estimates, for the different composite sequences of the bands for equivalent 
pressure and current density conditions are given below together with the impurity spectra. 
“The figures are based upon 100 for the (0,0) band in every case. 


i ‘Ue Helium Neon Argon Krypton 
st excitation : 
potential 19-77 16°58 11-57 9-98 
C, (Swan) sequences 
Av=+2 40 30? 120 50 
sl 200 300 380 310 
0 160 200 210 230 
—1 330 300 1 350 430 
“3 overlapped 
2 320 by Naive, 360 470 
Comet tail (CO*) prominent pene absent absent 
Triplet (CO) weak absent absent absent 
CH (4315) strong weak prominent weak 
H and He Ho gcebt Ali a few 
Atomic lines lines fair N ee £ ble =a - oy prominent 
in intensity si eone irene ah Kr lines 


It appears from these observations that, among the inert gases investigated, krypton acts 
as a more efficient source than helium or neon for excitation of Swan bands, and is more 
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-or less comparable with argon, if not a little better. Xenon (V,.=8-39), which is expected 
to be still better, is under investigation in this laboratory. In this connection it may be 
‘noted that there is evidence of selective excitation of the bands at a definite region of 
pressure characteristic of each gas. It appears that this excitation is a consequence of the 
-energy of electrons in the discharge (depending on the excitation and ionization potentials 
_of the gases), metastable states and the probability of the collisions of the second kind; 
‘this is discussed in a paper to be published in Current Science. A more detailed paper 
giving quantitative estimates with relevant aspects is under preparation. 
The authors are indebted to Dr. N. R. Tawde for his helpful guidance and discussion. 


Spectroscopic Laboratories, G. K. MEHTA. 
The Institute of Science, V. RAJESWARI. 
Bombay, India. 
24th November 1950. 


The ’C(n, 2n)"'C Reaction in an Anthracene Crystal 


When the 110-inch cyclotron at the Atomic Energy Research Establishment, Harwell, 
began operating recently, it was necessary to check that in those parts of the building which 
were occupied for any length of time the neutron flux (of maximum energy 175 Mev.) was 
less than the specified biological tolerance Jevel. 
For neutrons of energy greater than 20-4 Mev. this has been carried out satisfactorily 
by the use of the #2C(n, 2n)C reaction, the threshold for which occurs at the above 
neutron energy (Sherr 1945). The “C decays with a half-life of 20-35--0-08 minutes, 
-emitting a positron of maximum energy 0:97 Mey. (Mattauch and Flammersfeld 1949). 
The cross section for the reaction has been calculated theoretically by Heckrotte and 
Wolff (1948). It is, to a rough approximation, constant, and the calculated value at 90 Mev. 
(0-012 barn) is in substantial agreement with the value obtained experimentally by McMillan 
-and York (1948) at 90 Mev. (0:022+0-004 barn). 
The maximum permissible flux of fast neutrons of energy greater than 5 Mev. has been 
-calculated by Nearyand by Tait (A.E.R.E. private communication) as 30 neutrons cm~? sec™?. 
It is desirable to be able to measure a fraction of a tolerance dose, and so a sensitive method 
«of detection is necessary. 


15 2:0 25 
Total Energy, € (72,¢’ units) 


High sensitivity has been achieved by irradiating an anthracene crystal (mass 4:48 gm.) 
-and counting, with an E.M.I. photomultiplier tube, the scintillations produced in the crystal 
itself by the positrons from the decaying “C nuclei. It is possible to detect a neutron flux 
of 15 neutrons cm~? sec! by this method. It involves an irradiation time of 30 minutes 
andasimilar counting time. The method is not affected by the presence of gamma radiation. 
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The technique also provides an interesting example of the use of a scintillation counter 
for the study of a beta-ray spectrum. After the crystal had been irradiated for 30-5 minutes 
in a flux of 2X10® neutrons cm™ sec, a- pulse-size distribution was obtained using a 
single-channel ‘ kick-sorter’. The ‘ Fermi diagram’ derived from this distribution is 
illustrated in the Figure. The results compare favourably with those obtained with magnetic 
beta-ray spectrographs (Siegbahn and Bohr 1944). 

The scintillation counting apparatus and crystal were made available by Mr. R. B. Owen. 
The authors wish also to thank the Director, Atomic Energy Research Establishment, and 
the President, South African Council of Scientific and Industrial Research, for permission. 
to publish this note. 


Atomic Energy Research Establishment, J. SHARPE. 
Harwell, Didcot, Berks. G. H. STAFFORD. 
27th November 1950. (of South African C.S.I.R.). 
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The Reaction °Li(yn)*Li and Energy Levels of the *Li Nucleus 


In the course of photo-disintegration experiments using lithium loaded emulsions 
(Titterton 1950 a, b) evidence was sought for the reaction 


SLit+hy=n+5Li—Q,; *Li='H+*He+Q, © ...... (1) 
and while some events were observed which were compatible with this disintegration scheme, 
they were too few to be sure of the assignment. 

Recently Ilford E, emulsions loaded with the separated isotope *Li have been irradiated 
with gamma-rays from the 440 kev. *Lipy resonance and have yielded definite evidence 
for the reaction. The recoiling >Li nucleus (having a lifetime of about 10-*'second) breaks 
up into a proton and an alpha-particle to form a V-track in the emulsion. In some cases. 
the discrimination is good enough to identify the proton, in others itis not. Three typical 


events are shown in Figure 1 (a) and (4), where two correspond to break-up from the ground 
state of *Li and the third to break-up from an excited state. 


Figure 1 (a). V-tracks resulting from Figure 1 (6). V-track resulting from an 
break-up of *Li nuclei in ground state. excited state of *Li. 


The energy release Q on the break-up of a compound nucleus into fragments of mass. 
m, and mz moving in paths diverging by @° in the laboratory frame is given by 


1 
Q= POTS {m E.+ mE, —2(mym.E,E,)'? cos 8}, wa Geta 


ey 


— 
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where E, is the energy of the particle of mass m, and E, that of mass ms. For the ®Li case 


this reduces to | 
O=KE, 1 4E,— ELE, puizens 0}. alt lio Vs Boke (3) 


A second relation can be derived from the conservation conditions and, for the 17-6 Mev. 
y-ray, is 
(E,,+ Ep) =5Q/6+-2°3, RA (4) 


where all energies are measured in Mev., the y-ray momentum being neglected. 

Each event resulting from this y-ray must then satisfy (3) and (4) simultaneously. 
Some events are due to y-rays from the broad 14-8 Mey. group. However, as the y-rays 
of this group have only half the intensity of those due to the sharp line at 17-6 Mev. (Walker 
and McDaniel 1948) and the energy resolution of the experiment is not better than 300 kv., 
it is reasonable to apply equation (4) to all events. 

Seventy events have been observed, and in about 50% of the cases it is clear by inspection 
which track of the V is the proton. In these cases the appropriate values of (E,+E)) and 
O are calculated directly. For the remaining events each track is assumed to be first a _ 
proton and then an a-particle and the two possible values of Q and (E,+E,) czlculated. 
The assignment giving values simultaneously satisfying equations (3) and (4) is accepted. 
The histograms of Figure 2 show (a) number of events plotted against measured energy 
telease (E,+E,), (b) number of events plotted against Q. 

In each case definite groups a and a’ are found corresponding to break-up of *Li from 
the ground state. A value for the instability of the state is found to be O2>1 6+0-2 mMev., 

_ and from this Q, can be derived and is 5-4+0-2 mev. 


Ea+ Ep ( Mev.) Q (Mev 


Figure 2(a). Energy distribution of (E,+ Ep). Figure 2 (6). Distribution of O values 


The group b of Figure 2 (a), if real, should lead to a group in position b’ in Figure 2 (8), 
and, in fact, the corresponding Q values do fail within 0-5 Mev. of this expected position. 
While the evidence is by no means as strong as in the case of the ground state, it suggests 
a broad level in the *Li nucleus at 250-5 Mev. ‘The events marked c and c’ satisfy 
equations (3) and (4) simultaneously and could possibly be due to still higher levels, but 
further evidence will be required to confirm this. It may be that higher levels will become 
so broad that there will be considerable overlapping. A continuum of O values would then 
be expected. 

The possibility that such V tracks should result from the reaction 


ST i+ n=*He+°H+465 Mev. 2 nn wee (5) 
has been considered. 

It can be shown that the geometry of all except two of the events is incompatible with the 
dynamics of reaction (5) for any neutron energy. This is to be expected since, with such a 
low proton bombarding energy and an analysed beam, the neutron production will be by 
(yn) reactions only. 

A preliminary value for the mean cross section of the reaction *Li(yn) for the *Lipy 
radiations can be calculated and is found to be (5+2) x 10-78 cm?. 
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ABSTRACTS FOR SECTION B 


An Improved Beam Divider for Fizeau Interferometers, by H. HSS and 
J. S. PRESTON. 


ABSTRACT. A transparent material, presumed to be titanium oxide, is suggested for 
replacing aluminium or silver as a coating for glass beam dividers in Fizeau interferometers, 
particularly those adapted for measuring length or displacement in the reflected fringe 
system. ‘The method of depositing the films, and their optical properties, are described; 
the coatings are very resistant chemically and mechanically. Evidence produced in 
photographs and microphotometer records, supplemented by calculations of fringe intensity 
distributions, shows that the pertormance of the new material in the interferometer is 
superior to that of aluminium and of plain glass. 


Fresnel Diffraction at a Transparent Lamina, by R. C. Faust. 


ABSTRACT. Fresnel diffraction at the edge of a transparent lamina has been studied 
and an analogy drawn between the observed phenomena and interference fringes of the 
Fizeau and Edser—Butler types. The variation of the diffraction pattern with changing 
wavelength has been employed to determine accurately both the thickness and the refractive 
index of thin transparent sheets. The refractice index was found by an immersion method, 
the accuracy of which is here discussed. The Becke line which is often observed during 
such measurements is shown to be a strongly asymmetrical diffraction pattern occurring 
when the thickness of a lamina changes over a small but finite distance. 


The Refractive Index of Electron Optics and its Connection with the Routhian 
Function, by W. GLaseEr. 


ABSTRACT. ‘The derivation of the electron-optical refractive index from Hamilton’s 
Principle is discussed. It is shown that the criticisms made by Ehrenberg and Siday of 
our method are invalid, and that out derivation made in 1933 is correct. 'To prove that our 


method is generally applicable, the isotropic refractive index of an axially symmetrical field 
is obtained. 


Stress Birefringence in Polyethylene, S. M. Crawrorp and H. Kosky. 


ABSTRACT. ~The birefringence and strain produced in polyethylene by different stresses 
have been measured at temperatures between 20° and 90° c. It has been found that the 
relation between strain and birefringence is linear up to strains of 0:25 and for greater strains 
there is still no hysteresis. This indicates that no true flow occurs in this teniperature range. 
The stress—birefringence relation, on the other hand, rapidly departs from linearity and 
shows large hysteresis, so that the nature of the birefringence effects is in direct contrast 
with that observed in rubber. An expression has been derived for the birefringence in 
terms of the orientation of rod-like crystallites embedded in an elastic matrix and this was 
found to be in qualitative agreement with the experimental results. 


The Friction and Mechanical Properties of Powders Bonded by Synthetic Resins, 
by R. C. Parker and E. J. W. WHITTAKER. 


ABSTRACT. ‘The coefficient of friction and mechanical properties have been determined 
for a series of non-metallic materials formed by bonding mineral powders with resins. 
The manner in which the coefficient of friction varies with particle size has been studied, 
and a tentative explanation has been advanced involving certain geometrical parameters. 
A correlation was observed between the coefficient of friction and the hardness of the mineral 
powders expressed on Moh’s scale. Marked correlations were also obtained between 
compression strength, cross-breaking strength and the reciprocal of the particle size for 
most of the materials; consideration is given to the various mechanical and thermo- 
dynamical factors involved. 
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Magnetic Susceptibility and Anisotropy of Mica, by J. 'T. KenpaLi and D. YEo. 


ABSTRACT. The magnetic susceptibility and anisotropy of natural muscovite and 
synthetic fluorphlogopite have been measured. These minerals contain small quantities of 
iron which cause paramagnetic and ferromagnetic effects. To a fair approximation the 
mean susceptibility is found to vary linearly with the total iron content, and the paramagnetic 
anisotropy to be proportional to the ferrous iron content. At low field strengths ferro- 
magnetic impurities (probably sub-microscopic inclusions of magnetite) cause anomalous 
results. It is shown that such inclusions are present in all the micas examined, and account 
for the apparent paramagnetism of synthetic mica. The inclusions are probably orientated 
with their magnetic axes in a preferred direction relative to the crystallographic axes of the 
mica. 


‘The Scattering of Radio Waves, by W. DIEMINGER. 


ABSTRACT. A review of the recent work on scattering draws attention to a number of 
discrepancies in the interpretation of the data. Most workers believe that the scattering 
is due to clouds in the E layer. 

The experiments described in this paper indicate that there are several types of 
scattering phenomena. It is concluded that the most common source of scatter is the 
reflection of waves from irregularities on the surface of the earth. In this case the wave 
path is from the sender to F layer—earth—F layer and hence to receiver, i.e. is analogous. 
to a 2XF reflection. The characteristics of the different types of scattering and their 
effects on the propagation of short waves are discussed. 


Pinch Effect Oscillations in a High Current Toroidal Ring Discharge, by S. W. 
Cousins and A. A.’ Ware. 


ABSTRACT. A toroidal ring discharge with currents greater than 10+ amperes has been 
studied with a rotating mirror camera. ‘The photographs show a lateral oscillation of the 
discharge, which consists of an alternate contraction and expansion of the current filament. 
The frequency of the oscillations increases with decrease in pressure or increase in voltage 
and is inversely proportional to the square root of the atomic weight of the gas. The 
oscillations are thought to be due to plasma ion waves, which are excited by the large 
‘pinch’ effect forces. 


Turbulence in Gaseous Conductors, by K. G. EmgLtus. 


ABSTRACT. It is suggested that turbulent flow of electricity may occur under certain 
conditions in an ionized gas, even if the gas is not in a magnetic field. Examples are given of 
high and low frequency disturbances which may be of this character. The difficulties in the 
way of developing a quantitative theory, and the possibility of the coupling of electro- 
magnetic and hydromechanical phenomena, are pointed out. 
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Figure 1 (6). 


Two examples of o-stars showing at least one short-range proton (marked with an arrow) and 
other more energetic particles. According to Menon et al. (1950), these are typical of disintegrations 


Figure 1 (a). 
Figure 1. 


produced by the capture of negative 7-mesons in light elements. ‘The mesons were travelling 
upwards, as shown. These events have therefore been classified as 7-meson disintegrations. 
Observers; Mrs. M. H. George and J. Evans. 
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0 > 301 
Figure 3. An example showing a 1-prong o-star observed below ground. This class of o-star is 


a by far the most frequent. In most cases the single heavy particle leaves the emulsion. 
This is one of the three cases in which the heavy particle remains in the emulsion. 


ea Observer: R. H. Creamer. 
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